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Abstract 

This report reviews the properties of Higgs bosons in the Standard Model (SM) and its various 
extensions. We give an extensive overview of the potential of the ILC operated at center-of- 
mass energies up to 1 TeV (including the 77 collider option) for the determination of the Higgs 
boson properties. This includes the measurement of the Higgs boson mass, its couplings to SM 
fermions and gauge bosons, and the determination of the spin and the CP quantum numbers of 
the Higgs. We also discuss extensions of the SM, including heavy SM-like Higgs bosons, heavy 
Higgs bosons in the framework of Supersymmetry (SUSY) and more exotic scenarios. We review 
recent theoretical developments in the field of Higgs boson physics, and the impact of Higgs boson 
physics on cosmology in several SUSY frameworks is considered. The important questions as to 
what the ILC can contribute to Higgs boson physics after the LHC, the LHC/ILC interplay and 
synergy, are addressed. The impact of the accelerator and detector performance on the precision of 
measurements are discussed in detail and we propose a strategy to optimize future analyses. Open 
questions arising for the various topics are listed, and further topics of study and corresponding 
roadmaps are suggested. 
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Chapter 1 

Introduction 



Elucidating the mechanism of clectroweak symmetry breaking (EWSB), by which the electroweak 
gauge bosons and the fermions acquire masses, wiU be one of the main tasks of future cohider 
experiments. Ultimately, all theoretical models rely on some form of the Higgs mechanism to 
explain EWSB through spontaneous symmetry breaking. In the Standard Model (SM) [1], this is 
postulated by introducing one fundamental scalar doublet field whose potential leads it to acquire 
a non-zero vacuum expectation value. Three of the real scalar components of the doublet are eaten 
by the gauge bosons, leaving one physical scalar particle referred to as Higgs boson [2]. Although 
this particle has eluded experimental observation so far, high precision measurements of electroweak 
observablcs provide indirect indications for a light Higgs boson. Updated electroweak data set an 
upper limit on the Higgs mass of about 200 GcV at 95% confidence level [3]. The preferred value 
of the Higgs boson mass as extracted from the global fits of electroweak observables is 91 GeV, 
with an experimental uncertainty of +45 and —32 GeV [3]. This result strongly depends on the 
assumption that the SM describes all physics up to high energy scales, and thus is not a proof that 
the SM Higgs boson actually exists. Rather, the results of electroweak fits serve as guidelines as 
to what mass range is to be expected. The results of experimental searches for the Higgs boson at 
LEP set a lower limit on the SM Higgs boson mass of about 114 GeV [4]. 

If a SM-like Higgs boson exists, it is expected to be seen at the Large Hadron Collider (LHC) 
at CERN [5]. To be certain the state observed is the Higgs boson, it is necessary to measure the 
couplings of this state to the W and Z gauge bosons, and to fermions such as the top and bottom 
quarks and the tau leptons. Progress along these lines will be achieved at the LHG [5-7], but the 
international linear e"'"e~ -collider (ILC) is absolutely essential for making precise measurements of 
a full set of these couplings and for establishing the nature of EWSB [8-11]. 

Many models beyond the SM, e.g., supersymmetric extensions of the SM, predict extended 
Higgs sectors containing more than one physical scalar state and therefore enriching Higgs boson 
phenomenology. A common feature of the models with such extended Higgs sectors is the existence 
of additional scalar bosons, such as charged Higgs bosons and CP-odd Higgs boson(s). Discovery 
of these extra scalar particles would be direct evidence for physics beyond the SM and be vital for 
our understanding of the structure of EWSB. Once these additional Higgs bosons are found, we 
would be able to test each new physics model by measuring details of their properties accurately at 
ILC. A 77 collider option provides further opportunity for the heavy Higgs search, where various 
properties of the heavy Higgs bosons, such as the CP-parity, can be explored. 

In this report we address the following questions: 

• What are the most important measurements that the ILC should perform in the subject of 
Higgs boson physics? 
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• What are the key measurements by which the ILC will add to what will be known already 
from the LHC? 

• For each of these measurements, what criteria for the detectors are necessary to allow the 
measurements with the required precision? 

• What are the open questions in the fields of theory, experimental analyses and detector 
developments that have to be addressed for a successful ILC program? 

• For these open questions, what is the best way to answer them in a comprehensive and 
systematic way? 

Consequently, we review the properties of Higgs bosons in various models in Sect. [21 This includes 
the possible mass ranges, as well as a general overview about production and decay modes. In 
Sect. 131 we give an extensive overview about the potential of the ILC operated at center-of-mass 
energies up to 1 TeV (including the 77 option) for the determination of the Higgs boson properties 
within the framework of the SM and its various extensions. This comprises the measurement of the 
Higgs boson mass, its couplings to SM fermions and gauge bosons, and the determination of the spin 
and the CP quantum numbers of the Higgs. The extensions of the SM that are analyzed in more 
detail are heavy SM-like Higgs bosons, heavy Higgs bosons in the framework of supersymmetry 
(SUSY) and further exotic scenarios. We review in Sect, jll recent theoretical developments in the 
field of Higgs boson physics. This includes higher-order corrections to Higgs boson masses and 
production cross sections in the SM and the Minimal Supersymmetric Standard Model (MSSM), 
the possibility of CP-violation, as well as extensions of the MSSM such as the NMSSM, gauge- 
extended models, and Fat Higgs. We also summarize work on non-supersymmetric models, such 
as Little Higgs models, models with extra dimensions, and dimension-six Higgs operators. A list of 
Higgs related computer tools (spectrum generators, decay packages, event generators) is provided 
in Sect. O The important question of what the ILC can contribute to Higgs boson physics after 
the LHC, the LHC/ILC interplay and synergy, is discussed in Sect. (HI It is emphasized that if a 
Higgs- like state is discovered at the LHC and the ILC, independent of the realization of the Higgs 
mechanism, important synergistic effects arise from the interplay of LHC and ILC [12]. Higgs boson 
physics can also have an important impact on our understanding of cosmology. This is analyzed 
for the MSSM and the NMSSM case in Sect. 13 The above listed issues and analyses impose 
certain detector requirements to achieve the necessary precision. The impact of the accelerator 
and detector performance on the precision of measurements are analyzed in detail in Sect. |H1 We 
propose a strategy to optimize future analyses. Open questions arising for the topics in the various 
chapters are listed, and further topics of study and corresponding roadmaps are suggested. 
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Chapter 2 

Production and Decays of Higgs 
Bosons 

In this section we review the experimental status of Higgs boson production and decay at the 
ILC and the photon colhder (7C). For the corresponding theoretical results and developments, see 
Sect. H 

2.1 The Higgs boson mass 
2.1.1 Limits from theory 

The Higgs boson h is yet to be discovered, and its mass rufi remains unknown. In the SM, ruh is 
a parameter which characterizes the property of the Higgs dynamics. Since oc (A is the 
quartic term in the Higgs potential, v denotes the vacuum expectation value), a light h means that 
the Higgs sector is weakly-interacting, while a heavy h corresponds to strong coupling. Although 
nifi is an unknown parameter, by imposing the requirement that the theory must be consistent up 
to a given value of the cutoff scale A, the allowed region of nih can be predicted as a function of 
A [13]. Based on this requirement, a renormalization group analysis for the coupling constant A 
gives upper and lower bounds on ruh for a given A. In the SM, for A = 10^^ GeV, the allowed 
region of rrih is evaluated as about 135 < ruh < 180 GeV, while for A = 10^ GeV it is about 
nih < 500 GeV. The results are summarized in Fig. 12.11 [14], evaluated for rrit = 175 GeV. 

Hence, as long as the SM Higgs sector is assumed, a light Higgs boson would indicate a weakly 
coupled theory with a high cut-off scale. Scenarios based on grand unified theories might correspond 
to this case. In such cases, supersymmetry would necessarily be required to reduce the problem of 
the large hierarchy between the weak scale and the scale A. On the contrary, a heavy Higgs boson 
with mh ~ several hundred GeV would imply a strongly-coupled Higgs sector with a low cutoff 
A ~ 0(1) TeV. In such a case, the Higgs sector should be considered as an effective theory of a new 
dynamics at TeV scales. Therefore, from knowing the mass of the Higgs boson, a useful hint for new 
physics beyond the SM can be obtained. Figure shows the mass range of the (lightest) Higgs 
boson under the assumption of A = 10^^ GeV in the SM [13], the minimal supersymmetric SM 
(MSSM) [15-17] (see also Sect. 14. l|) . some extended SUSY models [18] and the general two Higgs 
doublet model [19]. Due to its strong dependence on the other model parameters, the lightest Higgs 
boson mass in the MSSM and its extensions is one of the most important observables in the Higgs 
sector. 
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Figure 2.1: Summary of the uncertainties 
connected to the bounds on nih. The up- 
per sohd area indicates the sum of theoreti- 
cal uncertainties in the m/^ upper bound for 
mt = 175 GeV [14]. The upper edge cor- 
responds to Higgs masses for which the SM 
Higgs sector ceases to be meaningful at scale 
A, and the lower edge indicates a value of nih 
for which perturbation theory is certainly ex- 
pected to be reliable at scale A. The lower 
solid area represents the theoretical uncer- 
tainties in the nih lower bounds derived from 
stability requirements using nit = 175 GeV 
and Qo = 0.118. 
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Figure 2.2: Predictions of the mass 
range of a Higgs boson for various mod- 
els [20]. The upper and lower bound of 
the Higgs boson mass is derived from an 
assumption that each model is valid up 
to the cut-off scale of the theory, which 
is taken to be lO^^ GeV. For the MSSM 
case, the mass bound is obtained with- 
out reference to the cutoff scale. 



2.1.2 Experimental limits on the SM Higgs boson mass 

The mass of the SM Higgs boson is bounded from below by direct searches at LEP to be Mh > 
114.4 GeV at the 95% C.L. [4]. An upper limit can be obtained from a global fit to all electroweak 
precision data [3]. Including the new experimental nit value of mt = 172.7 GeV [21], the fit yields 
Mh = 91^32 GeV with an upper limit of Mh < 186 GeV at the 95% C.L. The fit was performed 
under the assumption that the SM including the Higgs sector gives a correct description of the 
experimental data. Therefore it does not confirm the SM Higgs sector. Higgs masses larger than 
about 200 GeV are only possible if other new physics effects would compensate the effect of the 
heavy Higgs on the electroweak precision data, see e.g. [22,23]. In most cases, some hints of these 
new physics effects are expected to be visible at either the LHC or ILC (or both). 

2.2 Standard Model Higgs Production and Decays 

In e~^e~ collisions, the SM Higgs boson is predominantly produced through the Higgs-strahlung 
process, e+e~ — > HZ, and through vector boson fusion processes e^e~ — > VeV(.{&'^&~)H . The cross 
section of the main production mechanisms as a function of the Higgs boson mass is presented in 
Fig. 12.41 for the three representative centre-of-mass energies, 350, 500 and 800 GeV. 
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Figure 2.3: Global fit to all electroweak pre- 
cision data within the SM [3] (including the 
new rrit value of = 172.7 GeV). The fit 
yields Mh = 91^32 GeV with an upper limit 
of Mh < 186 GeV at the 95 % C.L. The 
fit was performed under the assumption that 
the SM including the Higgs sector gives a 
correct description of the experimental data. 



The Higgs boson can be detected and its profile can be studied in detail over a wide mass range 
by exploiting various decay modes. Being responsible for generation of fermion and weak boson 
masses, the Higgs boson decays preferentially into the heaviest kinematically accessible final states. 
The dependence of the branching ratios of the SM Higgs boson on its mass is illustrated in Fig. 12.51 




m„(GeV) 



Figure 2.4: SM Higgs boson production cross 
sections as a functions of the Higgs boson 
mass. 



Figure 2.5: SM Higgs boson branching ratios 
as a functions of the Higgs boson mass. 



2.3 Higgs Bosons in General Two Higgs Doublet Models 

Among models with an extended Higgs sector, two Higgs doublet models are of particular interest, 
since this structure of the Higgs sector is required in MSSM [24,25]. In CP-conserving two Higgs 
doublet models, spontaneous symmetry breaking gives rise to five physical scalar states: two CP- 
even Higgs bosons, the lighter one denoted h and heavier H, one CP-odd boson A and two charged 
bosons, H^. At tree level in the MSSM, the mass spectrum is determined by tan/3 and the mass 
of CP-odd boson rriA- The mass of the light Higgs boson is constrained to be smaller than the 
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mass of the Z boson, m/j < Mz- Higher order corrections, in particular from loops involving 
third generation fermions and their superpartners, move the upper limit on mh up to about 135 
GeV [16,17]. 

In addition to supersymmetry, there are lots of motivations to consider models with two Higgs 
doublets at the electroweak scale, such as top-color [26], Little Higgs models [27], some extra 
dimension models, and the model of gauge-Higgs unification [28]. Tiny neutrino masses (e.^., 
the Zee model [29]), and extra CP violating phases which may be required for the realization of 
electroweak baryogenesis [30] also can be studied by introducing multi scalar doublets (plus singlets) 
at the electroweak scale. The analysis of phenomenology in the framework of the general two Higgs 
doublet model can make it possible to distinguish these nonstandard models. 

In the CP-conserving MSSM and more general two Higgs doublet model, the cross sections 
of the Higgs-strahlung and weak boson fusion processes involving CP-even bosons scale with the 
coupling of the appropriate Higgs to gauge bosons relative to that in the SM. These couplings are 
given by 

5feZZ,hH/vl/ ~ sin(/3 - q), 5//zz,/fTyvi/ ~ cos(/3 - a), (2.1) 

where a is the mixing angle in the CP-even sector. The set of tree level couplings between neutral 
Higgs particles and weak bosons is extended by two additional couplings, 

~ cos(/3 - a), ffHAZ ~ sin(/3 - a). (2.2) 

As a consequence, in collisions, the Higgs-strahlung and fusion processes will be complemented 
by the associated Higgs boson pair production mechanisms, e^e~ h{H)A. Neutral Higgs bosons 
can be also accessed via Yukawa processes, e~^e~ — > H(A)bb, e^e~ — > H{A)ti. The former reaction 
is of particular importance, since it allows to extend the mass reach for neutral Higgs bosons up to 
values close to the collision centre-of-mass energy. 

The MSSM exhibits a so-called decoupling limit as tha becomes large, niA ^ 200 GeV. In 
this scenario, h approaches the properties of the SM Higgs boson, while H decouples from the 
weak bosons. As a consequence the H boson production via Higgs-strahlung or weak boson fusion 
processes gets suppressed, whereas the cross section of the process e'^e~ — > HA reaches its maximal 
value. Another distinct feature of the decoupling limit is a tiny mass splitting between H and 
A^ making it practically impossible to distinguish them by mass. For a large portion of MSSM 
parameter space, decoupling is realized and the decay properties of h are close to those of the 
SM Higgs boson, while H and A decay predominantly into the heaviest kinematically accessible 
fermions, t~^t~ ^ bb or ti. However, scenarios are possible in which decay rates of heavy neutral 
Higgs bosons into supersymmetric particles become significant. 

Charged Higgs bosons can be produced in e~^e~ collisions in pairs, e^e~ H^H~. In this 
process, the mass reach for the charged Higgs boson is limited to half of centre-of-mass energy. 
Hence, it would be desirable to investigate also the rare processes of single charged Higgs boson 
production. The dominant processes of the single charged boson production are: e~^e~ — > H~^bt, 
e^e~ H^T^Vr, e~^e~" H^W~ , e~^e~ H^e~i>e- Recent theoretical calculations [31-34] 
showed that, in general, the parameter regions for which single charged Higgs boson production 
cross sections exceeds 0.1 fb are rather small beyond the pair production threshold. The main decay 
channels of the charged Higgs bosons are — > cs, — > t'^Vj- and, if kinematically allowed, 
H+ tb. 

Figure 12.61 shows contour plots of the production cross-sections in the rriyi-tan [3 plane for the 
following processes: e+e" ZH, Ah, AH, H+H', W^H^, bbA, bbH, tlA, tiH, and uVH [35]. 
In Figures 12.61 (a), (c), (d), (e) and (f), the mass reach for A, H, and H^ at the ILC is determined 
by half of ^/s. In Figure I2.6f b). if the sensitivity reaches to O.lfb, the discovery contours for 
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the bbA and bbH modes go beyond y/s/2 for large tan/3. For tan/? < 10, the ZH (Ah) mode 
is available above niA > \/s/2. In Figure ITBT cVff). cross-section contours of e+e~ — > ttA, Mi?, 
bbA, and 56/? exhibit a dependence on tan/5 for 350 GeV < ^ \/s/2. These processes include 
e"'"e~ — > AH followed by ^ or decaying into a bb or quark pair. The tan /? dependence shown 
in Figure I^BT cVff) comes from the branching ratio of the heavy Higgs bosons. Figure [Tnjc)-(e) 
also shows that the ILC will cover the region of moderate tan/3 < 10 and Ma ^ where the 

detection of the heavy Higgs bosons at LHC is expected to be difficult. If kinematically allowed, 
the heavy Higgs bosons are expected to be found in several modes at the ILC. More details can be 
found in Sect. 14.21 
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Figure 2.6: The cross-section contours: (a) for ^/s = 500 GeV and o" = 1 fb, (b) for ^/s = 500 GeV 
and a = 0.1 fb, (c) for ^ = 1.0 TeV and a = 1 fb, (d) for ^ = 1.0 TeV and a = 0.1 fb, (e) for 
= 1.5 TeV and a = 1 fb, and (f) for ^ = 1.5 TeV and a = 0.1 fb. GRACE/SUSY [36] is used 
to calculate the production cross-sections of the tree-level processes. One-loop induced production 
of W^H~ is calculated as in Refs. [31,32]. Masses of the Higgs bosons and the mixing angle of the 
neutral Higgs bosons are obtained by using FeynHiggs [37], where we assume the diagonal masses to 
be (1 TeV)2 in the stop mass matrix and maximal stop mixing. HDECAY [38] is used to calculate 
decay widths of the Higgs bosons. 



2.4 Higgs Boson Production at a 77 collider 

At a 77 collider, neutral Higgs bosons can be produced resonantly in the s-channel. Although 
Higgs bosons do not couple to photons at tree level, Higgs boson production in 77 collisions 
becomes possible due to loop-induced couplings as illustrated in Fig. 12.71 As compared to the e+e" 
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collision option, the 77 collider, operated at the same centre-of-mass energy, has higher mass reach 
for the Higgs bosons, as the whole collision energy is available for the resonant production of the 
Higgs particles in the s-channel [39-42]. In addition, we can determine the CP parity of the heavy 
Higgs boson through the process 77 — > tt by measuring the helicity of the top quark [43]. Along 
with other nice features, described later on in this report, the high-mass discovery potential of the 
77 collider makes it a particularly attractive running option of the International Linear Collider 
facility. 

Technical aspects related to the 77 collision option are highlighted in References [44-46]. 

7 

< Figure 2.7: Loop-induced cou- 
plings of Higgs bosons to photons. 
Feynman diagram involving W bo- 
son loop (left plot). Feynman di- 
agram involving charged fermion 
loop (right plot). 
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2.5 The role of polarization 

For the ILC the possibility to have a polarized e~ beam is in the baseline design [47]. A polarization 
of up to 80% is foreseen. The possibility to have polarized e"*" as well is still under discussion. Here 
in principle a polarization of ~ 60% seems to be achievable. Polarization can be helpful in various 
aspects of Higgs boson physics. 

• One of the major physics goals at the ILC is the precise analysis of all the properties of 
the Higgs particle. For a light Higgs the two major production processes, Higgs-strahlung, 
e~^e~ HZ, and WW fusion, e~^e~ Hvv, will have similar rates at \/s = 500 GeV; see 
Sect. 12.21 Beam polarization will be important for background suppression and a better sep- 
aration of the two processes. Furthermore, the determination of the general Higgs couplings 
is greatly improved when both beams are polarized [48]. 

• Searches for heavy SUSY Higgs particles can be extremely challenging for both the LHC and 
the ILC. Exploiting single Higgs-boson production in e^e~ — > vvH extends the kinematical 
reach considerably. However, in the decoupling region. Ma ^ Mz, the suppressed couplings 
of the heavy Higgs boson to SM gauge bosons lead to very small rates. This difficulty could 
be attenuated by accumulating a very high integrated luminosity, together with a further 
enhancement of the signal cross section by polarizing both beams, see also Sect. 14.21 Polar- 
ization can be decisive to discover MSSM Higgs bosons with low production cross sections [49]. 
Polarization of both the e~ and e"*" beams is important in this case. 

• Compared with the case where only the electron beam is polarized, the polarization of both 
beams leads to a gain of about one order of magnitude in the accuracy of the effective 
weak mixing angle, sin^ ^off [50]. Within the SM, this has a dramatic effect on the indirect 
determination of the Higgs-boson mass, providing a highly sensitive consistency test of the 
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model that may possibly point towards large new-physics scales [51]. Within the MSSM, 
the large increase in the precision of sin^ ^cflF will allow to obtain stringent indirect bounds on 
SUSY parameters [48]. This will constitute, in analogy to the SM case, a powerful consistency 
test of supersymmetry at the quantum level and may be crucial to constrain SUSY parameters 
that are not directly experimentally accessible [52,53]. 

In order to fully exploit the possibility of beam polarization, several questions should be ad- 
dressed in the above scenarios: 

• What are the detector requirements for these measurements? 

• What are the accelerator requirements? 

• What are the theory requirements? 
Partial answers are already available in Ref. [48] . 
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Chapter 3 

Measurements in the Higgs Sector at 
e^e~ and 77 Colliders 



In the following we review the Higgs boson sector measurements that can be performed at the ILC 
and the 7C. The description of a SM-like Higgs boson (with mass below < 140 GeV) also applies 
to the lightest MSSM Higgs boson for large parts of the MSSM parameter space. The issue of 
heavy MSSM Higgs bosons is addressed separately. 



3.1 Model Independent Determination of the Higgs-strahlung cross 
section 

One of the nice features of the linear collider is its capability to detect the Higgs boson independent 
of its decay mode. Even if the Higgs boson decays into invisible particles^, it can still be detected by 
exploiting the Higgs-strahlung production mechanism with subsequent decay of the Z into electron 
or muon pairs. The signal manifests itself as a peak in the distribution of invariant mass of the 
system recoiling against the electron or muon pair stemming from Z boson decay. This is illustrated 
in Fig. O 

The results of the dedicated study [55] on model-independent measurement of the Higgs- 
strahlung cross-section, exploiting the HZ — > X(.'^i~ channel, are presented in Table 13.11 Com- 
bining the two final states HZ Xe^e~ and HZ X^^ ^~ , the cross sections can be measured 
with statistical errors of 2.6 to 3.1% for Higgs masses from 120 to 160 GeV. 



ruH (GeV) a (fb) Z ^ e+e" a (fb) Z 



120 5.26±0.18±0.13 
140 4.38±0.18±0.11 
160 3.68±0.17±0.09 



5.35±0.21±0.13 
4.39±0.17±0.10 
3.52±0.15±0.08 



Table 3.1: The results for the Higgs- 
strahlung cross-section measurement in 
the HZ Xe+e- and HZ Xfi+i^- 
channels at ^/s = 350 GeV with an in- 
tegrated luminosity of 500 fb~^. The 
first error is statistical and the second 
is systematic. 



'^Scenarios of the Higgs boson decays into invisible particles are discussed in Sect. 13. ID. 21 
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Figure 3.1: Distribution of the invariant mass of the system recoiling against a pair leptons in the 
channel HZ Xe+e" [54] (left plot) and in the channel HZ Xfi+ii- [55] (right plot). The 
Higgs boson mass is 120 GeV. The left plot corresponds to the centre-of-mass energy of 340 and 
the right plot to 350 GeV. Assumed integrated luminosity is 500 fb~^ for both cases. 



3.2 Determination of Higgs Boson Mass 

The SM Higgs boson mass is a free parameter of the model and of great importance for the 
exploration of the EWSB mechanism. The prospects of Higgs mass measurement at a linear e~^e~ 
collider are investigated in References [56,57]. The study presented in [56] focuses on the low and 
intermediate Higgs boson masses in the range from 120 to 180 GeV. The analysis utilizes Higgs 
boson decays into 5-quark and M^-boson pairs and Z boson decays into quark and charged lepton 
pairs. Considering hadronic decays of W bosons, four distinct topologies are covered: 2jet+2£ final 
states resulting from HZ qq£^£~ , 4jet from HZ bbqq, 4jet+2£ from HZ WWi~^i~ and 
6jet from HZ WWqq. The study is performed for centre-of-mass energy of 350 GeV, assuming 
an integrated luminosity of 500 fb^^. As an example. Fig. 13.21 shows the reconstructed Higgs boson 
mass spectra in the HZ WWqq channel for mn = 150 GeV and the HZ WWi~^i~ channel 
for ruH = 180 GeV. 

From the fit of the resulting mass spectra, the Higgs boson mass can be extracted. Results of 
the analysis are summarized in Table EH 



A{mH) in MeV 



Decay mode 


120 


150 


180 


HZ qqi+r 


85 


100 




HZ bbqq 


45 


170 




HZ WW£+£- 




90 


80 


HZ WWqq 




100 


150 


Combined 


40 


65 


70 



Table 3.2: Uncertainties on the determination of mn for mn = 120, 150 and 180 GeV. For mn = 
150 GeV, the combination is done using the HZ — > WW£~^£~ and HZ — > WWqq channels. 
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Figure 3.2: Distribution of the reconstructed Higgs boson mass in the HZ WWqq channel for 
ruH = 150 GeV (left plot) and in the HZ WW£+£- channel for rriH = 180 GeV (right plot). 



The study presented in Reference [57] focuses on a heavy Higgs boson with mass greater than 
200 GeV. Exploiting 4jet+2^ final states stemming from HZ WW{ZZ)£~^£~ decay modes, the 
Higgs boson mass is extracted from a lineshape analysis of the invariant mass of the hadronic system 
stemming from Higgs decays to weak boson pairs. The relative accuracy on the Higgs boson mass 
varies from 0.11 to 0.36% for m// between 200 and 320 GeV. As discussed in Sect. 13.61 in the case 
of a heavy Higgs boson not only its mass but also its width can be determined with good accuracy 
from the mass lineshape analysis. 

3.3 Measurements of the Higgs Boson Couplings to Gauge Bosons, 
Fermions and Itself 

In the SM, the electroweak symmetry is spontaneously broken by introducing the Higgs doublet 
field. Its neutral component receives the vacuum expectation value {v). The weak gauge bosons 
then obtain their masses through the Higgs mechanism. At the same time, all quarks and charged 
leptons receive their masses from the Yukawa interactions with the Higgs field. Moreover, the Higgs 
boson (/i) itself is also given its mass {nifi) by v through the self-interaction of the Higgs boson. 
All these masses of the SM particles are expressed as multiplication of coupling constants with v. 
Therefore, there is a universality between masses and coupling constants: 

2mw V2mt V2mb \/2mr mn nAa to ^\ 
= = = = ... = —= = t) ~ 246 GeV, (3.1) 

a yt Vb Vt V2A 

where g is the weak gauge coupling, y/ is the Yukawa coupling constant to the fermion /, A is 
the self-coupling constant of h, and nii is the mass of the field i. Experimental validation of 
these relations would be a crucial step in establishment of the Higgs mechanism, the fermion mass 
generation mechanism, and the structure of the Higgs potential. This is visualized in Fig. 13.31 
where the coupling of SM particles to the Higgs boson is shown as a function of their respective 
mass. Also shown are the prospective ILC precisions as discussed in the next sections. Information 
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on the triple Higgs coupling is crucial for the reconstruction of the Higgs potential by which the 
breaking of the electroweak symmetries is induced. The accurate measurement of the self-couplings 
can also be important to discriminate models. 




10 100 
Mass (GeV) 



Figure 3.3: The couplings of the SM 
particles to the Higgs boson are shown 
as a function of their respective mass. 
Also shown are the prospective ILC pre- 
cisions as discussed in the next sections. 



3.3.1 Couplings to heavy SM gauge bosons 

Couplings to Z and W bosons can be directly accessed by measuring Higgs production rates in the 
Higgs-strahlung and W boson fusion processes. The prospects for measuring the Higgs-strahlung 
cross section at ILC was discussed in Sect. 13.11 The feasibility of measuring W fusion cross section 
at centre-of-mass energies 350 and 500 GeV is investigated in Ref. [58]. The analysis utilizes the 
H ^ bb decay mode and is based on the selection of final states characterized by two jets with the 
invariant mass compatible with mn and missing energy carried by a pair of neutrinos^. Fig. 13.41 
presents the spectrum of missing mass after applying a dedicated selection in the e~^e~ HV(,i>(. — > 
bbue^e channel. 

The selected samples contain a large contribution from the Higgs-strahlung process followed 
by Higgs decays into b quarks and Z decays into neutrinos. However, the reconstructed missing 
mass clearly distinguishes between Higgs-strahlung and VFVF-fusion processes. The missing mass 
spectrum is then fitted with two normalization factors, quantifying contributions from the Higgs- 
strahlung and fusion processes, as free parameters. Results of the fits are given in Table ESI 

Using a polarized electron beam will considerably enhance the Higgs production rate in the 
W^VK-fusion process and allow for improvement of the precision on the cross-section by a factor of 2 
or better for electron beam polarization of 80% [58]. Polarization of the positron beam could lead 
to a further improvement; see also Ref. [48]. 

3.3.2 Higgs-photon coupling 

The coupling of the Higgs boson to photons can be determined through the measurement of the 
resonant Higgs production rate in 77 collisions. Colliding photons of high energy are produced 
through the Compton back-scattering of laser light on high energy electron beams. The high 

^This analysis assumes that values for the Higgs boson mass and H bb branching ratio will be already known 
from other channels. 
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Figure 3.4: The missing mass spectrum in the e^e H{uei'e) hhvei^e channel at ^/s ■ 
(left plot) and 500 GeV (right plot). Assumed integrated luminosity is 500 fb~^. 



400 

(GeV/c^) 

350 GeV 



(GeV) 


ruH (GeV) 


AfJvi/iy-fusion 




120 


3.3% 


350 


140 


4.7% 




120 


2.8% 


500 


140 


3.7% 




160 


13.0% 



Table 3.3: Results on the accuracy on the PFiy-fusion cross-section measurements in the e^e~ — > 
HVf,i>f> — > bbveVe channel. An integrated luminosity of 500 fb~^ is assumed for each center-of-mass 
energy. 



energy photon beams resulting from this process are not monochromatic. A number of simulation 
tools have been developed, allowing for estimation of the two-photon luminosity spectrum [59,60]. 
As an example, Fig. 13.51 shows the resulting two-photon luminosity spectrum for a linear e~e~ 
collider operated at ^/s = 210 GeV as simulated using the CIRCE2 program [60]. A Higgs boson 
produced in 77 collisions can be identified via its decays into 6-quarks [61-63]. Fig. 13.61 presents 
the spectrum of the reconstructed Higgs boson mass in the channel j'j ^ H ^ bb for mn = 120 
GeV and assuming integrated photon-photon luminosity of 410 fb"^ collected at a linear e~e~ 
collider operated at center-of-mass energy 210 GeV [63]. Dedicated analyses [61-63] showed that 
the quantity T(H — > 77)BR(i7 — > bb) can be measured with an accuracy of 2.1—7.8% for niH = 
120—160 GeV, respectively. Combining this with the measurements of H ^ bb in e^e~ collisions, 
the partial decay width T{H — > 77), and therefore the photonic coupling, of the Higgs boson can 
be extracted. 

One potential caveat is that the reach analyses performed to date have assumed that the QCD 
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Figure 3.5: Luminosity spectra for the spin 
of the two cohiding photons = and Jz = 
2. Distribution is obtained with the CIRCE2 
program for a hnear e~e~ colhder operated at 
center-of-mass energy of 210 GeV. 
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Figure 3.6: Reconstructed Higgs boson mass 
in the process ^ H ^ bb assuming inte- 
grated photon-photon luminosity of 410 fb~^ 
collected at a linear e^e^ collider operated at 
center-of-mass energy 210 GeV. 



background is well-described by current theoretical estimates. A study during this Workshop 
(see [64] for details) of the QCD structure in the resolved portion of the background in Fig. I3.6l has 
demonstrated a factor of 5 uncertainty in that term, which raises the overall background uncertainty 
by a factor of 2. The origin of this large uncertainty is that the gluonic component of the photon 
is nearly unconstrained at these large energies. The gluon contribution will be well-measured by 
an in situ determination once the collider is running. Hence, at worst, a slightly larger luminosity 
may be required to obtain the high precision measurements of the photon coupling to the Higgs. 



3.3.3 Higgs couplings to SM fermions 

Information on Higgs couplings to SM fermions can be extracted through the measurements of the 
Higgs boson decay branching fractions. Two approaches for the measurements of Higgs branching 
fractions can be used: 

• Extract topological cross sections, a{e~^e~ i?z^ei^e)BR(i? X) or a{e~^e~ — > HZ)IiIl{H — > 
X) from the event rate measurement and then divide by the total Higgs-strahlung or WW- 
fusion cross-section as obtained from corresponding studies. 

• Select unbiased inclusive samples of Higgs decays, e.g., events in the recoil peak in the HZ — > 
X£~^£~ channel, and then determine from this sample the fractions of events corresponding 
to a given decay mode H ^ X. Since in this direct method binomial (or generally multi- 
nomial) statistics is employed, the errors on the branching ratios are reduced by a factor of 
^1-BR{H ^ X). 

Examples of the first approach are given by a number of analyses examining the potential of 
the linear collider for the measurement of the Higgs boson branching fractions into vector bosons 
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and hadrons [65-68]. These analyses exploit either lyW^-fusion or the Higgs-strahlung process 
with various possible Z boson decays. The approach based on binomial statistics is described in 
Reference [69]. Although relying only on Higgs-strahlung events with Z — > ^ this method 
yields errors very similar to the analyses based on the first approach, exploiting a variety of Z 
boson decays following e^e~ HZ. Results on Higgs branching ratio measurements at ^/s = 800 
GeV and below are given in Table ITU 



ruH (GeV) 


120 


140 


160 


180 


200 


220 


Decay 




Relat 


ive precision (%) 






bb 


2.4(a)/1.9(e) 


2.6(a) 


6.5(a) 


12.0(d) 


17.0(d) 


28.0(d) 


cc 


8.3(a)/8.1(e) 


19.0(a) 












5.0(a)/7.1(e) 


8.0(a) 










99 
WW 


5.5(a) /4.8(e) 
5.1(a)/3.6(e) 


14.0(a) 
2.5(a) 


2.1(a) 








zz 






16.9(a) 
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Zl 


23.0(b)/35.0(e) 


27.0(c) 











Table 3.4: Precision on the Higgs boson branching ratio measurements [65-69]: (a) for 500 fb ^ 
at 350 GeV; (b) for 500 fb"^ at 500 GeV; (c) for 1 ab'^ at 500 GeV; (d) for 1 ab'^ at 800 GeV; 
(e) as for (a) but using binomial statistics method. 

Recently, a dedicated analysis based on detailed flavor-tagging tools has been performed [70] 
to evaluate the linear collider potential for the measurement of hadronic decays of the light Higgs 
boson, niH = 120 GeV. The study is based on selection of an inclusive sample of hadronic Higgs 
decays using the Higgs-strahlung process and exploiting all possible decay modes of the Z boson. 
Each selected event is assigned a quantified probability to contain b- or c-jets, referred to as 6- 
and c-tag variables. The branching ratios H ^ bb, H ^ cc and H ^ gg are determined from 
the fit of a two-dimensional distribution of 6-tag versus c- tag variables (Fig. 18. 5() with three free 
normalization parameters quantifying the fractions of — > 66, -ff — > cc and H ^ gg events in the 
final selected samples. Results of this study are presented in Table 031 



Decay 



z e+t 



Z 



vv 



Z ^ qq Combined 



Relative precision (%) 



H ^bb 
H ^cc 
H^gg 



3.0 
33.0 
18.5 



2.1 
20.5 
12.3 



1.5 
17.5 
14.4 



1.1 

12.1 
8.3 



Table 3.5: Relative precision on 
the measurement of the hadronic 
branching ratios of Higgs for mn 
= 120 GeV. Analysis is performed 
for ^/s = 350 GeV, assuming an 
integrated luminosity of 500 fb~^. 



One can improve the precision of branching ratio measurements at higher e~''e~ collision energies. 
This is done taking advantage of the increase in the Higgs production rate with energy due to a 
log ^/s dependence of M^14^-fusion cross-section. Moreover, the specific luminosity is expected to 
increase with increasing center-of-mass energy. Both these factors will enhance the signal statistics, 
allowing not only to improve the precision on Higgs branching ratios for the major decay modes, 
but also to measure rare Higgs decays, e.g., H — > [71]. An example of this signal is shown in 

Fig. Em 
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Figure 3.7: Invariant mass of the two 
muons in the selected sample of H ^ 
decays following VFVF-fusion pro- 
cess. Higgs boson mass is 120 GeV. 
Center-of-mass energy is 800 GeV. An 
integrated luminosity of 1 ab~^ is as- 
sumed. 



3.3.4 Top Yukawa Couplings of the Higgs Boson 

The top quark is exceptionally heavy as compared to the other quarks, with its mass being mt ~ 
v/V2. The size of the top Yukawa coupling is nearly 1. The fact that the mass of the top quark 
is the same size as the scale of electroweak symmtry breaking may suggest that the physics of 
the top quark may be strongly related to the dynamics of the electroweak symmetry breaking. 
Therefore, the experimental determination of the top Yukawa coupling is expected to be a very 
sensitive probe of possible new physics that might be responsible for generating the top quark mass. 
For example, alternative models of electroweak symmetry breaking with new strong interactions, 
such as Technicolor and Topcolor models, substantially modify the top quark interaction with the 
Higgs sector and give rise to new signals that could be studied at the ILC. 

If the Higgs boson mass is below 200 GeV, the top Yukawa coupling could be measured in 
associated ttH production, with an error of the order of 15% at the LHC [7] (with "mild" theory 
assumptions). At the ILC [72] in the SM the production cross section becomes maximal around 
^/s = 700 GeV for nih = 120 GeV. The expected accuracy for the top Yukawa coupling in the SM is 
4.2% at ^ = 700 GeV with / £ = 500 fb'^ An error of 5.5% can be achieved for = 800 GeV, 
/£ = 1000 fb"^ and rrih = 120 GeV [8]. However, recent studies show that for this Higgs mass and 
integrated luminosity even at ^/s = 500 GeV a measurement at the 10% level might be possible [73]. 
For a combination of LHC and the 500 GeV ILC, see Sect. 16.2.11 



3.3.5 Higgs Boson Self-Couplings 

Information about the structure of the Higgs potential will be obtained by measuring the triple 
Higgs boson coupling. Accurate information on this coupling is important to discriminate models 
beyond the SM. Even when the other Higgs boson couplings such as hZZ, hWW and hff are in 
good agreement with the SM predictions, the Higgs self-couplings can significantly deviate from 
the SM prediction due to non-decoupling quantum effects of heavy particles. In the two Higgs 
doublet model, radiative corrections of O(100)% on the hhh coupling are possible if the additional 
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heavy Higgs bosons show a non-decouphng property [74]. Such a non-decoupHng phenomenon is 
known to be required for a scenario of electroweak baryogenesis [30]. Therefore, the measurement 
of the triple Higgs boson couphng may be used to test such a model of the baryon asymmetry of 
the universe [75]. 

The hhh coupling can be measured via the double Higgs boson production processes e^e~ — > 
Zhh and also e~^e~ — > uVhh [76]. In these processes, in addition to the diagram involving the 
hhh coupling, there are double Higgs radiation diagrams which depend only on the hZZ or the 
hWW coupling. Hence, the production cross section is not simply proportional to the square of the 
self-coupling constant. Systematic studies on the self-coupling measurement have been performed 
for various values of ruh and y/s [77-79]. As seen in Figure (left) . the hhZ mode is sensitive to 
the deviation 6X3 of the hhh coupling from the SM value in the low values of the invariant mass 
Mhh of the hh system [78] . A cut on Mhh is therefore useful to improve the sensitivity. For ^/s > 
1 TeV, the sensitivity is further improved by using initial electron (and positron) polarization. The 
production cross section of the WW fusion process can be increased by up to a factor of 2 (4) in 
principle with electron (electron and positron) polarization. 

For y/s = 500 GeV, e~^e~ Zhh has much larger cross section than the e~^e~ vvhh process. 
For a higher energy, the production cross section of the latter process is enhanced because of its 
t-channel nature, whereas e+e" — > Zhh decreases as 1/s. For ^/s > 1 TeV, the WW fusion process 
e^e~ vvhh becomes important to determine the hhh coupling. The statistical errors in the hhh 
measurement are shown in Figure 13.81 (right) as a function of the Higgs boson mass from 100 to 
200 GeV, assuming an integrated luminosity of 1 ab~^ [78]. Dashed (dotted, sold) lines are the 
results for the hhZ mode {hhvv mode, both modes combined). At ^/s = 500 GeV, where the hhZ 
mode is dominant, the hhh coupling may be measured with about 20% accuracy for relatively light 
Higgs bosons with ruh ^ 150 GeV. For ^/s ^ 1 TeV, where the hhvv mode is dominant, higher 
sensitivity ((5A3/A3 < 10%) can be expected by using the invariant mass cut M^h < 600 GeV for 
the hhZ mode and a 100% polarized electron beam for the /i/iz^P mode. 
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Figure 3.8: (Left figure) The hh invariant mass dependence of the hhZ mode for several values 
of (5 A3. (Right figure) The A3 measurement sensitivity; hhZ (dashed line), hhuu (dotted line) and 
combined results (solid line) [77]. 
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The full calculation of the SM electroweak radiative corrections to the double Higgs-strahlung 
process e"'"e~ — > Zhh has been presented in Ref. [80]. 



3.4 Higgs Boson Width Measurement 

At Higgs masses below approximately 2My/, the detector resolution on the Higgs boson mass for 
the main decay channels, H ^ bb and H — > WW, is significantly larger than the natural width 
of the Higgs boson. Th- As a consequence, no direct measurement of the Higgs width is possible 
through the analysis of the reconstructed Higgs boson mass lineshape. However, using the relation 



T(H 



X) 



(3.2) 



with X = WW or 77, the width can be determined indirectly. For example, the partial width 
T[H 77) can be measured with the 77 collider whereas the branching fraction BR(i? — > 77) 
is accessible in e~^e~ collisions from the diphoton invariant mass spectrum. The same procedure 
can be performed for the VFH^-fusion process and the H WW decay mode. An accuracy in the 
determination of Th between 4% and 15% can be achieved [58] for mn up to 160 GeV. 

The natural width of heavy Higgs boson with mass above 2Mz becomes as large as few GeV 
and can be directly measured from the mass lineshape analysis [57]. Fig. I3.9l presents the spectra of 
reconstructed Higgs mass in the HZ WWl~^i~ and HZ — > ZZi~^£~ channels with subsequent 
hadronic decays of weak bosons for niH = 200 GeV. From the fit of the spectra with convolution 
of the detector resolution functions and Breit-Wigner function, the Higgs width can be determined 
with the relative precision of about 30% for mn = 200—320 GeV. 
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Figure 3.9: Reconstructed Higgs boson 
mass lineshape in the HZ WW£^£^ 
and HZ ZZi^i~ channels with sub- 
sequent hadronic decays of weak bosons 
for niH = 200 GeV. Selected sample cor- 
responds to 500 fb~^ of luminosity col- 
lected at center-of-mass energy of 500 
GeV. 



3.5 Spin and CP Quantum Numbers of Higgs Boson 

Full establishment of the Higgs mechanism implies measurements of Higgs spin and CP quantum 
numbers. The SM Higgs boson (like the CP-even MSSM Higgs bosons) has spin and CP-quantum 
numbers J^'-^ = 0^+. Information on CP-quantum numbers can be extracted by studying angular 
dependences of the Higgs-strahlung process. The analysis for the SM Higgs boson can also be 
apphed to a SM-hke MSSM Higgs boson. 
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The following quantities are sensitive to the CP quantum numbers of the Higgs: the polar 
production angle of the Higgs boson, 9; the angle between the Z momentum vector and the mo- 
mentum vector of the fermion stemming from Z decay, 0*\ and the angle between the production 
and decay planes of the Z boson, cj)* . The latter two angles are defined in the Z boson rest frame. 
The definitions of these angles are illustrated in Fig. 13.101 

Studies presented in Refs. [81,82] investigated the general case of a mixed scalar state $ con- 
taining both CP-even and CP-odd components. In the presence of a CP-odd admixture in the 
scalar field, the amplitude for the Higgs-strahlung process reads: 

Mz^ = MzH + iV-MzA, (3.3) 

where M.zH is the CP-even SM-like s-wave amplitude, and MzA is a new CP-odd p-wave amplitude 
contributing with strength r/ to the Higgs-strahlung process. The presence of the latter modifies the 
three-fold differential cross section d^a/d{cos 6)d{cos 9*)d{cos cj)*) with respect to the SM prediction. 
The parameter 77 can be determined using the method of optimal observables [81] . Neglecting terms 
quadratic in 77, the optimal observable is chosen to be 

^^rnHMi^^ (3.4) 

The choice of the observable is very close to optimal for small rj. At the next step, a "gauge" curve 
is produced, which quantifies the dependence of the mean value of the optimal observable and 
its error on the parameter rj. This is done using reference distributions of the optimal observable 
corresponding to different values of rj. Reference distributions are obtained for the selected samples 
in the HZ Xi'^i~ channel. The parameter rj can then be read from this "gauge" curve including 
the confidence band. The analyzing power of this method can be enhanced by taking into account 
the dependence of the Higgs-strahlung cross section on rj. Technically, the parameter rj is calculated 
by minimizing the x^j 

where {Omeas) and N^eas are the measured mean value of the optimal observable and the number 
of observed signal events in the hypothetical data sample, E{{0)){r]) and Nexp{rj) are the expected 
mean value of the optimal observable and expected number of signal events, and dgxpiv) is the 
expected error on E{{0)){r]). E{{0)){r]), Nf^xpij]) and a^xpiji) are determined from the "gauge" 
curves shown in Fig. 13.111 For small values of t], t] < 0.1, this parameter can be measured with a 
precision of about 0.03. 

An alternative method [82] consists of generating three-dimensional distributions in cos 6, cos 9* 
and cos 0* for various values of ry. These distributions are generated for the selected HZ — > Xi'^i~ 
signal samples corresponding to 500 fb~^ of data collected at ^/s = 350 GeV. In the next step, the 
likelihood 

Naata{i,j,k) ^^^^^ 

c= n t (■ -i-M (^-^^ 

is maximized, where Ndataihj, k) is the number of events of the hypothetical data sample and fUjk 
is the expected number in the ijk-th bin. /Xj^fc is calculated assuming a linear combination of the 
number of events of three Monte Carlo samples, corresponding to the production of scalar Higgs 
(MC{ZH)), pseudoscalar Higgs {MC{ZA)) and events for the interference term (MC(/A^)): 

fiijk = aMC{ZH)ijk + mCiIN),jk + -fMCiZA)ijk, (3.7) 
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Figure 3.10: Definitions of the Z production 
and decay angles 6*, 0* and 0* in the Higgs- 
strahlung process. 



Figure 3.11: Gauge curve reflecting de- 
pendence of the mean value of the optimal 
observable and the Higgs-strahlung cross- 
section on the parameter r]. Bands indicate 
Icj confidence level intervals as obtained for 
a data sample corresponding to 500 fb~^ at 
^ = 350 GeV. 



where a, 7 and /3 quantify contributions from the production of the scalar state, pseudoscalar state 
and their interference, respectively. The likelihood is then maximized with respect to a, /5 and 7. 
The quantity 77 and its experimental error can then be determined by parameterizing the fraction 
of pseudoscalar component 7 as a function of r\. It has been shown that with this method, the 
parameter r\ can be measured with an accuracy of 0.015 for \vj\ < 0.1. 

The Higgs boson parity can be established also by analyzing the spin correlation effects in the 
Higgs boson decay into fermion pairs [83]. The spin dependence of the decay probability is given 



where s and s denote spin vectors of fermion / and anti-fermion / in their respective rest frames 
(with the z axis oriented in the fermion flight direction). The positive sign in the transverse spin 
correlation term holds for the scalar particle (H) and the negative sign for the pseudoscalar (A). 
Although the dominant fermionic decay mode of the Higgs boson is into b quark pairs, it can hardly 
be utilized for measurement of the Higgs parity, as hadronization of quarks will inevitably dilute 
the spin information. The most promising channel is the Higgs decay into tau leptons, with their 
subsequent decay into p meson and neutrino [84]. The angle between the planes spanned by 
the momentum vectors of the products of the vt'^tt" and p~ — > 7r~7r^ decays in the p~^p~ rest 

frame can be used to distinguish between the scalar and pseudoscalar states. Recently, a dedicated 
study [85] has been done which investigated the feasibility of measuring Higgs parity at ILC in 
the H T^T^ p'^i7t-p~Vt tt'^ tt^ i/t-it' 1'^ channel. The analysis is performed at the level of 
detailed detector simulation and includes accurate estimates of the SM background. Higgs boson 
decays into tau leptons are selected using the Higgs-strahlung process with subsequent decays of 
the Z boson into electron and muon pairs and quarks. Even in the presence of SM backgrounds, 
the acoplanarity angle clearly discriminates between CP-even and CP-odd states, as illustrated 



by 



r{H,A^ff)r.l-s,s,±s^s^, 



(3.8) 
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Figure 3.12: The distribution of the acopla- 
narity angle ^* as defined for the H — > 
r^T~ p^u^-p^Ur — > TT^ Ut-tt^ t:^ v-r decay 
chain (see text) for the scalar state (error bars) 
and pseudoscalar state (histogram) in the se- 
lected sample of HZ — > r^r^qq events. Se- 
lected sample corresponds to 1 ab~^ of data 
collected at ^/s = 350 GeV. Higgs boson mass 
is 120 GeV. 
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Figure 3.13: The acoplanarity angle ^* in the 

decay chain. Gaussian smearing of <I>* and 
Higgs boson momenta are applied according 
to anticipated resolutions of the ILC detector. 
Thick line corresponds to a pure scalar state, 
while the thin line to a mixed one. 



in Fig. rrni 

With 1 ab^^ of data collected at y/s = 350 GeV, scalar and pseudoscalar states can be distin- 
guished at ~ Act level. Furthermore, the method using the acoplanarity angle <I>* can be applied in 
the general case of the mixed Higgs state, comprising both CP-even and CP-odd components.^ The 
presence of a CP-odd admixture results in a phase shift in the distribution of the ^* angle as shown 
in Fig. 13.131 From the measurement of this shift the mixing angle between scalar and pseudoscalar 
components can be determined [86] with an accuracy of about 0.2 rad if one assumes 1 ab~^ of 
data collected at y/s = 350 GeV. 

Recently, a new method to measure CP violation in a light Higgs boson decay from r-spin cor- 
relations has been proposed [87]. The method has been devised to be insensitive to beamstrahlung 
and exploits the information encoded in a unit vector a, the r polarization analyzer. The lat- 
ter can be computed [88] from the measured decay products for the three hadronic decay modes, 
r — s- TTv, pu^aiv, which constitute 55% of all decay channels of the r lepton. For example, in the 
case of T ^ vrz^ the polarization analyzer is represented by the pion direction. For t ^ pv and 

a'=M(2{q-p,)q'-{q-qy^), (3.9) 

where is a normalization factor, Pu = Pt — Pp is the neutrino momentum and q = Pj^± — p^o is 
the difference between the momenta of the charged pion and neutral pion. Once the vector a is 

^This scenario requires redefinition of acoplanarity angle resulting in coverage of [0, 2tv] range. For details see 
Reference [86]. 
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reconstructed, the information on the CP phase in the light Higgs state can be accessed through 
the distribution of the azimuthal angle 

dN/d(A(j)) ~ 1 cos(A0 - 2iP), (3.10) 

16 

where ip is the CP phase of the Higgs boson and A(j) is the difference in azimuthal angles of the 
vectors in the H t~^t rest frame. The reconstruction of A(j) is possible in the Higgs- 
strahlung process. It requires information on the visible decay products of the tau leptons, Z 
boson momentum, primary vertex position and the Higgs boson mass. The entire algorithm of 
reconstruction of A(j) is described in detail in Ref. [87]. Visible decay products of tau leptons are 
directly measured using only tracking in the case of the t ^ ttv decay mode or using calorimeter 
and tracking information in the case of r — > pi'. The momentum of the Z boson can be determined 
from the Z e~ , fx'^ fi~ , qq decays. The primary vertex position in rep plane is constrained 

by the beam spot, whereas the z coordinate is determined from charged decay products of the Z 
boson. Finally, it is assumed that the Higgs mass will be determined with high precision prior to 
this analysis. Distributions of A(j) when only beamstrahlung is taken into account are shown in 
Figure 13.141 The sensitivity of the tp measurement is expressed as 




Figure 3.14: The distributions of A^: a) at the generation level; b), c) and d) after dedicated 
reconstruction [87] for the three channels ttvkv, iriypiy and pi'pu at center-of-mass energy of 350 
GeV. Beamstrahlung effects only are taken into account. The histograms are normalized to the 
number of generated events and multiplied by the number of bins. The full lines correspond to a 
pure scalar Higgs state {ip = 0), the dotted to a mixed state with ip = vr/S. 

= 1/a^VN, (3.11) 
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where is the error on ip expected from a maximum hkehhood fit of the distribution for a 
sample of N events. Defined in such a way, the sensitivity measures the information per event on ip 
contained in the distribution ()3.1U|) . Anticipated sensitivities in various decay channels for different 
center-of-mass energies are given in Table EEl 



y/s (GeV) Sensitivity (5*^) 





TTVTTV 




pvpv 


230 


0.92 


0.88 


0.83 


350 


0.91 


0.73 


0.66 


500 


0.88 


0.64 


0.55 



Table 3.6: Sensitivity for the determina- 
tion of CP phase V of the Higgs boson for 
the three different channels ■kvkv^ irvpu and 
pvpu at center-of-mass energies 230, 350 and 
500 GeV. 



To evaluate detector effects, a semi-realistic simulation oi H ^ tt ^ TTi/iriy was performed. For 
the charged tracks an independent Gaussian smearing is performed on the five parameters: 9, 4>, 
1/pt and the two components of impact parameter resolution. Assumed widths of the Gaussians 
are: 

a{e) = a{(p) = 0.1 mrad 
a{l/pT) = 5 X 10-5 GeV"^ 

^m. 

The jet energy resolution is assumed to be 0.3/ y^Ej^. The reconstructed distributions of Acp after 
inclusion of the detector effects are shown in Figure 13.151 for the decays of Z both into /x^/i^ and 
qq. The corresponding sensitivities are given in Table 1X71 




Figure 3.15: The distributions of A(j) in the t^vku channel for the a) Z — > /i^/i^ and h) Z ^ qq 
decays at center-of-mass energy 350 GeV. All experimental effects are included. The full lines 
correspond to a pure scalar Higgs state (■0 = 0), the dotted to a mixed state with ip = vr/S. 



^/i (GeV) 


Sensitivity {S^) 






Z —>■ p~^p~ Z ^ qq 


Table 3.7: Sensitivity of the determination 


230 


0.69 0.71 


of CP phase V of the Higgs boson for the 


350 


0.60 0.61 


TivTiv channel when all experimental effects 


500 


0.58 0.58 


are included. 



The analyses of other decay modes, including the study of detector effects and vr reconstruction, 
have yet to be performed; nevertheless it appears that a reasonable goal for the measurement of 
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the phase V' should be to use ah the above mentioned channels and reach a sensitivity better than 
0.5, i.e., < 0.6■^■/^/I%^. 

Finally, Higgs boson parity can determined in 77 collisions of linearly polarized photon beams [83] . 
The matrix elements for the production of scalar and pseudoscalar Higgs bosons in 77 collisions 
can be expressed as: 

7W(77 ^ ^r[0++]) ~ el • 6*2, 7W(77 ^ ^[0^+]) ~ e*i x 62 ■ (3.12) 

where ei and 62 denote polarization vectors of colliding photons and is the momentum vector 
of one of the colliding photons. Hence, the parity of the produced Higgs state can be extracted 
by measuring the dependence of the production rate on the relative orientation of the polarization 
vectors of the linearly polarized photon beams. 

The measurement of Higgs spin can be performed by analyzing the energy dependence of the 
Higgs boson production cross section just above the kinematic threshold [89,90]. For a spin zero 
particle the rise of the cross section is expected to be ~ /?, where /? is the velocity of the boson in 
the center-of-mass system. For a spin-one particle the rise is ~ 0^ and for spin- two like ~ (3^. With 




Figure 3.16: The cross section of HZ 
qq£~^i~ just above the threshold assuming 
TTiH = 120 GeV. The dots correspond to a 
measurement and the curves are predictions 
for several spins. 



210 220 230 240 250 

Vs (GeV) 

a very small luminosity of about 20 fb~^ per energy point the scalar nature of the Higgs boson can 
be established and other spin hypotheses are disfavored, as shown in Fig. 13.161 It should be noted, 
however, that there are particular scenarios for spin-1 and -2 particles which show a threshold 
behavior similar in shape to the spin-0 case. This can be disentangled using angular information 
in addition. 

3.6 Heavy SM-like Higgs boson 

Most analyses for SM-like Higgs bosons have been performed for Higgs masses below ~ 150 GeV. 
However, higher masses are still allowed by the electroweak fits; see Sect. 12.1.21 Alternatively, 
new physics could compensate the effects of a heavier Higgs boson in the electroweak precision 
observables [22,23]. For heavy SM-like Higgs bosons, Mh > 150 GeV, the main decay channels 
are into WW and ZZ. Most analysis performed for lower masses can therefore not be taken over 
to this case. 
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3.6.1 Heavy Higgs Production 



In addition to the treatment of different decay modes, a very heavy SM-hke Higgs demands special 
theoretical treatment, because its width becomes non- negligible. For large M^, the width of the 
Higgs is dominated by decays into longitudinal gauge bosons, and grows as M^, 

where V = W, Z , and Qv is the electric charge of V. This behavior may be simply understood 
from the fact that the would-be Goldstone bosons are part of the Higgs doublet, and couple to the 
Higgs proportionally to the Higgs quartic A = Mh/v (at tree level). It is worth noting that the 
partial width into top quarks grows only as a single power of Mh, which explains why the decay 
into tops never dominates. 
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Figure 3.17: Cross sections for 
ZW~^W~ production function of 
for a variety of Higgs masses. The 
solid lines of each curve include the 
appropriate SM Higgs width effects, 
whereas the dotted lines indicate the 
nominal threshold for HZ production. 



Once the width is large, the Higgs production signal and the SM background processes can 
interfere in an important way. As a result, the Higgs contribution in a given channel may be 
important below the putative threshold energy for on-shell Higgs production, and the behavior 
of the signal slightly above threshold may be modified in an important way by the width [23]. 
This is illustrated in Figure l!^. 171 for the reaction e^e~ — > W~^W~ Z, including the Higgs resonant 
production via e'^e" HZ followed hy H ^ W~^W~, as well as all of the SM non-resonant 
production processes. For heavy Higgses, typically the associated production with a Z boson, 
e'^e~ — > HZ, and the I^-fusion process e~^e~ — > Ve^^eH have roughly similar rates. The exception 
is when the Higgs is heavy enough that on-shell H and Z cannot be produced, in which case the 
M^- fusion rate may dominate (see Figure [TIJ). 

3.6.2 Existing analyses 

A heavy SM-like Higgs boson decays predominantly into WW and ZZ. For a Higgs boson with 
mass 2Mz < mn < 2mt, these final states constitute almost 100% of all decay modes. This scenario 
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has been investigated separately [57]. In this study Higgs events are selected exploiting the Higgs- 
strahlung process with Z — > i'^i' . Only hadronic decays of W and Z bosons are considered. 
Separation of WW and ZZ decay modes is done on the basis of mass information. The hadronic 
system assigned to the decay of the Higgs is clustered into four jets. For each possible jet-pairing a 
kinematical fit is performed imposing four-momentum conservation, taking into account the center- 
of-mass energy {^fs = 500 GeV) and the measured momenta of the two leptons stemming from Z 
decay. An additional constraint is imposed forcing the invariant masses of the two dijet systems to 
be equal, leading to a five constraint (5C) fitting procedure. With this procedure, the WW and 
ZZ decay modes can be efficiently separated by the mass obtained with the 5C kinematical fit. 
This is illustrated in Fig. 1,3.181 
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Figure 3.18: The distribution of the 
reconstructed vector boson mass in 
the sample of selected Higgs-strahlung 
events with H WW or H ^ ZZ and 
Z — > Distribution is obtained 

with 5C kinematic fit as described in 
text. The WW and ZZ decays of the 
Higgs boson are clearly separated. The 
Higgs mass is 200 GeV. The sample cor- 
responds to integrated luminosity of 500 
fb-i collected at ^fs = 500 GeV. 
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With clear separation of the WW and ZZ decay channels, the measurement of the corresponding 
branching ratios of heavy Higgs can be performed. Results are given in Table 13.81 



niH (GeV) 


H WW 


H ^ ZZ 




Relative precision (%) 


200 


3.5 


9.9 


240 


5.0 


10.8 


280 


7.7 


16.2 


320 


8.6 


17.3 



Table 3.8: Precision on H ^ WW 
and H ZZ branching ratio measure- 
ments for a heavy Higgs boson [57]. Re- 
sults are obtained assuming 500 fb^^ of 
luminosity collected at ^/s = 500 GeV. 



The prospect of measuring Higgs branching ratios with 1 ab~ of e~^e~ collision data collected 
at 1 TeV center-of-mass energy has been investigated in Ref. [91]. The analysis is based on the 
measurement of topological cross sections c'"vKVF-fusionBR(-ff X) with X = bb,WW, ggj^yy and 
makes use of results for the Higgs branching ratio measurements into bb and WW at lower center-of- 
mass energies. Electron and positron beam polarization of 80% and 50%, respectively, is assumed. 
The study makes use of the following relations: 

BR{H -^X) = [aww-insionBR{H ^ X)) (^aww -insionBR{H ^ WW)y^BR* [H ^ WW) 

= (aww-iusion^HH ^ X)){aww-insionB^{H ^bb)y^BR*{H ^bb), (3.14) 

where aww -insion^^iH X) is the topological cross section measured at 1 TeV and RR*{H — > 
W^VF) and RR*{H bb) are branching ratios measured at 350 GeV center-of-mass energies. Using 
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these relations, a least-squares fit is performed to obtain measurement errors for BR(// bb), 
BR(i? WW), BR(i? gg) and BR(i/ — > 77). This procedure provide significant improvement 
(by a factor of two or more) for decay modes with small branching fractions, such as BR(i/ — > bb) 
for 160 GeV < m// < 200 GeV, m.{H WW) for 120 GeV < < 140 GeV, and BR{H gg) 
and BR(ff — > 77) for all Higgs masses. More details can be found in Ref. [91]. 

Another study [92] , making use of the large number of Higgs bosons produced in the WW fusion 
channel, obtained the following results. For 1 ab~^of data at ^/s = 800 GeV a So" sensitivity to the 
bottom Yukawa coupling is achievable for Mh < 210 GeV. A measurement of the branching ratio 
BR(ff — > bb) is possible with the accuracies given in Tab. 13.91 



niH (GeV) 


H ^bb 




Relative precision (%) 


160 


6.5 


180 


12.0 


200 


17.0 


220 


28.0 



Table 3.9: Precision on H ^ bb branch- 
ing ratio measurements for a heavy 
Higgs boson [93]. Results are obtained 
assuming 1 ab~^ of luminosity collected 
at ^/i = 800 GeV. 



Another analysis concerns the top Yukawa coupling measurement for Higgs bosons heavier 
than 2mt. In this case the top Yukawa coupling could be determined from the Higgs decay into tt 
pairs. This decay can be studied at the ILC using the Higgs resonant contribution to the process 
W'^W" tt. While this electroweak tt production process cannot be observed at the LHC due to 
the huge QCD background of tt production by gluon fusion, it can be measured via a high energy 
e~^e~ linear collider. A complete simulation study of the W~^W~ H ^ tt process from e~^e~ 
collisions has been performed including realistic backgrounds and experimental effects [94]. The 
helicity cross sections for the e~^e~ ttvu signal process, including ISR and beamstrahlung effects, 
are given in Table IXTUl for a collider energy of ^/s = 1 TeV and different Higgs boson masses. 



100 


0.11 


0.22 


0.21 


0.65 


400 


1.89 


0.34 


0.33 


4.45 


600 


0.89 


0.22 


0.21 


2.21 


800 


0.33 


0.22 


0.21 


1.09 



niH ajLL) = a{RR) a{LR) a{RL) Total Table 3.10: Hehcity cross sections (in fb) 

for the e+e" — > tiui) process, at -y/s = 1 
TeV and for different Higgs boson masses 
(in GeV). a{X,\') denotes the cross section 
for production of a top quark with helicity 
A and an antitop with helicity A'. CP in- 
variance implies a{RR) = a{LL). 

The relevant signal is a tt pair with missing transverse momentum carried by the neutrinos or 
from beam electrons not observed by the detector. The largest background is direct tt production. 
It can be suppressed by requiring missing mass greater than some minimal value, and further 
by choosing the jet association that best reconstructs the t and W masses. Another significant 
background is e~^e~tt, which can be reduced by requiring the missing transverse energy to be 
greater than some value such as 50 GeV as used in Ref. [94]. The background from 7W can be 
mostly eliminated with an appropriate cut on the of the ti pair [95]. The expected number of 
reconstructed 6-jet events as a function of the visible mass at ^/s = 1 TeV with C = 1 ah~^ is 
shown in Fig. 13.191 (left), and the expected relative precision in the top Yukawa measurement is 
shown in Fig. 13.191 (right) [94]. The QCD corrections of 0{as) to the VlVl tt scattering have 
been studied in Ref. [96], where Vl represents the longitudinally polarized vector boson V = W^ 
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Figure 3.19: [Left figure] Expected number of reconstructed 6-jet events as a function of the visible 
mass, at ^/s = 1 TeV and with C = 1 ab~^ after all cuts, for a Higgs boson of 400 GeV (dots) and 
for the backgrounds (dashed). The expectations (including background) for a Higgs of 100 GeV 
(solid) are also shown for comparison. [Right figure] Expected relative precision in the top Yukawa 
coupling mesurement as a function of the Higgs boson mass, from 6-jet events. 

or Z. It has been found that corrections to the cross sections can be as large as 30% and must be 
accounted for in any precision measurement of VV — > tt. 

The second question about heavier Higgs bosons is whether the Higgs line-shape parameters 
(mass, decay width, Higgs-strahlung production cross section) can be measured. A complete study 
of the mass range 200 GeV < niH < 320 GeV has been performed [97]. The final state qqqq£~^i~ 
resulting from HZ — > ZZZ and from HZ W~^W~ Z is selected. A kinematic fit is used to assign 
the possible di-jet combinations to bosons (VF"'"VF~ or ZZ). The resulting di-boson mass spectrum 
can be fitted by a Breit-Wigner distribution convolved with a detector resolution function. A 
relative uncertainty on the Higgs mass of 0.11-0.36% is achievable from 500 fb~^at 500 GeV for 
Higgs masses between 200 and 320 GeV. The resolution on the total width varies between 22 and 
34% for the same mass range. Finally, the total Higgs-strahlung cross-section can be measured with 
3.5-6.3% precision. Under the assumption that only H WW and H ZZ decays are relevant, 
their branching ratios can be extracted with 3.5-8.6% and 9.9-17.3%, respectively (see Table I3TT|) . 
The expected mass spectra for mn = 200 GeV and niH = 320 GeV are shown in Fig. 13.201 

3.6.3 Open questions 

In the field of heavy SM-like Higgs bosons many questions remain for future studies. Besides 
extending the results quoted above to other (possibly) higher Higgs mass values the following 
issues should be addressed: 

• What are the chances to observe couplings of the heavy Higgs bosons other than those already 
investigated? 



32 




Di-boson mass (GeV) Di-boson mass (GeV) 



Figure 3.20: Expected reconstructed Higgs boson mass spectra for nifj = 200 GeV and tjih = 
320 GeV from 500 fb~^at 500 GeV (from [97]). 



mH (GeV) 


Aa (%) 


AruH (%) 


Ar^^ (%) 


200 


3.6 


0.11 


34 


240 


3.8 


0.17 


27 


280 


4.4 


0.24 


23 


320 


6.3 


0.36 


26 



Table 3.11: Expected precision on Higgs boson line-shape parameters for 200 < mn < 320 GeV at 
the ILC with = 500 GeV. 



• What can be learned about the quantum numbers and the spin of heavy (SM-like) Higgs 
bosons? 

• Do these additional measurements pose new requirements for the detector? 

3.7 Heavy Neutral Higgs Bosons in the MSSM 

Besides the light CP-even Higgs boson in the MSSM there are two other neutral Higgs bosons: 
the heavier CP-even H and the CP-odd A. For Ma ^ 150 GeV the two heavy Higgs bosons are 
very similar in mass. The couplings of the Higgses to gauge bosons is either strongly suppressed 
(for HVV) or even zero at the tree-level (for AW). Therefore the analyses performed for a heavy 
SM-like Higgs boson as outlined in the previous section cannot be taken over to the case of heavy 
MSSM Higgs bosons. 

Similar to the light CP-even Higgs boson, the heavy neutral CP-even Higgs boson of MSSM 
can be produced in e~^e~ collisions via Higgs-strahlung and weak boson fusion processes. Neutral 
Higgs bosons of MSSM can be also produced in pairs, e~^e~ — > hA and e~^e~ HA, and singly 
through Yukawa processes, e.g., e+e" h{A, H)bb; see Fig. 12.61 In the decoupling limit of the 
MSSM, the Higgs-strahlung and weak boson fusion processes involving H, as well as pair production 
process involving h, get strongly suppressed (see, however, Ref. [49]). Hence, in this scenario the 
main production mechanisms of H and A bosons at an e~^e~ collider will be Yukawa processes and 
e'^e~ HA. The pair production of heavy neutral Higgs bosons at an e~^e~ collider has been 
studied in Ref. [98] . The analysis is based on selection of bbbb and bbT~^T~ final states, which are 
expected to be dominant channels for e~^e~ — > HA process for a large part of MSSM parameter 
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space. For the reconstruction of these final states, efficient 6-tagging is important. The analysis 
employs the same flavor-tagging tools [99] as used in the study of the hadronic decays of the 
light Higgs boson [70]. The potential of the linear e~^e~ collider for detecting heavy neutral Higgs 
bosons and measuring their properties has been evaluated for various Higgs boson mass hypotheses. 
Although in large parts of the MSSM parameter space H and A are almost degenerate in mass, 
scenarios with large mass splitting between neutral Higgs bosons may occur. This feature is present 
in CP-conserving as well as in CP-violating scenarios. 

Studies have been performed for a linear collider operated at center-of-mass energies of 500 and 
800 GeV. An integrated luminosity of 500 fb~^ is assumed for each energy. The branching fractions 



e 

sh 
O 

s 

3 

e 




120 140 160 180 200 220 

reconstructed ditau mass [ GeV ] 



e 

> 

sh 
O 

U 

01 

s 

3 

e 



0, 



^ Simulated data 

□ HA^bbT*T" 
H HA^xVbb 

in„=200GeV m,=lSOGeV 




120 140 160 180 200 220 

reconstructed dijet mass [ GeV ] 



Figure 3.21: Distributions of the reconstructed di-tau mass (left plot) and di-jet mass (right plot) 
in the HA — > bbr^ t~ , t~ bb channels for {mH^friA) = (200,150) GeV at ^/s = 500 GeV. Assumed 
integrated luminosity is 500 fb~^. 

of the Higgs bosons into 6-quarks and r-lepton pairs are set to their typical values in the MSSM: 
BR(i/ ^ bb) = BR{A bb) = 90%, BR(iJ ^ rr) = BR(A ^ rr) = 10%. Furthermore, the 
maximal allowed cross section for e'^e~ — > HA is assumed, corresponding to | sin(a — /3)| = 1. These 
assumptions, however, do not restrict the generality of the study, and results of the analysis can 
be applied for arbitrary values of Higgs branching fractions by appropriate rescaling the expected 
number of signal events in each analyzed channel. As an example Fig. 13.211 presents reconstructed 
ditau and dijet mass spectra in the bbT^T~ channel for the Higgs boson mass hypothesis {mH^rnA) 
= (200,150) GeV. For the same mass hypothesis the distribution of the reconstructed Higgs boson 
mass sum and mass difference in the bbbb channel is presented in Fig. 13.221 Exploiting these two 
channels, the precision on Higgs boson masses and production cross sections is evaluated as a 
function of the Higgs masses. Higgs boson masses can be measured with accuracy ranging from 0.1 
to 1 GeV for Higgs pair production far above and close to the kinematic limit, respectively. The 
topological cross section a^+^-^^ji^^^i^i can be determined with relative precision of 1.5—7% and 
cross sections (T„4 



^HA^bbr+T- 



*HA^bbbb 

and (T„ 



*HA- 



-bb 



with a precision of 4—30%. Furthermore, 
Higgs boson widths can be determined with relative accuracy ranging from 20 to 40%, depending 
on the Higgs masses. In the case of mass degenerate H and A bosons, the 5a discovery reach 
corresponds to a common Higgs boson mass of about 385 GeV in the HA — > bbbb channel at t/s = 
800 GeV, assuming an integrated luminosity of 500 fb^^. The mass reach for heavy neutral Higgs 
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Figure 3.22: Distributions of the reconstructed Higgs boson mass sum, niA + mn, (left plot) and 
mass difference, \mH — tuaI, (right plot) in the HA — > bbbb channel for {mH,mA) = (200,150) GeV 
at ^/s = 500 GeV. Assumed integrated luminosity is 500 fb~^. 



bosons can be extended by operating the linear collider at higher centre-of-mass energies. Fig. 13.231 
shows as an example the signal in the HA bbbb channel at ^/s = 1 TeV for the Higgs boson 
mass hypothesis (m/f,mA)=(394.6,394.9) GeV, which corresponds to the SPS la benchmark point. 
With 1 ab^^ of data collected at ^/s = 1 TeV, uih and uia can be measured with an accuracy of 
1.3 GeV, and the topological cross section with precision of 9% [98]. 
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Figure 3.23: Distribution of the 

reconstructed Higgs boson mass sum, 
niH + ruA, for the SPS la benchmark 
point. Higgs boson mass hypothesis 
is (mH,mA)=(394.6,394.9) GeV. Distribu- 
tion corresponds to an integrated luminos- 
ity of 1 ab~^ collected at -y/s = 1 TeV. 



Alternatively, heavy neutral Higgs bosons can be produced resonantly in 77 collisions. The 
main advantage of the photon collider with respect to the e~^e~ machine is its higher mass reach 
for heavy neutral Higgs particles. Neutral Higgs bosons with masses up to 350 GeV can be detected 
with high statistical significance at a photon collider operated at a modest center-of-mass energy 
of -s/s^ of 420 GeV [100]. More generally, Higgs bosons can be produced with masses up to 80% 
of the total e~^e~ energy. An example of the signal in the ^ H,A ^ bb channel is presented 
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Figure 3.24: Distribution of the invariant 
mass of the hadronic system in the 77 
H,A bb channel. Distribution corre- 
sponds to 808 fb^^ of luminosity collected 
with photon collider operated at center- 
of-mass energy of ^/s^ = 419 GeV [100]. 
Higgs boson masses are mn = ruA = 300 
GeV. 



in Fig. 13.2^1 With this process a 7C can, e.g., cover the "LHC wedge" where only a SM-like Higgs 
boson is within the reach of other colliders. 

As was discussed above, in the decoupling limit of MSSM one of the most promising channels 
to detect heavy neutral Higgs bosons will be associated Higgs boson production, e^e~ — > HA. 
Unfortunately, the mass reach for neutral Higgs bosons in this channel is limited: the Higgs bosons 
can be produced only if mn+mA < \fs. H and A bosons can be also produced in Yukawa processes, 
such as e'^e~ — > H{A)bb. In this case the mass reach is extended almost up to collision energy. 
Although a number of theoretical studies have been performed which investigated the prospects of 
detecting heavy neutral Higgs states at e~^e~ collider in Yukawa processes, no related experimental 
analysis including realistic detector simulation and SM background estimates has been done so far. 
We hope that this issue will be addressed by future studies. 

3.8 Charged Higgs Bosons 

The charged Higgs boson (H^) is clear evidence of extended Higgs sectors beyond the SM. In the 
MSSM, the mass of the charged Higgs boson is related to that of the CP-odd Higgs boson by 
^loop' so t^^t and A are almost degenerate in the large mass limit. By 
measuring this mass relation, the MSSM may be distinguished from the general two Higgs doublet 
model in which no such stringent relationship holds. While the derivative couplings VH^H~, 
W^H^A and W^H^H appear at tree level, there is no tree level H^W^V = 7 or Z) coupling 
due to the C/(l)cm gauge invariance and also the global SU{2)v symmetry in two Higgs doublet 
models including the MSSM. 

At the ILC, a relatively light charged Higgs boson is produced in pairs, e^e^ H^H~, 
as long as the mass is smaller than one half of the collider energy, mjj± < \/s/2. A complete 
simulation of this process with the decay tb and H~ bt has been performed in Ref. [101] 

for = 800 GeV. The mass reconstruction is done by selecting the final state to be the four b 
jets plus four non-6 tagged jets from the two W bosons from top decay. The expected signal and 
background are shown in Fig. 13.251 for = 300 GeV, assuming the integrated luminosity to 
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Figure 3.25: Fitted charged Higgs boson 
mass for H~^H~ (tb){tb) with Tn^± = 
300 GeV. The histogram is normahzed to 
an integrated luminosity of 1 ab~^, with a 
100% branching ratio in the analyzed de- 
cay mode. The contribution of the tbtb 
background events is shown by the dark 
histogram. 

"200 220 240 260 280 300 320 340 360 380 400 

m (H) in GeV 

rec^ ^ 

be 1 ab-i. The mass resolution is approximately 1.5%. A 5a level discovery will be possible for 
m^± < 350 GeV for y/s = 800 GeV. 

While the mass reach of pair production is limited to \/^/2, the rare processes of single charged 
Higgs production may be useful for the search of heavier charged Higgs bosons. The dominant 
processes are e^e~ — > btH~^, e^e~ — > tv't-H'^, e^e~ W~H~^ and e^e~ H^e~i>e (one-loop 
induced). Their production rates have been calculated at leading order in Refs. [31-34]. QCD 
corrections to e"'"e~ — > btH~^ have been studied in Ref. [32], and have turned out to be sizable. For 
more details about the progress in the cross section calculations, see Sect. 14.21 In Ref. [102], the 
simulation for the process e~^e~ — > tbH~ + c.c. — > 4:b+ii+i+p^^^^ {I = e, fi) has been performed for 
a charged Higgs boson mass above half of the collider energy, where the pair production channel 
e'^e~ H~H~^ is no longer available. It has been shown that one can establish a statistically 
significant signal only over a rather limited mass region, of 20 GeV or so beyond AIjj± ~ \fs/2, 
for very large or very small values of tan/3 and provided high 6-tagging efficiency can be achieved; 
see Fig. HTIEl 

Figure 3.26: (Solid and dashed lines) Dif- 
ferential distribution in the reconstructed 
Higgs mass from both 5-jets not generated 
in top decays and the two top systems for 
the signal e+e" biH+ + tbR- Mb 
and the background e'^e~ tig* ttbb, 
yielding the signature, after the kinematic 
cuts, including all decay BRs. (Dotted and 
dot-dashed lines) Same as above, but using 
only the 6-jet with highest transverse mo- 
mentum of the two not generated in top 
decays. See Ref. [102] for details. 

The situation is different for charged Higgs masses smaller than mt. A precise measurement 
of Mjj± in that range is a challenging task for any collider. In this case the charged Higgs decays 
dominantly to a tau lepton and neutrino {H^ — > t^v), and it is impossible to reconstruct directly 
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the invariant mass of the di-tau final state. However, due to the polarization of the leptons, 
the energy of the decay products depends strongly on Mf^±, a feature that can be exploited 
to extract Mjf± indirectly. The main background to a signal in the decay mode comes 
from the W-boson decays — > t^u. However, thanks to the different structure of the and 
electroweak interactions to r^-leptons, the tau polarization is very different which will allow 
a separation of and ^ T^u on a, statistical basis. More specifically, the decay 

products {t^ vr+z?, or r"*" — > p^v, etc.) have strikingly different topologies according as to 
whether they originate from a parent or H^. 

The key point is that leptons arising from t^u decays are almost purely right-handed, 

in contrast to the left-handed leptons which arise from decays. This contrast follows from 
the helicity flip nature of the Yukawa couplings of the Higgs fields, and the helicity conserving 
nature of the gauge interactions. The most dramatic difference is seen in the energy distribution 
for the single pion channel {H^ /W~^ t~^v tt'^uv) in the rest frame of the parent boson {W^ 
or H^) [103]. In practice, it is nearly impossible to reconstruct the rest frame of the or 
bosons because the momentum of the neutrinos cannot be measured. Instead, one can use the top 
quark pair production channel with top quark decays {t — ^ H'^h) and reconstruct the top quark 
rest frame. 

The distributions for the energy of the single pion in the top rest frame has the following form 
for and cases, respectively [104-106], 



I dV 1 \ {l-Pr)logf^ + 2Pry^{-^-^), < < 

r dy^ ~ Xmax - Xmrn 1 (1 " PT)log^ + 2P^(1 - ^), X^in < 



(3.15) 



where 1I — r ■ — =^ r — =^ j^min _ _™s_ jpmax _ mop y pr^ boson 

wneie (/tt — j^^^^ , Xmin — Mtop •^'^"■^ ~ Mtop ' T ~ 2Mtop' " 2 " "^"^ ^^^^ OObOn, 

Ft = —1, and for the charged Higgs boson, Pr = 1. 

The energy distribution for a pion from decay has a maximum at the point -B(7r^) = 
Mfj±^ /{2Mt). This dependence allows one to extract Mjj± from the shape of the spectrum. 

The QCD corrected top-quark decay width (in terms of qH+ = m'jj+/mf) is given by [107], 

r^Soit - hH+) = -^^^M^ - <lH^?m,\mj) tan^ p x 



TT 



' - ^ - 21og(l - g^+) + 2(1 - 
+ (^ + \ log(l - qH+)] (1 - QH+f 



(3.16) 



The Amf, corrections have been included using CPsuperH [108]. The numerical simulations have 
been performed assuming a center of mass energy = 500 GeV and a total integrated luminosity 
C = 500 fb"^. We have performed detailed computations and Monte Carlo simulations for different 
sets of MSSM parameters, all leading to light charged Higgs bosons [103]. Here we present two 
scenarios, based on the mass parameters Mq = Mu = Md = 1 TeV, M~ = M2 = 1 TeV, At = 500, 
II = ±500 GeV, and tan/3 = 50, which give Mh± = 130 GeV. The branching ratio BR{t H^b) 
will be enhanced or suppressed depending on the sign of /j,. The branching ratios for these examples 
are 

(i) Mh± =130 GeV, /x < and tan/3 = 50: BR{t ^ H+b) = 0.24 

(ii) Mh± = 130 GeV, /x > and tan (5 = 50: BR(t H+b) = 0.091 

The couplings have been implemented in CompHEP [109], which has been used to compute the 
cross sections for signal and background processes, including decays of top to and which 
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subsequently decay to polarized r leptons. CompHEP was also used to generate events, and ef- 
fects from initial state radiation and Beamstrahlung were included. Polarized decays have been 
simulated using TAUOLA [110] interfaced to CompHEP. Hadronization and energy smearing in the 
final state are accounted for by means of PYTHIA [111] using the CompHEP-PYTHIA [112] inter- 
face based on Les Houches Accord [113]. Effects from final state radiation have been implemented 
using the PHOTOS library [114]. 

A brief description of the fitting procedure is given in Ref [103]. We use the method of maximum 
likelihood to fit a spectrum created from simulated signal and background events, and obtained 
the following results for the cases described above: 

(i) Mh± = 129.7 ± 0.5 GeV, 

(ii) Mh± = 129.4 ± 0.9 GeV, 

Figure l5.27l shows results of the simulations and fits. The charged Higgs mass may be determined, in 
both cases, with an uncertainty of the order of 1 GeV. This study is a theoretical level analysis, and 
no systematics or detector eff^ects have been included. A full detector simulation will be necessary 
to test the robustness of this result. 
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Figure 3.27: Generated vr^ energy spectra (left plots) and the fit (right plots) for the two sets of 
MSSM parameters described in the text. 



3.9 Measuring tan f3 in 77 Collisions 

The measurement of the mixing parameter tan /?, one of the fundamental parameters in the Higgs 
sector of supersymmetric theories, is a difficult task, especially if tan /? is large. In this section two 
methods are described to measure tan/? at the 7C. 
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rr fusion to Higgs bosons at a photon collider [115] can provide a valuable method for measuring 
this parameter. The couplings in the limit of large tan (3 (and moderate values of Ma in the case 
of the h) are given by (see e.g. Ref. [116]) 

g^T-T = tan fj for = A 

g^rr ^ tan/3 iov<^ = h,H (3.17) 

[normalized to the Standard Model Higgs coupling to a tau pair, mr/v]. Thus the entire Higgs 
mass range up to the kinematical limit can be covered for large tan /3 by this method. 
We consider the process 

'y-f ^ TT TT ^ TT + h/H/A . (3.18) 

For large tan/3 all the Higgs bosons decay almost exclusively [80 to 90%] to a pair of b quarks, and 
the final state consists of a pair of r's and a pair of resonant b quark jets. The main background 
channels are r+r" annihilation into a pair of 6-quarks, via s-channel "f/Z exchanges, naturally 
suppressed as a higher-order electroweak process, and diffractive 77 {T~^T~){bb) events, reduced 
kinematically by choosing proper cuts. The SUSY parameters have been chosen according to the 
SPS lb scenario [117], but Ma has been varied. 




Figure 3.28: The cross sections for the production of the H/A Higgs bosons in the rr fusion process 
at a 77 collider for tan /? = 30. The other parameters have been chosen according to the SPS lb 
scenario (but with Ma varied). Also shown is the background cross section for properly chosen 
experimental cuts, ^/s denotes the 77 collider cm. energy, corresponding to approximately 80% of 
the e^e~ linear collider energy. 

The cross sections for the signals of H and A Higgs-boson production in the rr fusion process, 
together with the background processes, are presented in Fig. I3.28[ including proper cuts. The 
cross section for rr fusion to the light Higgs boson h is of the same size. 

The statistical accuracy with which large tan/? values can be measured in rr fusion to h/H/A 
Higgs bosons is exemplified for three tan/3 values, tan/3 = 10, 30 and 50, in Table Results 
for scalar H production are identical to pseudoscalar A in the mass range above 120 GeV. The 
two channels h and A, and H and A are combined in the overlapping mass ranges in which the 
respective two states cannot be discriminated. In Table IXT^ the relative errors Atan/3/ tan/3 are 
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Table 3.12: Relative errors Atan/3/tan/? on tan/3 in measurements for tan /? = 10, 30 and 50, 
based on: combined A®h [I,IV] and A®H [II,III,V-VIII] production, assuming [E^^ = 400 GeV, 
C = 100 and [E^^ = 600 GeV, C = 200 fb~^]. Proper cuts and efficiencies are applied on the 
final-state r's and h jets. 

presented. Since in the region of interest the rr fusion cross sections are proportional to tan^ j3 
and the background is small, the absolute errors A tan (5 are nearly independent of tan (3, varying 
between 

Atan/3~0.9 and 1.3 (3.19) 

for Higgs mass values away from the kinematical limits. 

This analysis compares favorably well with the corresponding 6-quark fusion process at the 
LHC [118]. Moreover, the method can be applied readily for a large range of Higgs mass values 
and thus is competitive with complementary methods in the e^e~ mode of a linear collider [119]. 

An alternative method to measure tan/3 is provided by the associated tH^ production in 77 
collisions [120]. The subprocess 67 — > H~t utilizes the 6-quark content of the photon, tan/3 enters 
through the thH^ vertex, so that the amplitude squared for the subprocess 67 — > tH~ is given by: 

^|M(67 oc V2GF^a{mltiin^ + (3). (3.20) 

Fig. 13.291 shows the cross section as a function of tan /3 with the measurement precision superim- 
posed. It has been found that 77 tH^ + X can be used to make a good determination of tan (3 
for most of the parameter space with the exception of the region around tan /? ~ 7 where the cross 
section is at a point of inflection. This measurement provides an additional constraint on tan /3 
which complements other processes. 

3.10 Measurements in Exotic Scenarios 

An interesting possibility is that the Higgs decays into particles not present in the Standard Model. 
In particular, if there are exotic light objects which did not obtain their mass directly from the 
Higgs mechanism, their couplings to the Higgs may be large despite their small masses, and they 
may dominate over the standard Higgs decay modes. Such a scenario depends very strongly on how 
these new objects manifest themselves inside the detector, and can play a very significant role in how 
well both the LHC and the ILC can determine the Higgs properties. In fact, the tiny expectation 
for the SM Higgs width (provided Mh ^ 150 GeV) implies that the Higgs is particularly sensitive 
to new decay modes, even if they are not very strongly coupled. 
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3.10.1 Higgs decaying into unflavored jets 

The first such scenario arises when the Higgs decays primarily into light hadrons, without any 
particular bottom or charm content. This actually occurs in some limits of the minimal super- 
symmetric standard model [121]. If, for example, the pseudo-scalar Higgs A is light enough, the 
dominant decay mode of the SM-like Higgs may he H ^ AA. If the A mass is less than 2mb, and 
tan (3 is larger than 10 or so, the dominant A decay will be into strange quarks, with separation 
too small to be resolved as separate jets. Thus, the decay H — > AA — > ssss results in two hard jets 
of energy roughly Mh/'2' in the Higgs rest frame^. 

A second possibility is that the Higgs decays into down-type squarks [123]. The most likely 
candidate is the scalar bottom 6, which could further decay into a pair of anti-quarks through a 
baryon-number violating interaction of the type UDD. Again, if the mass of the b is relatively 
small, the two quarks may not be sufficiently separated to be resolved as individual jets, and the 
Higgs decays effectively into two jets. Such a light bottom squark was postulated in [124] to address 
an excess in the open bottom quark production cross section at run I of the Tevatron. The MSSM 
Higgs width into a pair of scalar bottoms can be compared with the SM width into bb as, 



where fj, is the higgsino mass parameter, tan/? is the ratio of Higgs VEVs, and 6), is the sbottom 
mixing angle. The decay into b dominates whenever /itan/3/m/i is large, and the mixing angle is 
not too small and the b mass not too large. 

However the new decay mode arises, this scenario can be modelled by assuming that the Higgs 
receives a new contribution to its width into jets (and assuming the width into other SM particles 
remains SM-like). For definiteness, we discuss the case with the light b, but the results are largely 
insensitive to this choice. As the width into the new decay mode becomes large, it eventually 
overwhelms the SM decay modes, driving their branching ratios to be small. Figure 13.301 shows the 
branching ratios for all of the important decay modes as a function of /x tan (3/mh- One should note 
that the light b have been included in the loop-generated H-g-g and H-^-^ couplings, and thus 
there is some model-dependence in those curves. As can be seen from the figure, once /itan/3/m/i 
is larger than about 10, the decay into b completely dominates, with the canonical SM decay into 
bb having a branching ratio of less than 10%. 

The primary feature of the new decay mode is that, by reducing the decay into SM particles 
while leaving the Higgs production cross sections unchanged, it reduces the statistics available to 

"^This situation is also likely to arise in the NMSSM for regions of parameter space with small fine-tuning [122]. 




(3.21) 
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Figure 3.30: Branching fractions for various 
Higgs boson decay channels as a function of 
the ratio /xtan/3/m/i, with mh = 120 GeV. 
The sbottom mass is fixed to mr = 5 GeV. 



measure the Higgs interactions with SM particles and thus limits the precision with which they 
may be measured. At the same time, the decay into jets is experimentally very challenging, and can 
prevent the LHC from discovering the Higgs by depressing the decay into observable SM particles 
below LHC sensitivities and opening a new decay channel which cannot be seen at the LHC because 
of huge QCD jet backgrounds. A careful analysis shows that if the width into the new hadronic 
state is more than a few times the width into bottom quarks, the LHC completely loses the ability 
to detect the light CP-even Higgs [123]. This could lead one to the perplexing situation in which 
supersymmetry and perhaps one or more of the MSSM heavy Higgses are discovered, but the Higgs 
actually responsible for electroweak symmetry breaking is too difficult to extract from the large 
hadronic backgrounds. 

The ILC will easily detect such a Higgs, because its primary search mode is not sensitive to the 
Higgs decay mode (see Sect. IXT|) . At the same time, the high sensitivity expected at the ILC will 
help allow it to make important measurements of the Higgs. As the decay width into jets increases, 
the statistics for SM decays decreases, and some measurements are degraded. However, since the 
decay into jets can in general be extracted from the background, measurements of the couplings 
associated with Higgs production remain accessible and precise. Figure 13.311 plots the projected 
uncertainties in the Higgs couplings as a function of the ratio of widths T{H — > jj)/T{H bb). 
Note that H — > jj includes the SM decay into bottom quarks as well as the new contribution into 
unflavored jets. 

3.10.2 Invisibly decaying Higgs 

The Higgs might also decay predominantly into a neutral exotic particle, such as the lightest 
neutralino of a supersymmetric theory. In such a case, the Higgs becomes effectively invisible, 
because the decay products do not interact with the detector material. However, even in this 
case a precise mass measurement is possible with the recoil method down to the level of about 
~ 100 MeV [8]. Concerning the other (now rare) decay channels, the existing studies apply with 
appropriately scaled signal rates; see Sects. 13.31 13.61 
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Figure 3.31: Expected accuracy in the mea- 
surements of the bb and jet-jet branch- 
ing fractions, the hZZ and hWW coupling 
strengths, and the total width of the Higgs 
boson, as a function of the ratio of the jet-jet 
and the bb widths, assuming the Higgs boson 
couplings to bb, ZZ, and WW* are standard. 
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Chapter 4 

Theoretical Developments 



The ILC does not only offer an even greater physics potential, but in turn represents a great 
challenge for theorists to understand phenomena at the experimentally achievable level of precision. 

At energies exceeding the reach of LEP2, many new processes will be accessible, such as top- 
quark pair production, Higgs production, or reactions with new-physics particles, as predicted, 
e.g., by SUSY models. To exploit the potential of the ILC, predictions for such reactions should be 
based on full transition matrix elements and improved by radiative corrections as much as possible. 
The higher level of accuracy at the ILC requires a corresponding level of precision in the theoretical 
predictions. 

In the following we focus directly on Higgs physics, but there are other important areas (related 
to EWSB) where higher-order corrections are indispensable. This comprises precision calculations 
for M\Y (i.e., for decay); precision observables on the Z resonance (e.g., for the effective leptonic 
weak mixing angle, sin^^eflf); radiative corrections to 2 — 3, 4, . . . processes in general; single W 
production; etc. 

4.1 Radiative corrections to Higgs masses and couplings in the 
MSSM 

Within the MSSM the masses of the CT^-even neutral Higgs bosons as well as the masses of the 
charged Higgs bosons are calculable in terms of the other MSSM parameters. The status of the 
available corrections to the masses and mixing angles in the MSSM Higgs sector (with real pa- 
rameters) can be summarized as follows. For the one-loop part, the complete result within the 
MSSM is known [15,125-127]. The by far dominant one-loop contribution is the 0{at) term due 
to top and stop loops {at = /if/(47r), ht being the superpotential top coupling). Concerning the 
two-loop effects, their computation is quite advanced and has now reached a stage such that all 
the presumably dominant contributions are known. They include the strong corrections, usually 
indicated as 0{atas), and Yukawa corrections, 0{af), to the dominant one-loop 0{at) term, as 
well as the strong corrections to the bottom/sbottom one- loop 0{ai,) term (ctt = /i^/(47r)), i.e., 
the 0{a{,as) contribution. The latter can be relevant for large values of tan /3. Presently, the 
0{atas) [16,128-131], O(a^) [128,132,133] and the 0(a(,as) [134,135] contributions to the self- 
energies are known for vanishing external momenta. In the (s)bottom corrections the all-order 
resummation of the tan /3-enhanced terms, 0(ab(as tan /?)"), is also performed [107,136]. Recently 
the 0{ataf)) and 0{a'^) corrections became available [137]. Finally a full two-loop effective po- 
tential calculation (including even the momentum dependence for the leading pieces) has been 
published [138]. However, no computer code is publicly available. 
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In the case of the MSSM with complex parameters (cMSSM) the higher order corrections have 
so far been restricted, after the first more general investigations [139], to evaluations in the effective 
potential (EP) approach [140, 141] and to the renormalization group (RG) improved one-loop EP 
method [142,143]. These results have been restricted to the corrections coming from the (s)fermion 
sector and some leading logarithmic corrections from the gaugino sector. A more complete one-loop 
calculation has been attempted in Ref. [144]. More recently the leading one-loop corrections have 
also been evaluated in the Feynman-diagrammatic (FD) method, using the on-shell renormalization 
scheme [145]. The full one-loop result can be found in Ref. [146,147]. 

The upper limit is estimated to be rrih ^ 135 GeV [16,17]. The remaining theoretical uncertainty 
on the lightest CT^-even Higgs boson mass has been estimated to be below ~ 3 GeV [17,52]. The 
above calculations have been implemented into public codes. The program FeynHiggs [16, 17, 37, 
148] is based on the results obtained in the FD approach, it includes all available corrections. The 
code CPsuperH [108] is based on the RG improved effective potential approach. For the MSSM 
with real parameters the two codes can differ by up to ~ 4 GeV for the light CP-even Higgs boson 
mass due to subleading higher-order corrections that are included only in FeynHiggs. 

4.2 Single Higgs production at one-loop in 2 — ^3 processes 

Recently some one-loop calculations of electroweak radiative corrections have been presented for 
2 — > 3 processes that are interesting at the ILC: e^e~ — > vvH [149,150] and e'^e~ tiH [151-153]. 
The results of Refs. [149,152] were obtained with the Grace-Loop [154] system (see below). In 
Refs. [150,151,153] the technique [155] for treating tensor 5-point integrals was employed. While 
Refs. [149,151,152] make use of the slicing approach for treating soft-photon emission, the results 
of Refs. [150, 153] have been obtained by dipole subtraction and checked by phase-space slicing for 
soft and collinear bremsstrahlung. 

In e~^e~ annihilation there are two main production mechanisms for the SM Higgs boson. In 
the Higgs-strahlung process, e~^e~ ZH, a virtual Z boson decays into a Z boson and a Higgs 
boson. The corresponding cross section rises sharply at threshold {^/s ^ Mz + M}{) to a maximum 
at a few tens of GeV above Mz + Mh and then falls off as s~^, where \fs is the CM energy of 
the e~^e~ system. In the VK-boson-fusion process, e~^e~ v^VeH, the incoming e"*" and e~ each 
emit a virtual W boson which fuse into a Higgs boson. The cross section of the VF-boson-fusion 
process grows as Ins and thus is the dominant production mechanism for ^/s ^ Mh- A complete 
calculation of the 0{a) electroweak corrections to e"'"e~ — > vvH in the SM has been performed 
in Refs. [149,150]. The results of Refs. [149,150] are in good agreement, see, e.g., [156]. The 
agreement of the correction is within 0.2% or better with respect to the lowest-order cross sections. 

Sample results taken from Ref. [150] are shown in Fig. 14.11 for two Higgs boson masses, Mh = 
115, 150 GeV. The higher-order corrections are of the order of a few percent if the tree-level result 
is expressed in terms of [150]. 

The Yukawa coupling of the top quark could be measured at a future ILC with high energy and 
luminosity at the level of ~ 5% [8] by analyzing the process e^e~ tiH. A thorough prediction 
for this process, thus, has to control QCD [157] (possibly SUSY QCD [158]) and electroweak 
corrections [152,153]. The results of Refs. [152,153] have been compared with other calculations, 
see e.g. Ref. [156]. Agreement within ~ 0.1% is found. This also holds for other energies and 
Higgs-boson masses. The results of the previous calculation [151] roughly agree with the ones of 
Refs. [152,153] at intermediate values of ^/s and Mh, but are at variance at high energies (TeV 
range) and close to threshold (large Mh)- 

Sample results taken from Ref. [152] are shown in Fig. l4.2l for Mh = 120 GeV. The higher-order 
corrections are of the order of 10 percent if the tree-level result is expressed in terms of [150]. 
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Figure 4.1: Relative electroweak corrections to the complete process e~^e~ — > vi/H and to the 
contributions from Zi?-production and TVPF- fusion channels for Mh = 115 GeV and Mh = 
150 GeV [150]. 
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Figure 4.2: (a) Total cross section as a function of the CMS energy for Mh = 120 GeV. Shown 
are the results for the total cross sections for the tree level, full 0{a) and 0{as) level in (a). The 
relative corrections are shown in (b). Solid lines are tree level, dashed lines are the full 0{a) and 
dotted lines are the 0{as) corrections. In addition, the genuine weak correction 6w and the relative 
correction 6-^^ in the scheme are presented. For more details see Ref. [152]. 
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More higher-order corrections possibly have to be evaluated in order to reach the required ILC 
accuracy. 

We now turn to the corresponding predictions within the MSSM At the one-loop level, first the 
contributions of fermion and sfermion loops in the MSSM to e~^e~ vuh have been evaluated in 
Ref. [49,159]. The corrections to the light Higgs boson production have been found to be at the 
percent level [49]. 

Since at the tree level and in the decoupling limit the heavy neutral MSSM Higgs bosons 
decouple from the Z, the mass reach for their discovery at the ILC is limited to approximately 
\/s/2 from the pair production process. It has been investigated how single production mechanisms 
could extend the mass reach of the ILC. In particular, the VFiy-fusion process e^e~ — > u^VeH has 
been investigated [49]. Its tree level cross-section is proportional to cos(/3 — a). Depending on 
the SUSY parameters, radiative corrections might increase the cross-section for e'^e~ — > Vei'eH, 
possibly allowing discovery beyond the pair production kinematic limit for certain choices of the 
MSSM parameters. Using left-polarized electron beams and right-polarized positron beams the 
cross-section can further be enhanced. A particular scenario where this is the case has been chosen 
in [49] (MsusY = 350 GeV, /i = 1000 GeV, M2 = 200 GeV and large stop mixing). Cross-section 
contours for this scenario are shown in Fig. 14.31 




400 500 500 700 800 400 500 500 700 80 



Ma Ma 

Figure 4.3: Cross-section contours for e^e^ Huv for a particular MSSM scenario (see text) in 
the ruA — tan (3 plane for ^/s = 1 TeV. The different shadings correspond to: white: a < 0.01 fb~^, 
light shaded: 0.01 fb-^< a < 0.02 fb~\ dark shaded: 0.02 fb-^< a < 0.05 fh'^ black: a > 0.05 
fb^Hfrom [49]). The left fi gure is for unpolarized beams, the right figure for an electron (positron) 
polarization of 0.8 (0.6). 

In a CP violating scenario the three neutral Higgs bosons, Hi, H2, H3, are mixtures of the CP 
even and CP odd Higgs fields. Consequently, they all couple to the Z boson and to each other. 
These couplings may be very different from those of the CP conserving case. In the CP violating 
scenario the Higgs-strahlung processes e~^e~ — > HiZ {i = 1,2,3) and pair production processes 
e"'"e~ — > HiHj [i 7^ j) may all occur, with widely varying cross-sections. 
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4.3 Higgs production at one-loop in 2 



2 processes 



The most promising channels for the production of the CP-even neutral MSSM Higgs bosons in 
the first phase of an ILC are the Higgs-strahlung processes [160], 

e-^e-^ZHi, (4.1) 

{Hi^2 = h, H) and the associated production of a scalar and a pseudoscalar Higgs boson, 

e+e-^AHi. (4.2) 

The full one-loop corrections [161] and the leading two-loop corrections [162] are available. They 
have been combined to obtain the currently most accurate results for these production cross sections. 
The results have been implemented into the code FeynHiggsXS [162]. 

In Fig. 14.41 the dependence of azh and aAh on are shown in the "no-mixing" benchmark 
scenario [163] (MsusY = 2 TeV, Xt = 0, At = Ab, /i = 200 GeV, Ms = 200 GeV, m-g = O.SMsusy) 
for tan (3 = 5 and tan (3 = 50. The solid (dot-dashed) line represents the two-loop result including 
(excluding) box contributions, the dotted line shows an effective coupling approximation and the 
dashed line shows the one-loop result. For azh, the A boson decouples quickly; the dependence 
on M4 becomes very weak for Ma ^ 250 GeV, when azh is already practically constant. In the 
same limit, a Ah goes quickly to zero due to suppression of the effective ZhA coupling, which is 
~ cos(aeff — also the kinematical suppression plays a role, but this becomes significant only 
for sufficiently large Ma, Ma > 350 GeV. For large tan/3 the decoupling of Ma is even more 
rapid. The differences between the FD two-loop result and the effective coupling approximation 
for the Higgs-strahlung cross section tend also to a constant, but they increase with Ma for the 
associated production. The latter can be explained by the growing relative importance of 3- and 
4-point vertex function contributions compared to the strongly suppressed Born-like diagrams. As 
can be seen from Fig. 14. 4( for tan /3 = 50 and Ma > 300 GeV the two-loop result is almost an order 
of magnitude larger than the result of the effective coupling approximation, and starts to saturate. 
This can be attributed to the fact that the (non-decoupling) vertex and box contributions begin to 
dominate the cross section value. However, such a situation occurs only for very small aAh values, 
aAh ~ 10~^ fb, below the expected experimental ILC sensitivities. 

Comparing the one-loop and two-loop results an uncertainty larger than ~ 5% can be attributed 
to the theory evaluation of the production cross sections. Another type of possibly large corrections 
in the MSSM are the so-called Sudakov logs (see Ref. [164] and references therein). They appear in 
the form of \og{q'^ / M^^jgy) (where q is the momentum transfer) in the production cross sections of 
SUSY particles at e~^e~ colliders. In order to reach the required level of accuracy of 1-2% further 
two-loop corrections as well as the corresponding Sudakov logs will have to be calculated. 

4.4 Double Higgs production processes at one-loop order 

Full electroweak radiative corrections to the double Higgs-strahlung e~^e~ — > ZHH has been calcu- 
lated in Refs. [80,165]. The two results are in good agreement, within 0.1%, for ^/s ^ 800 GeV for 
niH = 115, 150, 200 GeV. In Ref. [80], the computation is performed with the help of GRACE-loop. 

It has been shown that the QED corrections are large especially around the threshold, so that 
a proper resummation of the initial state radiation needs to be performed in this case. At energies 
where the cross section can be substantial, however, the corrections are modest, especially for a 
Higgs boson mass of Mh — 120 GeV. The genuine weak corrections near the peak of the cross 
sections are also not large, not exceeding ~ 5%, and therefore are below the expected experimental 
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Figure 4.4: azh and a Ah as a function of Ma at ^/s = 500 GeV, shown for tan /3 = 5 and tan (3 = 50, 
in the no-mixing scenario (see text), as evaluated with FeynHiggsXS. The sohd (dot-dashed) hne 
represents the two-loop result including (excluding) box contributions, the dotted line shows an 
effective coupling approximation and the dashed line shows the one- loop result. 



precision. As discussed in Ref. [77,78], the invariant mass Mhh distribution can be useful to extract 
the triple Higgs boson coupling. It is found that the genuine weak corrections, contrary to the QED 
corrections, hardly affect the shape of the distribution at least for energies where this distribution 
is to be exploited [80]: see Figure EH and also Figure EiHl (left). Therefore, an anomalous triple 
Higgs boson coupling could still be distinguished, if large enough, in this distribution provided that 
a proper inclusion of the initial QED corrections is allowed in the experimental simulation. 
Recently, the corrections to the double Higgs production were also calculated in Ref. [166]. 
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Figure 4.5: da/dMuH for Mh = 120 
GeV at -y/s = 1 TeV. We show the tree- 
level (full curve), the effect of includ- 
ing only the genuine weak corrections 
(dotted curve) and the effect of includ- 
ing the full 0{a) corrections (dashed 
curve) [80]. 
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4.5 Heavy MSSM Higgs bosons with CP violation at the 77 col- 
lider 

The 2-doublet Higgs sector of the MSSM is automatically CP conserving at the tree level. However, 
CP violation can be induced by radiative corrections transmitting CP violating phases from the 
soft SUSY breaking Lagrangian to the Higgs system, in particular the relative phases between the 
higgsino mass parameter and the trilinear sfermion-Higgs parameters Af. The Higgs system 
in 2-doublet models is generally described by a 3 x 3 complex mass matrix, composed of a real 
dissipative part and an imaginary absorptive part. The real dissipative part includes the parameters 
of the Higgs potential extended by loop corrections. The imaginary part includes the corresponding 
absorptive contributions of the loops, i.e., the widths in the diagonal elements. 

CP asymmetries are naturally enhanced [167, 168] in the decoupling regime where Ma — Mjj- 
Near mass degeneracy of the scalar and pseudoscalar states in CP violating theories allows for 
frequent mutual transitions which induce large CP-odd mixing effects. In this limit the CP violation 
effects are located in the heavy Higgs sector and the light Higgs field can be ignored^. The mass 
matrix of the heavy Higgs sector can be reduced to a 2 x 2 form: 

The CP violation effects arc encoded in A|^^ connecting the scalar H and the pseudoscalar A 
fields through dissipative and absorptive contributions. Diagonalizing this mass matrix leads to 
the CP-mixed mass eigenstates H2 and i?3. The rotation is described by the complex mixing 
parameter 

1 

X = - tan 20 = _ ^2 _ ^^MhVh - MaTa] ^^'^^ 

The shift in the masses and widths and the CP coupling of the states are quantitatively accounted 
for by the mixing parameter X: 

1 rr 1 r rxVl + iX^ 



(4.5) 



^In general, however, it would be desirable to include the full 3x3 mixing (with complex entries). 
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For a typical SUSY scenario, 



SUSY 



0.5 TeV ^ Ma, \At\ = 1.0 TeV, 



1.0 TeV, 



; tan (3 = 5 



(4.6) 



the mixing parameter X and the shifts of masses and widths are illustrated in Figs. I4.6r a). (b) and 
(c). (pA denotes the CP-violating phase of At. 
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Figure 4.6: The dependence on the phase of At, (pA, of (a) the mixing parameter X and of the shifts 
of (b) masses and (c) widths for the parameter set (|4.6I) : in (b,c) the mass and width differences 
without mixing are shown by the broken lines. 3^e/9mX(27r — c^a) = +3^e/— 9mX((/)^) for angles 
above tt (taken from Ref. [167]). 

Several asymmetries can be defined which reflect the presence of CP violating mixing effects 
either indirectly or directly. Direct CP violation can be studied in a classical way by measuring 
asymmetries of inclusive cross sections between left- and right-polarized photons [167]. At a photon 
collider, the circular polarization of laser photons is transferred completely to the high energy 
photons in Compton back-scattering of laser light near the maximum of the photon energy [8, 169, 
170]. 

The CP violating asymmetry is defined by 



A 



hel 



(4.7) 



where the subscripts denote the helicities of the two colliding photons. On top of the H2, Higgs 
states the asymmetries are parameterized by the complex mixing parameter. 



2Im(cos^sin0*) 
I cos^l^ + I sin^l^ 



(4.8) 



in the asymptotic decoupling limit. Even though corrections due to non- asymptotic Higgs masses 
are still quite significant for the SUSY scenario introduced above, remarkably large asymmetries 
are predicted nevertheless in 77 fusion, shown in Fig. 14.71 as a function of the CP violating phase 
of the trilinear stop-Higgs coupling At- 

Thus a photon collider provides an instrument that allows us to perform a unique classical 
experiment on CP violation in the Higgs sector of supersymmetric theories. 



4.6 Charged Higgs Bosons 

Charged Higgs bosons can be pair-produced at the ILC via e+e~ H^H'^ if m/^± < ^/s/2 [171]. 
A complete simulation of this process for the decay — > th has been performed for -y/s = 800 
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Figure 4.7: The 4)a dependence of the 
CP-violating asymmetry ^hei on top of 
the H-i Higgs states, respectively. 



GeV, 1 ab~^, and mij± = 300 GeV [101]. The mass resolution is approximately 1.5%. A 5o" 
discovery will be possible for mfj± < 350 GeV. 

Since in pair production the mass reach for charged Higgs bosons is limited to \/s/2, also the 
rare processes of single charged Higgs production may be considered. The dominant processes 
for single charged Higgs production are e~^e~ btH~^,e'^e~' t"i7t-H^, and e^e~ — > W^H^. 
Their cross-sections have been calculated at leading order in Refs. [31,32]. Although the H^W^ 
production is a 2 — > 2 process [33], the cross section is suppressed because of the absence of the 
tree level H^W^Z coupling in multi Higgs doublet models. QCD corrections to e~^e~ — > btH~^ 
have recently become available [32] and are sizable. In general, parameter regions for which the 
production cross-section exceeds 0.1 fb are rather small for charged Higgs masses beyond the pair 
production threshold. Cross-section contours for ^/s = 500 GeV and 800 GeV are shown in Fig. 14.81 
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Figure 4.8: Cross-section contours for the processes e^e btH^ (blue/dark), e^e r Ut-H^ 
(green/light grey), and e+e" — > W~H~^ (red/medium grey) at -y/i = 500 GeV (left) and at 800 
GeV (right) [93]. 

At a 77 collider, production of a heavy charged Higgs boson is typically one order of magnitude 
larger than the production of the same final state at an e"*"e~ collider with the same collision 
energy due to the t-channel nature. For ^/s^^ > 27X1 jj± , pair production is dominant similarly to 
the electron-positron collision [171, 172]. When the mass of the charged Higgs boson is larger than 
the threshold of pair production, single production processes become important. Cross sections 
of various single production processes have been calculated in Ref. [173]. Moreover, a polarized 
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77 collider can determine the chirality of the Yukawa couplings of fermions with the charged 
Higgs boson via single charged Higgs boson production and, thus, discriminate models of new 
physics [174]. 

4.7 NMSSM 

4.7.1 Overview 

One attractive super symmetric model is the Next to Minimal Supersymmetric Standard Model 
(NMSSM) [175, 176] which extends the MSSM by the introduction of just one singlet superfield, 
S, in addition to the Higgs up- and down-type superfields if„ and if^- [Hatted (unhatted) cap- 
ital letters denote superfields (scalar superfield components).] When the scalar component of S 
acquires a TeV scale vacuum expectation value (a very natural result in the context of the model) , 
the superpotential term XSHuH^ generates an effective fiHuHd interaction for the Higgs doublet 
superfields with ji = //eflf = M^) = Such a term is essential for acceptable phenomenology. No 
other SUSY model generates this crucial component of the superpotential in as natural a fashion. 
At the same time, the NMSSM preserves all the most attractive features of the MSSM. First, 
since the NMSSM introduces only an additional singlet superfield, the SU{3), SU{2)l and ^7(1) 
coupling constants continue to unify at the unification scale, Mu ^ few x 10^^ GeV, just as in 
the MSSM. Further, radiative electroweak symmetry breaking remains a natural possibility within 
the NMSSM. The electroweak symmetry breaking generates non-zero vevs for all the Higgs fields, 
{Hu) = hu, {H(ii = h(i and {S) = s. A possible cosmological domain wall problem [177] can be 
avoided by introducing suitable non-renormalizable operators [178] that do not generate danger- 
ously large singlet tadpole diagrams [179]. Alternatively, it has been argued that the domain walls 
are themselves the source of dark energy [180]. 
The trilinear term in the Higgs potential reads 

XAxSHuHd + '^A^S^ . (4.9) 

Aside from A, k, /leS, and A\, the final Higgs sector parameter is tan/3 = hu/hd- Thus, as 
compared to the three independent parameters needed in the MSSM context (often chosen as ^, 
tan/3 and m^), the Higgs sector of the NMSSM is described by the six parameters 

X , K , Ax , A^, tan/3 , fies ■ (4-10) 

A possible sign convention [181] is such that A and tan/3 are positive, while k, A\, A^ and fj,eS 
should be allowed to have either sign. 

The particle content of the NMSSM differs from the MSSM by the addition of one CP-even 
and one CP-odd state in the neutral Higgs sector (assuming CP conservation), and one additional 
neutralino. Altogether, in the Higgs sector we have the CP-even Higgs bosons hi,h2,h3, the CP- 
odd Higgs bosons ai,a2 and the usual charged Higgs pair H"^. The five neutralinos (Xi,2,3,4,5) 
typically mixtures of the gauginos, higgsinos and the singlino. In some limits, the singlino can be 
a fairly pure state and it can either be the heaviest or the lightest of the neutralinos. 

4.7.2 The Higgs boson mass spectrum 

The mass of the lightest CP-even Higgs boson at the tree-level is given by 

mluee < Ml [ cos^ 2/3 + sin2 2/3 ] , (4.11) 
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where A should be below A < 0.7 to remain perturbative up to the GUT scale. The evaluation 
of higher-order corrections to Higgs boson masses and mixing angles is less advanced than in the 
MSSM. No full one-loop calculation is available so far although there has been much work [182- 
187]. Most recently the leading two-loop Yukawa corrections have become available [188]. All 
available corrections have been implemented into the code NMHDECAY [181]. No genuine higher- 
order corrections to Higgs boson production and decay are available. However, in many cases 
(especially if the new singlet mass scale is large) the results from the MSSM can be taken over to 
the NMSSM case. 

An important issue for NMSSM Higgs phenomenology is the mass and nature of the lightest 
CP-even and CP-odd Higgs bosons. In particular, if the ai is very light or even just moderately 
light there are dramatic modifications in the phenomenology of Higgs discovery at both the LHC 
and ILC [121,181,187,189-194]. A light oi is natural in the context of the model, see Ref. [121] 
for a discussion. A very light ai is experimentally not ruled out, especially if it has a large singlet 
component. In some regions of the NMSSM parameter space, one can also get the lightest CP-even 
state, hi, to be very light as well. 

Another interesting feature of the NMSSM is that the standard measure of fine-tuning 



F = MaXpFp = Maxp 



dlogMz 



dlogp 



(4.12) 



where the parameters p comprise all GUT-scale soft-SUSY-breaking parameters, can be very low 
for Higgs boson mass values above the LEP limit [122, 176, 195, 196]. 



4.7.3 LEP Higgs searches 

Within the NMSSM it is also possible to have a Higgs boson with a mass below the LEP limits [122]. 
It escapes LEP constraints by virtue of having unusual decay modes for which LEP limits are 
weaker. In particular, parameters for which hi aiui decays are dominant over hi bb decays 
must be carefully considered^. Particular care is required if there is a mixture of hi aiUi and 
hi bb decays, and the dominant ai decay is to bb. In this situation, the single channel Z2b and 
Z4:b rate limits from LEP are less sensitive than the combination of both channels, since in the 
Z2b interpretation the Z46 final state is ignored, and vice versa, which leads to a loss in sensitivity 
compared with the full combination. 

This has particular relevance to the NMSSM scenarios with the very lowest fine-tuning, F < 10, 
as examined in Ref. [199]. Two types of scenarios emerge for F < 10. In both, BR(/ii — > bb) ~ 
0.1 0.2 and BR(/ii aioi) ~ 0.85 0.75. In scenarios of type I (II), > 2m;, {nia^ < 2mh) 
and BR(ai bb) ~ 0.92 (0). More details can be found in Refs. [198, 199]. It has been found that 
the only F < 10 scenarios consistent with the full LEP Higgs Working Group (LHWG) analysis 
are of type II. Type I scenarios are typically excluded at about the 98 99% CL after data from 
all experiments are combined. 

A very important feature of the type-II -F < 10 scenarios is that a significant fraction of 
them can easily explain the well-known excess in the Zh — > Zbb final state in the vicinity of 
ruh ~ 100 GeV. This is illustrated in Fig. El [199], where all F < 10 NMSSM parameter choices 
(from a lengthy scan) with niai < ^m^ are shown to predict nih^ ~ 96 105 GeV along with 
a Ceff = [dzzh/ 9zzhsM^^^^^ ~^ ) would explain the observed excess with respect to the 
expected limit. The statistical significance of this excess is of order 2a. The properties of the four 
most ideal points for describing the excess are given in Table H?T1 Note that for point 3, rua^ < 2mr 

^This is also possible in the MSSM [197]; corresponding holes from the LEP Higgs searches in the MSSM parameter 
planes exist [198]. 
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implying that oi quarks and gluons. A general feature with F < 10 is that they predict the 
Higgs boson properties needed to describe the LEP Zh Z + 6's excess. Another possibility, 
occuring at higher F values, is that the hi has a ZZhi coupling that is reduced in strength relative 
to the SM ZZ-Higgs coupling. This can give rise to other scenarios with light Higgs bosons that 
are consistent with the LHWG limits on the Z + 6's final state. 
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Figure 4.9: Expected and 

observed 95% CL limits on 

Cos = [glzJahhsM^^m - bb) 

from Ref. [4] are shown vs. ruh. 
Also plotted are the predictions for 
NMSSM parameter choices in our 
scan that give fine-tuning measure 
F < 25 and niai < 2mij with fixed 
tan /? = 10 and gaugino masses of 
Mi,2,3(Mz) = 100,200,300 GeV. [199] 
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Armax 


(GeV) 


(GeV) 


hi bb 


hi oioi 


ai — > r+r 


300 fb"^ 
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3.6 
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0.938 


2.4 
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3.1 


0.141 


0.832 


0.0 


2.5 


102.2 


3.6 


0.177 


0.789 


0.814 
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Table 4.1: Some properties of the hi and oi for four points with F < 10 and < 2mf, from 
a tan/3 = 10, Mi,2,3(Mz) = 100,200,300 GeV NMSSM scan. iV^^^(300 fb"^) is the statistical 
significance of the best "standard" LHC Higgs detection channel after accumulating integrated 
luminosity of L = 300 fb~^. 



4.7.4 ILC measurements 

In this section we analyze the implications for the ILC of an NMSSM scenario in which there is a 
Higgs with fairly SM-like couplings, but that decays to two lighter Higgs bosons, emphasizing the 
most probable case of /i — > ao and nih near 100 GeV. 

To set the stage for the ILC, an important question is the extent to which the type oi h ^ aa 
Higgs scenario (whether NMSSM or other) described here can be explored at the Tevatron and the 
LHC. This has been examined in the case of the NMSSM in [189-192, 194], with the conclusion 
that observation of any of the NMSSM Higgs bosons may be very difficult at hadron colliders if 
the decay mode hi aioi is dominant. For example, the largest statistical significance in any 
of the "standard" search channels at the LHC, given in the last column of Table W?^ for the four 
sample type-II scenario points, is below ba even after L = 300 fb~^ of integrated luminosity is 
accumulated. In this case, if any signals are observed at the Tevatron or LHC, it seems very 
possible that they will be quite tenuous, perhaps difficult to interpret, and possibly of quite low 
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statistical significance. In such a circumstance, tlie ILC becomes essential. It is the ideal machine 
for detecting a mh ~ 100 GeV Higgs boson with unusual decays. One simply employs the search 
channel of inclusive e~^e~ — > Z* — > Zh — > i~^i~X in which Mx is computed from the incoming 
and outgoing lepton information. An unmistakable peak in Mx is predicted, and it is completely 
independent of the h decay mode. Detection of 77 ^ /i —> aa is also predicted to be completely 
straightforward (regardless of the dominant a decay mode) at the 77 collider facility that will 
probably accompany the ILC [200]. Typical signal and background plots are shown in Ref. [200] 
assuming that the 77 collider is operated at an electron beam energy of 75 GeV and with laser and 
electron polarizations chosen to produce the general purpose broad spectrum in E^y. 

It should be noted that much of the discussion above regarding Higgs discovery is quite generic. 
Whether the a is truly the NMSSM CP-odd ai or just a lighter Higgs boson into which the SM-like 
h pair-decays, hadron collider detection of the h in its h ^ aa decay mode will be very challenging 
— only the ILC can currently guarantee its discovery. 

4.7.5 Open questions 

As mentioned above, the NMSSM studies are in general less advanced than the corresponding 
MSSM studies. The following open questions should be answered so that the ILC results can be 
fully exploited in the context of the NMSSM: 

• Calculations for the Higgs masses, mixing angles and couplings have to be refined to reach the 
required level of accuracy. This task is similar to the MSSM case, where, however, already 
more precise calculations exist. Also the renormalization in the NMSSM Higgs sector is less 
advanced than in the MSSM. More precise calculations are also needed for the NMSSM Higgs 
boson production cross sections and branching ratios. 

• In the important low fine-tuning case that hi aioi with ai t^t^ is dominant, excellent 
r-tagging capabilities at the ILC are necessary. Further experimental studies are needed. In 
a more general case where hi aiai and oi bb, excellent 6-tagging capabilities would be 
crucial. This includes the distinction between b and b in order to reduce the number of b-jet 
combinations in multi-6 final states. 

• As in the MSSM, also in the NMSSM the lightest neutralino is a good candidate for cold dark 
matter, see e.g. [201,202] and references therein. It is possible that Higgs physics plays an 
important role here if the LSPs annihilate via a Higgs boson to SM particles. More precise 
calculations for the Higgs boson couplings are necessary to reach the required precision. More 
experimental efforts have to be undertaken to extract the corresponding Higgs boson masses 
and widths with the required accuracies. This is discussed further in Sect, d There, the 
difficult scenario of annihilation of the lightest neutralino via a rather light singlet-like CP- 
odd Higgs boson (as possible in the low-F NMSSM scenarios) is discussed. 

• There are other, non-Higgs-sector-related issues. These comprise the question of baryogenesis 
within the NMSSM, the investigation of the LSP quantities at the ILC, and the determina- 
tion of neutralino sector observables in order to reconstruct all corresponding fundamental 
parameters. 

4.8 Gauge Extensions of the MSSM 

Another way to increase the mass of the lightest SUSY Higgs boson can arise from theories with 
extended gauge interactions. If the Higgs is charged under one or more of the new gauge groups. 
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there are additional D-teim contributions to its quartic term in the Higgs potential, and thus the 
light CP even Higgs may be considerably heavier than the tree-level MSSM expectation. Early 
efforts in this direction included additional U{1) gauge groups [203]. These are effective, but the 
possibility of large Higgs masses is limited if one requires perturbativity up to the grand unification 
scale. More recently, non-Abelian gauge groups which are asymptotically free have appeared which 
allow for large Higgs masses, but whose interactions are free from Landau poles to arbitrarily high 
scales [204,205]. 

A simple example is provided by supersymmetric topflavor, which has gauge structure SU (3)c x 
SU{2)i X SU{2)2 X C/(l)y. The MSSM Higgses and the third family are charged under 5C/(2)i 
whereas the first two generations are charged under SU {2)2- Precision EW constraints require the 
masses of the extra gauge bosons to be greater than a few TeV. The Higgs mass at tree level may 
be expressed as, 

mlo < ^[g^A + g^)v^cos^(3 (4.13) 

where g and gy are the usual MSSM SU{2)l and U{1)y gauge couplings, and A parameterizes the 
new contribution from the extended SU{2) structure, 

^ - J ? (4-14) 
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where m? is a soft mass for the Higgs which breaks SU {2)i x SU {2)2 SU{2)l, u is the breaking 
scale, and gi and 52 are the underlying SU{2) couplings, satisfying l/gl + 1/(72 — ^Idh- Tree level 
Higgs masses on the order of 350 GeV are relatively easy to obtain in this model, which clearly 
evades the LEP bounds [4]. This motivates studies of heavier (SM-like) Higgs bosons, see Sect. 13.61 

The model contains the extra Z' and W's from the broken extended gauge symmetry. These 
couple more strongly to the third family than to the first two and the W is expected to be discovered 
at the LHC if this model is realized in nature [206]. The ILC could further potentially pin down 
the model and its parameters by indirectly measuring the properties of the Z' , and studying the 
neutralino and chargino sectors, which now include the partners of these new gauge bosons as well 
as the usual MSSM components. Detailed studies of the branching ratios, production cross sections, 
etc., would be interesting to pursue. 

The presence of extra asymptotically free gauge interactions may also result in an interesting 
modification of the NMSSM. The new strong gauge interactions may drive the coupling A small 
at high energies, effectively removing the perturbative unification bound on the value of A (see 
Sect. HT^ . and realizing an otherwise theoretically unmotivated region of parameter space [207]. 
A sample Higgs spectrum (at one loop) is presented in Figure B.lOl and shows several interesting 
features including a rather heavy set of CP even fields, and a charged scalar as the lightest of all 
of the Higgses. 



4.9 Fat Higgs Models 

A final possibility is to interpret the strongly coupled NMSSM as composite theory which confines 
at the strong coupling scale [208-210]. The original Fat Higgs [208] introduces an additional SU{2) 
interaction and a number oi preons charged under it. At an intermediate scale. A, the SU{2) theory 
confines, forming a singlet superfield S and Higgs-like superfields H and H which play the role of 
the two MSSM Higgses. There is a dynamical super-potential induced of the form, 

Wdyn = \s(hH-vI) (4.15) 
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Figure 4.10: Higgs spectrum in an NMSSM-like model with large A [207]. 



where naive dimensional analysis predicts A ~ 47r at the confinement scale, and Vq ~ mA/Air 
is related to the mass of the pair of underlying preon fields (m) which bind into S. The mass 
parameter m is ultimately what determines the electroweak scale. 

Notably different from the NMSSM, the contribution to the Higgs quartic A is much larger 
because the effective cut-off is much lower than the unification scale^. In addition, the electroweak 
symmetry is broken even in the supersymmetric limit. Higgs spectra similar to those expected 
from the gauge extensions with large A result, with similar phenomenology quite distinct from the 
MSSM and usual NMSSM; see Figure ETTUl for an example. 

Variations of the original Fat Higgs either restrict the composite sector to include only the 
gauge singlet [209] or generate the top quark from the strong dynamics as well as the Higgses [210], 
motivating the fact that the top mass is much larger than any of the other fermions by postulting 
that the large top Yukawa interaction is a residual of the force that bound a set of preons into both 
top and Higgs. 

4.10 The Significance of tan/3 in Two-Higgs-Doublet Models 

There are a number of discussions in this report about the measurement of the parameter tan/? 
at the ILC. A precise theoretical understanding of the significance of this parameter is important. 
Recent progress in this direction is reviewed. 

4.10.1 tan /3 as an (un)physical parameter 

In the two-Higgs-doublet model (2HDM), the Higgs sector of the Standard Model is extended to 
include two complex hypercharge-one Higgs doublets [25,211]. After the breaking of electroweak 
symmetry (preserving the electromagnetic U(1)em symmetry group), the Higgs fields' vacuum 
expectation values (vevs) can be chosen to take the following form: 




(4.16) 



^Unification of the MSSM gauge coupHngs remains possible (and calculable, despite the strong dynamics) provided 
some addtional matter is included to arrange for the unification. 
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where < ^ < 27r and the Vi are real and positive and normahzed such that v'^ = vf -\- V2 = 
Arny^/g^ = (246 GeV)^. In addition, the ratio of the magnitudes of the two vevs is denoted by 
tan/3 = V2/vi^ where < /? < 7r/2. 

The relative phase E, is clearly not a physical parameter of the model, since one is free to rephase 
the $2 field to absorb this phase. Perhaps more surprisingly, tan [3 is also not generally a physical 
parameter. Since the two Higgs fields and ^2 possess identical gauge quantum numbers, one 
is free to redefine the scalar fields via an arbitrary U(2) transformation [212,213], — > Uij^j. 
Different choices of the unitary matrix U correspond to different basis choices for the scalar fields. 
Physical predictions of the model must, of course, be basis independent. 

If tan P is not a physical parameter, what is the meaning of tan /? that appears so often in this 
report? In fact, tan/3 can be promoted to a physical parameter by imposing some restrictions on 
possible terms that appear in the 2IIDM Lagrangian. For example, by imposing supersymmetry 
on the model, the two Higgs fields and <I>2 are distinguishable in a particular basis, and this 
basis provides a physical definition of tan /?. However, if phenomena consistent with a two-Higgs- 
doublet sector are observed, it may not be immediately obvious what theoretical constraints govern 
the underlying 2HDM. Thus, in order to carry out a truly model-independent study of 2HDM 
phenomena, one cannot employ tan /? as defined below eq. (|4.16|) . 

Here we consider the most general 2HDM. One natural basis choice that exists is the so-called 
Higgs basis [213-218], consisting of the scalar fields and in which only one of the neutral 
Higgs fields (which by convention will be denoted by has a non-zero vacuum expectation value. 
This basis choice is not unique, since one is always free to independently rephase the individual 
Higgs fields. One linear combination of these two rephasing transformations (corresponding to 
$a e^'^^a, e^^^b) is the global hypercharge transformation U(1)y- This transformation 

does not modify any of the scalar potential parameters. However, we shall employ it simply to 
remove the phase from the vacuum expectation value of In this convention, the Higgs basis is 
parameterized by one angle x- 

$a = $1 cos /? + e-*^ ^>2 sin /? , (4.17) 
e'^^f, = -e*«^>isin/? + $2Cos/?. (4.18) 

One then obtains 

= ( ^ (. +T? + ^G^) ) ' = ( (X+ ^A) ) • (4.19) 

Indeed only develops a vacuum expectation value, which by convention is chosen to be real and 
positive. Of the original eight scalar degrees of freedom, three Goldstone bosons {G^ and G^) are 
absorbed ("eaten") by the and Z. In eq. ()4.19|) . ipi, are CP-even neutral Higgs fields, A is 
a CP-odd neutral Higgs field, and is the physical charged Higgs boson pair. If CP is conserved, 
then A is the physical CP-odd Higgs scalar and ip^ and (^2 to produce two physical neutral 
CP-even Higgs states h and H. If the Higgs sector is CP- violating, then ipi, ip2, and A all mix to 
produce three physical neutral Higgs states of mixed CP properties. 
In the Higgs basis, the scalar potential takes the following form: 

+ ^Al(cI>tcI>,)2 + iA2(cI>tcI>,)2 + A3(<I>i$a)(^Il>6) + A4(^I^6)(^3;i>a) 

+ I^Ag e2^x(cDt<|>^)2 + [Ag e'^i<fl<i>a) + A7 e'^i<fl<^>b)M^b + h-cj , (4.20) 
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M^-^ and fixed by the scalar potential minimum conditions, and is directly related to 

the physical charged Higgs boson squared-mass: 

ml± = + ^v^As . (4.21) 

The conditions for CP-invariance of the Higgs sector are very simple. The 2HDM is CP- 
conserving if and only if the Higgs basis scalar potential parameters satisfy [213, 217] 

Im (AgAg) = Im (A^A?) = Im (A^Ar) = . (4.22) 

If eq. ()4.22|) is satisfied, then there exists a choice of x such that all the Higgs basis scalar potential 
parameters are simultaneously real [219-221]. By working in the Higgs basis, it is clear that 
tan/3 never appears. Thus, the couplings of the Higgs bosons to the vector bosons and the Higgs 
self-couplings can depend only on the scalar potential parameters (in the Higgs basis), the gauge 
couplings and the mixing angles that determine the neutral Higgs boson mass eigenstates in terms 
of ipi, and A. 

In practice, the parameter tan /3 is most often associated with the Higgs- fermion Yukawa cou- 
plings. Thus, we now focus on the coupling of the 2HDM to three generations of quarks and leptons. 
The Higgs couplings to fermions are model dependent. The most general structure for the 2HDM 
Higgs-quark Yukawa couplings, often referred to as the type-Ill model [222], is given in the Higgs 
basis by (with an analogous form for the Higgs- lepton Yukawa couplings): 

-Cy = + Ql'^a^''''D% + Q° + qI'^,p''''D% + h.c. , (4.23) 

where <1> = ia2^*, is the weak isospin quark doublet, and ?7j^, are weak isospin quark 
singlets, k^'^ and p^''^ {Q = U , D) are Yukawa coupling matrices in quark flavor space. 

The fermion mass eigenstates are related to the interaction eigenstates by the following unitary 
transformations : 

Ul = VEuI, UR = V}iul, Dl = VEdI, Dr = v£d%, (4.24) 

and the Cabibbo-Kobayashi-Maskawa matrix is defined as = V^^V^K It is also convenient to 
define "rotated" linear combinations of the Yukawa coupling matrices: 

(4.25) 

The quark mass terms are identified by replacing the scalar fields with their vacuum expectation 
values. The unitary matrices V^, V^, Vj^ and are chosen so that and are diagonal with 
real non-negative entries. These quantities are proportional to the diagonal quark mass matrices: 

Md = -^K^ , Mu = -^K^ . (4.26) 

V 2 v 2 

However, the matrices p^ and p^ are complex non-diagonal matrices, which are not constrained in 
a general model. 

The Higgs-fermion interaction, expressed in terms of quark mass eigenstates, takes the following 
form [211,213,223]: 

-Cy = -DMoDifl + 1=D{p''Pr + p''^Pl)D^I + ^Dip^'Pn - p""^ Pl)DA + -DMdI.^DG'^ 
V y 2 y 2 v 

+-UMuU^l + -^Uip'^PR + p^^Pl)U^I - -^U{p^Pr - P^^Pl)UA - -UMdI^UG^ 
■y V2 V2 V 



U 



Kp^Pr - p^^KPl] DH+ + [KMdPr - MuKPl] DG+ + h.c. | , (4.27) 
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where Pr^l = |(1 ^Ts)- The physical couphngs of the quarks to the Higgs fields are then obtained 
by expressing ^p\, </?2 and A in terms of the neutral Higgs mass eigenstates. Once again, the 
parameter tan (3 does not appear. 

The most general 2HDM is phenomenologically untenable over most of its parameter space. 
In particular, since and are non-diagonal, eq. H4.27() [even when expressed in terms of the 
neutral Higgs mass eigenstates] exhibits tree-level Higgs-mediated flavor changing neutral currents 
(FCNCs). One way to avoid this problem is to go the region of parameter space where M22 is 
very large (or equivalently mjy± <C mz). In this decoupling limits the second doublet $2 is very 
heavy and decouples from the effective low-energy theory. To a very good approximation, Lp^ is a 
CP-even Higgs mass eigenstate with couplings nearly indistinguishable from those of the Standard 
Model Higgs boson [224,225]. In particular, tree-level Higgs-mediated FCNCs are suppressed, and 
the corresponding 2HDM is phenomenologically viable (at the expense of the fine-tuning necessary 
to guarantee that the second Higgs doublet is very heavy). 

The more common procedure for avoiding FCNCs is to impose a symmetry on the theory in order 
to restrict the form of the Higgs- fermion interactions in a particular basis [226]. This procedure 
then establishes the symmetry basis. If the scalar fields of the symmetry basis are denoted by 
$1 and ^2, whose relation to the Higgs basis is given by eqs. (|4.17)) - (|4.18j) . then the parameter 
tan /3 is promoted to a physical parameter of the theory. For example, a discrete symmetry can be 
imposed such that ^2 decouples from the fermions in the symmetry basis. This yields the so-called 
type-I Higgs-fermion couplings [227,228]. Another discrete symmetry can be found in which down- 
type right-handed fermions couple exclusively to $1 and up- type right-handed fermions couple 
exclusively to <I>2. This yields the so-called type-II Higgs-fermion couplings [215,228]. As a final 
example, the MSSM also exhibits type-II Higgs-fermion couplings [229]. 

For simplicity, consider a one- generation type-II CP-conserving model. It is straightforward to 
show that the Higgs basis Yukawa couplings satisfy the following equation [213]: 

K^K^ + p^p^ = 0. (4.28) 

Moreover, one can compute tan /3 (where /3 is the relative orientation of the Higgs basis and the 
symmetry basis): 

tan/3=^ = -^, (4.29) 
where = \/2Mq/v. Thus, it follows that: 

p = tan p , p = cot p . (4.30) 

V V 

Inserting this result into eq. ()4.27() . and expressing and P2 terms of the CP-even Higgs mass 
eigenstates h and H yields the well-known Feynman rules for the type-II Higgs-quark interactions. 



4.10.2 Interpretation of ILC measurements 

How should one interpret precision measurements of 2HDM Yukawa couplings? A priori, one does 
not know whether a symmetry basis exists. Thus, one should attempt to determine experimentally 
the p^ and p^ coupling matrices. Testing relations such as eq. (|4.28j) could reveal the presence 
of a symmetry basis, in which case one can determine tan/3 from eq. ()4.29l) . However, one could 
discover that the equality of eq. (|4.29|) does not hold. For example, in a one-generation type-Ill CP- 
conserving model, a more general procedure would consist of defining three tan /3-like parameters: 

_pD 

tan I3d = — ^ , tan/3„ = — , (4.31) 
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and a third parameter tan/3e = —p^/n^ corresponding to the Higgs-lepton interaction. The 
equahty of these three tan /3-hke parameters would imply the existence of a symmetry basis (and 
an underlying type-II Higgs- Yukawa coupling). Deviations from these equalities would imply a 
more complex realization of the electroweak symmetry breaking dynamics. 

We can illustrate this procedure in a very simple model approximation. In the MSSM at large 
tan/3 (and supersymmetric masses significantly larger than mz), the effective Lagrangian that 
describes the coupling of the Higgs bosons to the third generation quarks is given by: 

- £eff = h{qL^i)hR + ht{qL^2)tR + Ahb {qL^2)bR + h.c. , (4.32) 

where = {ul , di)- The term proportional to Ahb is generated at one-loop due to supersymmetry- 
breaking effects. The tree-level relation between and hb is modified [107,136]: 

rub = ^ cos P{l + Ab), (4.33) 

where Af, = (A/tft/Zift) tan/3. That is, A^ is tan /3-enhanced, and governs the leading one-loop 
correction to the physical Higgs couplings to third generation quarks. In typical models at large 
tan/3, Af, can be of order 0.1 or larger and of either sign. An explicit expression for A;, is given in 
Refs. [107,136,230]. In the approximation scheme above, ~ /it sin /3, ~ /if cos /3, and [213] 

~ hb cos /3(1 + Ab) , ^-hb sin /3 fl ^) ~ -hb sin /3 . (4.34) 

V tan^ /3 / 

It follows that: 

tan 3 

tan (5d — T— , tan /3u ~ tan /3 . (4.35) 

1 + Aft 

Thus, supersymmetry-breaking loop-effects can yield observable differences between tan /3-like pa- 
rameters that probe the underlying nature of supersymmetry-breaking. 



4.11 Lepton Flavor Violation in the Higgs boson decay at ILC 

Lepton flavor violation (LFV) in charged leptons directly indicates new physics beyond the SM. 
In the models such as based on SUSY, in addition to the gauge boson mediated LFV process, the 
LFV Yukawa couplings are naturally induced by slepton mixing [231]. The direct search for such 
LFV Yukawa couplings would be possible by measuring the Higgs decay process h, {H, A) — r^/i^ 
at collider experiments [232-234]. Such a possibility has been examined at the LHC in Ref. [235]. 
In Ref. [236] , this possibility has been discussed and the feasibility of this process at a ILC via the 
Higgs-strahlung process has been studied under the constraints from current data on LFV rare tau 
decay processes. A possibility of a search for the LFV Yukawa interaction via the deep inelastic 
scattering process e~N ^ t~X at the fixed target option of ILC has also been studied [237]. 



4.12 Dimension-Six Higgs Operators 

4.12.1 The SM as a low energy effective theory of new physics models 

Under the constraint from the current precision data, the SM with a light Higgs boson would be 
an effective theory to the cutoff scale which could be considered to be as high as the Planck scale 

without contradicting theoretical arguments such as vacuum stability and triviality. On the other 
hand, such scenario of the SM with very high cutoff scale has been considered to be unnatural. 
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because due to quantum corrections the Higgs boson mass is quadratically sensitive to the cutoff 
(new physics) scale (~ A^). Demanding relatively small (less than 90%) hierarchy, the scale of the 
new physics naturally would be expected to be of the order of 1-10 TeV. The effect of new physics 
would then enter into the low energy theory as higher dimensional operators. If a light Higgs boson 
is discovered at the LHC but no additional particles are seen at the LHC or the ILC, it is important 
to search for small deviations of the Higgs boson potential from the SM predictions to probe new 
physics scales. If the reason for such small deviations is beyond-the-SM physics at large scales 
A, the effective operator approach can be chosen to parameterize the low-energy behavior of such 
models. 

The effective Lagrangian of the SM with the cutoff A can be described as 

>Cefr = >CsM + E J^^n + 0{l/A^+% (4.36) 

n 

where On are dimension six operators, and /„'s are dimcnsionlcss "anomalous couplings". It is 
expected that /„'s are of order of unity at the cutoff scale A. There are eleven SM gauge invariant 
bosonic operators in dimension six (dim-6) terms [238] , and eight of them are relevant to the Higgs 
boson field: 



Obw 




(4.37) 




= (D^$)^$^$(£)^$), 


(4.38) 


Oww 




(4.39) 
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(4.41) 


Ob 




(4.42) 




= ^5^($t$)a^($t$)^ 


(4.43) 






(4.44) 



These operators are separated to the following three categories. 



[Class I] (Obm',0*,3) 

These operators contribute to the gauge boson two point functions {O^^^, S, 3 ^ T), so that 
they receive tight constraints from the current precision data: their coefficients are bounded at the 
95% C.L. as [239] 



[Class II] (0„,„.,0££.0„,, Ob) 

These four operators lead to interactions of the Higgs boson with gauge bosons such as 

+9'klzHZ^Z^^^ + 9^h\vw{W+W-^^9''H + h.c.) + g^^^^HW+W'^^r (4.46) 

Combinations of the coefficients of C^KW' C^^, and contribute to these V-H couplings 

9hvV'- Theoretical bounds from partial wave unitarity on these coefficients are rather weak. By us- 
ing the precision data for the LEP II Higgs search results and one-loop corrections, these coefficients 
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are constrained not so strongly as 

_4 < Im < 2 -6 < ^ < 5 

^ (A/TeV)2 < ^' < (A/TeV)2 < ^' 

-17 < ImM < 20 -5 < "^^^ < 6 (A 47) 

^ (A/TeV)2 < ^ < (A/TeV)2 < ^^'^^^ 

The measurement of the triple gauge boson couplings can also contribute to constrain these anoma- 
lous couplings, but the present results give weaker bounds. They are expected to be constrained 
more precisely at future collider experiments; for example, via WW — > WW [239] at the LHC, via 
e+e" W+W--f (ZZ-f) at the ILC [240], and 77 ^ ZZ at a 77 collider [241]. The effects of 
the anomalous couplings of this class modify the production and decay of the Higgs boson, so that 
they can be well tested via these processes at the LHC and the ILC after the Higgs boson is detected. 

[Class HI] {0^ „ 0^ 2) 

The last two are genuine dim-6 interactions of the Higgs boson, which are not severely constrained 
by the current precision electroweak data: see next subsection. 

4.12.2 Genuine Dimension Six Higgs Operators 

Recently, phenomenological consequences of the genuine dim-6 Higgs interactions have been stud- 
ied [242]. The operators ^ and 2 l^ad to a Lagrangian 

^' = E^O*,- (4-48) 

where ai = i v"^ / h?. This is the consistent formulation respecting the SM gauge symmetry. These 
operators contribute to the HVV and HHVV interactions and the Higgs boson self-interactions 
and H'^. In terms of the canonically normalized Higgs field and the physical Higgs boson mass, 
the interactions are expressed by 

^Hvv = (M^W;W-^ + \MlZ,Z^)({l-^)'^-^ + {l-a,)^\ (4.49) 
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= i^\{l- — + —^)H' ' , (4.50) 
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The measurement of HVV vertices can test ai while by measuring the Higgs self-coupling constants 
a2 (as well as ai) can be probed. For and interactions, the effect of 02 not only changes the 
Higgs potential, but also gives kinetic corrections which are enhanced at high energies. 

In [242], it has been shown that the parameter oi can be well measured at the ILC to an 
accuracy of 0.005 (0.003) corresponding to a scale A ~ 4 TeV, from 1 ab^^of data at 500 (800) GeV 
through the measurement of the production cross-sections from Higgs-strahlung and WW / Z Z- 
fusion for niH = 120 GeV. The parameter 02 modifies the form of the Higgs potential and thus the 
Higgs pair production cross-section. At the LHC the measurement of the HH production cross 
section is difficult because of the small rate with huge backgrounds [243]. At the ILC with the 
same integrated luminosity (1 ab~^), for mu = 120 GeV, 02 can be measured to 0.13 (0.07) at 500 
(800) GeV corresponding to a scale A « 1 TeV [79]; see Fig. HHl 



65 



Error 



500 GeV 

800 GeV - 0.25 



0.2 




0.15 



0.1 



0.05 



-0.4 



-0.2 







0.2 



0.4 







Figure 4.11: Combined statistical accuracy on 02 with an integrated luminosity of 1 ab for 
yfs = 500 GeV and 800 GeV for = 120 GeV, using the Higgs-strahlung and the VrVF-fusion 



4.12.3 Open questions 

In order to exhaust the possibilities of the physics of dimension six operators, the following items 
should be considered in future studies. 

• Fermionic operators should be included in the analysis. 

• It is also important to calculate the anomalous coupling in various new physics scenarios. 
By classifying models in terms of anomalous couplings, a systematic search for new physics 
beyond the SM is possible with the future precision data at the ILC, even when only a light 
Higgs boson is found at the LHC. 



Recently, a new class of models for electroweak symmetry breaking has been proposed, the "Little 
Higgs" models [244]. Similar to supersymmetry as a model-building principle, the Little-Higgs 
mechanism stabilizes the mass of scalar particles against its quadratic cut-off sensitivity. The Higgs 
boson appears as a member of the Goldstone boson multiplet pertaining to a global symmetry that 
is spontaneously broken at some scale A ~ 0(10 TeV). Naive dimensional analysis connects the 
cutoff A of such models with the decay constant of these bosons F by F ~ A/47r. This scale F is 
supposed to be above the electroweak scale v = 246 GeV, presumably in the TeV-region. 

Being Goldstone bosons, the scalars in Little-Higgs models have neither a mass nor a potential 
as long as the global symmetry is exact. However, the SM gauge and Yukawa interactions of the 
Higgs explicitly break the global symmetry, thereby radiatively generating a potential for the Higgs 
by the Coleman- Weinberg mechanism. This idea dates back to the 1970s [245], but at that time no 
mechanism was known that suppressed the one-loop corrections to the Higgs self-energy without 
fine-tuning among the couplings. 

This drawback is cured by the mechanism of collective symmetry breaking, originally observed 
in the context of deconstructed extra dimensions [246] . The structure of gauge and Yukawa inter- 
actions can be arranged such that if any one of their couplings is set to zero, the emerging global 



channels. 



4.13 Little Higgs 



66 



symmetry makes the Higgs an exact Goldstone boson. Only if all couplings are different from 
zero, a potential is generated. If the global symmetry is broken in this particular way, one-loop 
contributions to the Higgs potential depend only logarithmically on the cut-off, while a quadratic 
divergence appears only at two-loop order. This allows for a three-scale hierarchy without fine- 
tuning: V ~ F/Att ~ A/(47r)^. 

To fully implement collective symmetry breaking in all sectors and thus cancel all one-loop 
quadratic cutoff sensitivity, a complete set of new gauge bosons, scalars and fermions is needed. 
Just like the MSSM extends the SM by new weakly-interacting particles with a characteristic 
soft-SUSY-breaking mass scale, Little-Higgs models thus lead to an extended weakly-interacting 
effective theory, where most of the new particles are expected with masses of order F > 1 TeV. 

In order to successfully implement collective symmetry breaking, either the global symmetry 
representation must be reducible or the gauge group must be non-simple. Consequently, there are 
two basic lines of model biiilding. Particularly economic models emerge from a simple global symme- 
try group with an irreducible representation, such as the "Littlest Higgs" with its SU{5) S0{5) 
breaking chain [247] . The non-simple gauge structure of such models generically introduces W and 
7'/Z' recurrences of the SM vector bosons (for additional U(l) factors there always is the freedom 
not to gauge it, see the discussion at the end of this subsection). By contrast, simple- (gauge-) group 
models tend to have a proliferation of scalar states in multiple representations and typically lead to 
a two-Higgs-doublet model at low energies, while the extra gauge bosons corresponding to nondiag- 
onal generators of the extended EW gauge group are nontrivial to detect [248] . Finally, the "moose 
models" that were found first [249] establish a relation to higher-dimensional gauge theories [250] 
and feature factorized symmetries both in the scalar and vector sectors. 

Collective symmetry breaking in the fermion sector requires the introduction of extra fermions, 
most notably heavy top-like states that implement a see-saw mechanism for the top mass in con- 
junction with the cancellation of fermion loops in the Higgs potential. In simple-group models, the 
irreducibility of the gauge group calls for extra fermion partners in all three generations. 

By construction, Little-Higgs models involve tree-level mixing between the SM fields and their 
heavy partners. This leads to shifts in electroweak precision observables of the order v^/F^ which 
are at variance with data if F is below the TeV range. Detailed analyses [251] generically constrain 
the model parameters to _F ^ 2 — 3 TeV. An elegant way of removing some of the dangerous 
contributions is the implementation of a custodial SU{2) symmetry [252]. Incidentally, in some 
models the problem is also ameliorated if the Higgs boson is rather heavy while the hierarchy 
V F <^ A is kept [253]. The mixing of states is eliminated by the implementation of a discrete 
symmetry, T-parity [254]. Analogous to i?-parity in the MSSM and KK-parity in extra-dimension 
models, this leads to a one-loop suppression of shifts in precision observables and makes the lightest 
T-odd particle a suitable dark-matter candidate. Finally, embedding the little-Higgs mechanism in 
a supersymmetric model [255] can nicely remove the fine-tuning problems that are present both in 
the MSSM and in non-supersymmetric models. 

If the Little-Higgs mechanism is realized in nature, we need experimental means to verify this 
fact and to identify the particular model. If possible, we should observe the new particles, the 
symmetry structure in the couplings that is responsible for stabilizing the Higgs potential, and 
the non-linearity of the representation that is characteristic for Goldstone scalars at low energies. 
This can be achieved by a combined analysis of LHC and ILC data [251,256,257], while flavor 
data will play a minor role [258] . The prospects for uncovering the ultimate UV completion of the 
model [255,259] remain uncertain. 

The LHC provides us with a considerable detection reach for new vector bosons. If found, 
their properties (masses, mixing angles and branching ratios) can be determined with considerable 
accuracy. Simultaneously, new vector bosons can be studied by their indirect effect on contact 
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interactions at the ILC, where the search hmit extends into the 10 TeV range. New quarks can 
also be searched for at the LHC, although with less sensitivity, and they affect top-quark precision 
observables that can be measured at the ILC. 

For the scalar sector of Little Higgs models we have to consider a wide range of possibilities, 
essentially arbitrary combinations of singlets, doublets, triplets, and more. While such particles 
may be nontrivial to isolate in LHC data, at least the lighter part of the Higgs spectrum can be 
covered by searches at the ILC. Precision studies of vector-boson and Higgs anomalous couplings 
can even uncover traces of the nonlinear Goldstone-boson representation and thus establish the 
main ingredient of the Little-Higgs paradigm. 

Finally, a rather generic feature of Little-Higgs phenomenology originates from the enlarged 
global symmetry. This symmetry necessarily involves spontaneously broken U{1) subgroups. Con- 
sequently, the spectrum contains either Z' bosons (gauged U{1)) [251], or light pseudoscalar parti- 
cles otherwise [260]. The latter can be studied in diphotons at the LHC, in top-quark production 
at the ILC, and as resonances in 77 fusion at a photon collider. 

4.14 Models with Higgs Triplets 

Some of the "Little Higgs" models discussed in Sect. I4.1S1 postulate the existence of an SU(2)i 
triplet of scalar fields in addition to the conventional doublet. Unfortunately, the physical states 
associated with this triplet are expected to be too heavy for direct observation at a TeV linear 
collider. Nonetheless, "Little Higgs" models do not exhaust all possible models predicting triplets 
of scalar fields. These may occur, for example, in the models with extra spatial dimensions [261] 
(see also Sect. I4.15|) and in other scenarios. Even the Standard Model Higgs sector can be extended 
to contain both doublet and triplet fields without contradicting existing experimental data [262]. 
Models with Higgs triplets predict an enriched spectrum of physical states, including not only 
neutral and singly charged particles, but also doubly charged states. These can be directly produced 
and detected at the linear e+e~ collider. One shouldn't also exclude the possibility to detect a 
doubly charged Higgs boson, H , as a resonance at an e~e~ collider, though this will require 
an anomalously strong Hee coupling. The existence of scenarios in which the H state can be 
produced in e~e~ collisions with sufficiently high rate enabling its detection would emphasize the 
importance of a future ILC facility with several collision options. Clearly, rich phenomenology of 
the models with Higgs triplets invites investigation of the ILC potential for probing such a models. 

4.15 Higgs Bosons and Extra Dimensions 

In recent years, models with extra spatial dimensions have become very popular. They primarily 
address the electroweak hierarchy problem either by diluting gravity in a large extra dimensional 
volume that the SM interactions do not feel (as in the large or ADD extra dimensions [263]) or 
by explaining the apparent difference in scales by a warped metric in a small extra dimension (as 
in the warped or RS models [264]). Large extra dimensions necessarily have the Standard Model 
confined to a three dimensional sub-manifold, or brane, whereas RS exists in several versions which 
primarily differ as to how much of the SM is in the bulk. 

The 'classic' signatures involve the effects of the Kaluza-Klein (KK) tower of gravitons. In the 
large extra dimension case, the size of the extra dimensions is expected to be large enough that 
the KK states are extremely tightly spaced and thus appear as a continuum of objects, each one of 
which interacts only with ordinary gravitational strength. Relevant impact on collider observables 
occurs because there are so many states that they collectively have a large effect. In the case of 
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warped extra dimensions the states have masses with spacing on the order of TeV. However, each is 
much more strongly coupled than the ordinary graviton, and so individually they can have a large 
impact on processes at colliders. The upshot is that both theories can induce deviations of SM 
processes like e^e~ — > // and e~^e~ — > W^W~ from the virtual exchange of KK modes [265,266], 
or their real emission together with SM fermions or gauge bosons [267]. These modes have been 
studied with experimental models of detectors, see e.g. Ref. [8]. 

More recently, also the impact of extra dimensions on the Higgs boson phenomenology has 
been studied. In the ADD scenario, three effects have been analyzed (for another analysis, see 
Ref. [268]): 

1. A modification of the quasi-resonant W^W~ — > H production process through interference 
of the SM amplitude with the imaginary part of the graviton/graviscalar KK exchange am- 
plitude [269]. In order to yield a significant modification, a large total Higgs width is needed 
(i.e., large m//), which implies on the other hand a large center-of-mass energy. While the 
graviscalar contribution only modifies the normalization of the cross-section (by few percent 
for y/s = 1 TeV, itlh = 500 GeV and 2 extra dimensions at a fundamental Planck scale of 
1 TeV), a significant change of the angular distribution is expected from the spin-2 graviton 
exchange. 
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Figure 4.12: Effective Feynman diagrams from Equation (|4.52j) corresponding to single Higgs 
production with the subsequent Higgs decay into aa, where a is a heavy particle, such as W , Z or 
t [270]. 

2. A modification of the process e"'"e~ HHZ and the existence of the process e^e~ — > HH'j 
which is absent at tree level in the SM [271]. For a 1 TeV ILC and mn = 120 GeV, a sizable 
correction to e~^e~ HHZ both in normalization and angular distribution is expected for 
fundamental Planck scale up to a few TeV. Furthermore, the cross-section for e^e~ HH^ 
exceeds 0.1 fb for a fundamental Planck scale below approximately 2 TeV. In [271], expected 
5cr discovery limits on the fundamental Planck scale of 880-1560 (1640-2850) GeV have been 
derived at ^/s = 500 (1000) GeV for 6-3 extra dimensions. 

3. Direct Higgs pair production e+e^ HH [272,273]. The signature of this channel is two 
pairs of 5-jets with an invariant mass close to the Higgs mass. The main SM backgrounds 
are e~^e~ — > W~^W~, ZZ, ZH and bbbb. The differential cross section gives a characteristic 
signal of the spin 2 interaction. For an ILC operating at 1 TeV and a fundamental Planck 
scale of 2 TeV, the cross section is of the same order as the SM e~^e~ — > ZH for mn = 120 
GeV. Indeed, 700 Higgs pair events can be generated with an integrated luminosity of 500 
fb-i. 
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Figure 4.13: Feynman diagrams which contribute to the processes e~^e~ ZHH or e~^e~ 'jHH 
[271]. In this case, a five particle effective vertex generated by the virtual graviton exchange or 
branon radiative corrections has to be taken into account, which is also provided by the Lagrangian 




Figure 4.14: Feynman diagrams associated with the Higgs pair production e~^e~ HH [273]. 
The effective four particle vertex is provided by the exchange of virtual KK-graviton in the rigid 
ADD and by a branon loop for flexible branes (See Equation 1)4. 521) ). 



4.15.1 Radion Effects 

Extra dimensional models contain moduli which describe the size and shape of the extra dimensions. 
These scalar fields, usually called the radion, can be understood as a zero mass component of the 
higher dimensional metric, i.e., gxy where x and y are two of the compactified directions. In the case 
of a five dimensional model, there is a single scalar field which parameterizes the size of the extra 
dimension. In order to stabilize the size of the extra dimensions, these fields must acquire masses, 
and the specifics of the stabilization mechanism will in general play some role in determining the 
properties of the radion [274]. Generically, they are expected to couple to SM particles through 
the trace of the energy-momentum tensor, i.e., very similar to the SM Higgs boson, up to the trace 
anomaly of QCD. The lightest radion might in fact be lighter than the lightest graviton excitation 
and thus could be the discovery channel for the model. 
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In flat space, the radion is not usually relevant for colliders (see however Ref. [275]). The 
couplings are too small and not compensated by a large number of fields. They may nonetheless 
play an important role in cosmology and could even make up a fraction of the dark matter [276]. In 
warped backgrounds, the radion is much more relevant, because it picks up the same enhancement 
of its couplings as the graviton KK modes. Thus, the influence on the Higgs sector could be quite 
dramatic through mixing of the radion with the Higgs. This leads to mass eigenstates which are a 
mixture of both, and thus a strong modification of both the Higgs boson and radion properties, in 
particular their couplings to gauge bosons and fermions. For a review of the radion phenomenology, 
see e.g. [277]. Analyses of radion physics at the LHC and the ILC can be found in Ref. [12]. The 
radion sector is governed by 3 parameters: the strength of the radion-matter interactions described 
by an energy scale K^f), the mass of physical radion m^, and the radion-Higgs mixing parameter ^. 
In Fig. 14.151 the effective couplings squared of the Higgs boson and the radion (relative to those 
of a SM Higgs boson) are shown for the choice = 5 TeV, and three values of the radion mass 
(20, 55, 200 GeV) as a function of ^. Large deviations of the Higgs couplings from their SM values 
are expected if there is large radion-Higgs mixing present. The radion itself has couplings which 
are reduced by a factor v/ with respect to those of a SM Higgs in the case of no mixing, which 
requires high luminosity for direct discovery. The sensitivity of the trilinear Higgs coupling to 
radion admixtures has been studied as well in [277]. 
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Figure 4.15: Effective coupling of the Higgs boson (left) and the radion (right) to Z boson (from 
[277]). 



4.15.2 Branon Effects 

The previous discussion has assumed that the branes are perfectly rigid. However, it is expected 
that the dynamics responsible for generating the brane will in fact leave it with a finite tension 
which will allow it to oscillate. The brane fiuctuations are manifest as new pseudo-scalar particles, 
which are called branons [278]. For appropriate choices of parameters, branons are natural dark 
matter candidates, being stable and weakly coupled to Standard Model fields [279]. The collider 
signature for branons is missing energy; in particular, the mono-jet and single photon signals for 
hadron colliders [280] and the single Z and single photon channels for e^e~ colliders [281]. This last 
signature has been studied experimentally by LEP [282], finding the most constraining bound for 
the brane tension scale: / > 180 GeV, for light branons. These collider searches are complementary 
to astrophysical and cosmological ones [283] and they need to be completed with branon radiative 
correction analyses. Indeed, these quantum corrections are very similar to the virtual KK-graviton 
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effects in the rigid ADD model since the two effects predict the same effective Lagrangian [284]: 

A£efr = WlT^,T^'- + W2TliT^ , (4.52) 

where Wi and W-2 are two constants which depends on the renormahzation procedure. This means 
that branons lead to analogous four body interactions and Higgs boson phenomenology. 
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Chapter 5 

Analysis Tools 



In this section we give a brief overview about the state-of-the-art tools that are available for the 
evaluation of Higgs boson observables. The first subsection summarizes the computer codes that 
have been developed for the calculation of Higgs boson mass spectra and branching fractions. The 
second subsection briefly describes event generators available for Higgs boson physics. All codes 
are listed in Tab. 15.11 together with their corresponding web address (see also the Les Houches Web 
Repository for BSM Tools [285]). 



HDECAY people . web . psi . ch/ spira/hdecay/ 

FeynHiggs www.feynhiggs.de/ 

CPsuperH www . hep . man . ac . uk/ u/ j slee/CPsuperH . html 

NMHDECAY www . th . u-psud . f r/NMHDECAY/nmhdecay . html 

HERWIG hepwww . rl . ac . uk/ theory/ seymour/herwig/ 

ISA JET www.phy.bnl.gov/~isajet/ 

MADGRAPH madgraph . hep . uiuc . edu/ 

PYTHIA www . thep . lu. se/~torbjorn/Pythia.html 

SHERPA www . physik . tu-dresden . de/~krauss/hep/ 

WHIZARD www-ttp . physik . uni-karlsruhe . de/ whizard/ 



Table 5.1: Web links for Higgs-related tools. 



5.1 Higgs spectra and branching fractions 
5.1.1 Standard Model Higgs 

The Standard Model Higgs decay branching fractions and total width are computed by the packages 
HDECAY [38] and FeynHiggs [37,148,286]. We outline here the known radiative corrections and 
sources of theoretical uncertainty in the SM Higgs branching ratios. We focus here on HDECAY: for 
some of the decays to bb, cc, ss, we believe that HDECAY is slightly better than FeynHiggs in terms 
of containing SM higher order radiative corrections to the decay widths; however, the other decays 
are treated to comparable accuracy. 

The radiative corrections to Higgs decays to fermion and boson pairs have been reviewed in 
Ref. [287]; we give here a brief sketch of the known corrections and refer to Ref. [287] for references 
to the original calculations. The full QCD corrections to the Higgs decay to qq are known up to 
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three loops neglecting the quark mass in the kinematics; at the two loop level they are known for 
massive final-state quarks in the leading-A'^p approximation. The electroweak corrections to the 
Higgs decay to quark or lepton pairs are known at one-loop; in addition, the QCD corrections 
to the leading top-mass-enhanced electroweak correction term are known up to three loops, to 
order Gi?m|a^. All of these corrections to the Higgs partial widths to fermions are included in a 
consistent way in the program HDECAY [38]. 

For the Higgs masses below the WW threshold, Mh ^ 160 GeV, decays into off-shell gauge 
bosons {WW, ZZ) are important and affect the total Higgs width, thus feeding in to the Higgs 
branching fractions. Off-shell decays to W and Z bosons are taken into account. One-loop elec- 
troweak corrections to Higgs decays to WW and ZZ arc known, together with the QCD corrections 
to the leading 0{GFmf) result up to three loops. These corrections to Tw^z amount to less than 
about 5% in the intermediate Higgs mass range [287] (translating to less than roughly 2% in 
BR(i7 hb) for Mh — 120 GeV) and have been neglected in HDECAY, although their inclusion 
would seem straightforward. 

The Higgs decay into gluon pairs, H — > gg, is of order 0{GFal) leading one-loop order 

and therefore suffers a large scale dependence. In the SM, the state of the art is the 0{GFOt^) 
contribution in the large mt limit computed in Ref. [288]. This correction is of order 20% for a 
light Higgs, increasing the total Higgs hadronic width by only about 1%, and the remaining scale 
uncertainty from varying the renormalization scale between Mj^/2 and 2Mh is roughly ±5%. This 
correction appears not to be included in HDECAY, although its inclusion would be trivial. 

The SM H^jj decay partial width receives QCD corrections, which of course only affect the 
quark loop diagrams. Because the external particles in the 'j'jH vertex are color neutral, the virtual 
QCD corrections are finite by themselves. Since no real radiation diagrams contribute, the QCD 
corrections to — 77 are equivalent to those to the inverse process 77 — H. This is in contrast 
to the QCD corrections to the ggH vertex. The QCD corrections to F-y in the SM are known 
analytically at the two- loop [OlGpas)] order [289] and as a power expansion up to third order in 
Mn/mt at three-loop [©(G^a^)] order [290]. They are small for Higgs masses Mh < '^nit, the 
0{as) corrections are only of order 2% for Mh < 2M]y, and the 0(0^) corrections are neghgible, 
demonstrating that the QCD corrections are well under control. The SM — >^ 77 decay partial 
width also receives electroweak radiative corrections. The electroweak corrections are much more 
difficult to compute than the QCD corrections and a full two-loop calculation does not yet exist. 
The electroweak correction due to two-loop diagrams containing light fermion loops and W or Z 
bosons (with the Higgs boson coupled to the W or Z boson, because the light fermion Yukawa 
couplings are neglected) was computed recently in Ref. [291] and contributes between —1% and 
—2% for Mh ^ 140 GeV. The leading 0{GFm^) electroweak correction due to top- mass-enhanced 
two-loop diagrams containing third-generation quarks was also computed recently in Ref. [292] as 
an expansion to fourth order in the ratio Mfj / {2M\y)^ The expansion appears to be under good 
control for Mh ^ 140 GeV, where this correction contributes about —2.5% almost independent of 
Mh- The leading 0{GfM'^) correction was computed in Ref. [294] for large Mh', however, this 
limit is not useful for a light Higgs boson. We conclude that the electroweak radiative corrections 
to — >^ 77 appear to be under control at the 1-2% level. 

The Higgs branching ratios in the SM also have a parametric uncertainty due to the experimental 
uncertainties in the SM input parameters. The largest sources of parametric uncertainty are the 
bottom and charm quark masses and (to a lesser extent) the strong coupling , which contributes 
via the QCD corrections to the Hqq couplings. These parametric uncertainties were evaluated in, 
e.g., Refs. [295,296]. Ref. [296] found a parametric uncertainty in BR{H bb) of about 1.4% 

^The 0{GF'mt) electroweak correction was also considered in Ref. [293], whose results disagree with that of 
Ref. [292]. The source of this disagreement is addressed in Ref. [292]. 
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for Mh = 120 GeV, using the standard ag = 0.1185 it 0.0020 [297] and a somewhat optimistic 
Tnb{mb) = 4.17 ± 0.05 GeV (MS) [298]. The parametric uncertainty in BR(// bb) is suppressed 
due to the fact that makes up about 2/3 of the Higgs total width at Mh = 120 GeV, leading to 
a partial cancellation of the uncertainty in the branching ratio; the uncertainty instead feeds into 
the other Higgs branching ratios to W^W~, tt, etc. We also expect the parametric uncertainty 
in BR(i? bb) to be somewhat larger at higher Higgs masses, where no longer dominates the 
total width. 

The best measurements of as come from LEP-I and II; the Tevatron and LHC are unlikely 
to improve on this. The ILC is expected to improve the precision of as by about a factor of 
two [53,299]. The bottom quark mass is extracted from heavy quarkonium spectroscopy and B 
meson decays with a precision limited by theoretical uncertainty. There are prospects to improve 
the bottom quark mass extraction through better perturbative and lattice calculations [300] and 
more precise measurements of the upsilon meson properties from CLEO [301]. 

5.1.2 MSSM Higgs 
CP-conserving MSSM 

The three major codes used to compute the Higgs spectrum and branching fractions in the CP- 
conserving MSSM are HDECAY [38], FeynHiggs [37,148], and CPsuperH [108]. 

HDECAY computes the MSSM Higgs spectrum and mixing angles with the user's choice of pack- 
ages: Subh [302] and Subhpole [303] (now subsumed into CPsuperH [108,143,304]) as well as Haber 
et al. [305] , which use a renormalization-group-improved effective field theory approach to compute 
Higgs masses and mixing angles; and FeynHiggsFast^ [37], which contains the two-loop diagram- 
matic computation (the most recent version uses FeynHiggsFast by default). A reconciliation 
between the two-loop diagrammatic computation of the mass of the lightest CP-even Higgs boson 
and the effective field theory computation was performed in Ref. [306]. The code contains, in ad- 
dition to the SM modes, Higgs decays to SUSY particles and SUSY contributions to loop- induced 
decays. In the case of Higgs decays to SM particles through tree-level couplings, the MSSM rates 
are computed by scaling the SM couplings by the appropriate Higgs mixing angles. 

FeynHiggs, in addition to computing the MSSM Higgs mass spectrum and mixing angles using 
the on-shell Feynman-diagrammatic approach at the two-loop level, also computes the Higgs decay 
branching fractions and the Higgs boson production cross sections at the Tevatron and the LHC for 
all relevant channels (in a simplified parametric approach). In some of the SUSY decays, FeynHiggs 
includes corrections that are not in HDECAY, e.g., the A;, corrections in the charged Higgs decay to 
tb. Details and references were given in Sec. 14.11 

CPsuperH contains a computation of the MSSM Higgs boson pole masses and mixing angles 
based on the renormalization-group-improved effective potential computation in the MSSM [304]. 
Since the case with real parameters is only a special case of the complex MSSM, more details about 
the two codes will be discussed further in the next subsection. 

CP-violating MSSM 

The two major codes used to compute the Higgs spectrum and branching fractions in the MSSM 
with explicit CP violation are FeynHiggs [37,145,147,148] and CPsuperH [108]. 

The program FeynHiggs [37,148,286] is based on the results obtained in the Feynman-diagrammatic 
(FD) approach [16,17,129,137]. It computes the MSSM Higgs spectrum and mixing angles, taking 
into account the full phase dependence at the one-loop and a partial phase dependence at the 

■^FeynHiggsFast, though still available, has been superseded by FeynHiggs2 . 2. 
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two-loop level. The code CPsuperH [108] is based on the rcnormalization group (RG) improved 
effective potential approach [303,306-308]. For the MSSM with real parameters the two codes can 
differ by up to ~ 4 GeV for the light CP-even Higgs boson mass, mostly due to formally subleading 
two-loop corrections that are included only in FeynHiggs. For the MSSM with complex parame- 
ters the phase dependence at the two-loop level is included in a more advanced way [143, 309] in 
CPsuperH, but, on the other hand, CPsuperH does not contain all the subleading one-loop contri- 
butions that are included [145-147] in FeynHiggs. In the case of real parameters, the remaining 
theoretical uncertainty on the light CP-even Higgs boson mass has been estimated to be below 
^ 3 GeV [17,52], if FeynHiggs is used. In the case of complex parameters, no detailed estimate 
for the additional theory uncertainty due to the complex phases has been performed. However, the 
intrinsic uncertainties should be larger than for real parameters by at most a factor of two. 

5.1.3 Beyond the MSSM 
The NMSSM 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) contains a Higgs sector enlarged 
by a singlet chiral superfield compared to the MSSM. The Higgs masses, couplings and decay 
widths in this model are computed by the Fortran code NMHDECAY [181,310]. The computation 
of the Higgs spectrum includes leading electroweak corrections and leading two-loop terms. The 
computation of the decay widths is carried out as in HDECAY [38] (temporarily without three-body 
decays); Higgs-to-Higgs decays are included and can be important. Each point in parameter space 
can be checked against low-energy observables, as well as negative Higgs boson searches at LEP and 
negative sparticlc searches at LEP and the Tevatron, including unconventional channels relevant 
for the NMSSM. Version 2.0 [310] allows to compute the dark matter relic density via a link to a 
NMSSM version of the MicrOmegas code [311]. 

Other extended models 

There are many extended Higgs models in the literature for which no public codes exist. The Higgs 
couplings have been studied in detail using private codes for some of these models, as follows: 

• Randall-Sundrum model: mixed radion and Higgs [312,313], 

• Randall-Sundrum model: Higgs with a profile in the warped extra dimension [314], 

• Universal extra dimensions and Kaluza-Klein effects [315], 

• Littlest Higgs model [256,257,316,317]. 

5.2 Event generators 

Here we give brief descriptions of Monte Carlo event generators that produce Higgs events for 
the ILC. Many of these rest upon the spectrum and branching fraction calculators like FeynHiggs 
and/or HDECAY discussed in the previous sections.^ 

• HERWIG 

^Tho depth of the descriptions does not reflect the importance of the respective generator, but the information 
provided to the authors. 
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HERWIG [318—322] is a general-purpose Monte Carlo (MC) event generator for high-energy 
collider processes within the SM and MSSM [323,324]. The HERWIG source codes and related 
information can be found in [325]. 

Higgs production and decay is available within the SM and MSSM. While the SM imple- 
mentation is built-in, the MSSM generation relies on input information provided by other 
programs. As the HERWIG generator does not contain any built-in models for SUSY-breaking 
scenarios, in all cases the general MSSM particle spectrum and decay tables must be pro- 
vided just like those for any other object, with the following caveats: (i) (coloured) SUSY 
particles do not radiate (which is reasonable if their decay lifetimes are much shorter than 
the QCD confinement scale); (ii) CP-violating SUSY phases are not included. A package, 
ISAWIG, has been created to work with ISAJET [326] to produce a file containing the SUSY 
(and Higgs) particle masses, lifetimes, couplings and mixing parameters. This package takes 
the outputs of the ISAJET MSSM programs and produces a data file in a format that can 
be read into HERWIG for the subsequent process generation. Of particular relevance to Higgs 
production and decay is also the interface to the HDECAY program [327]. As of 2005, also a link 
to FeynHiggs is available. However, the user can produce her/his own file provided that the 
correct format is used; see the corresponding web page, where examples of input files can be 
found (the SUSY benchmark points recommended in Rcf. [117] are also uploaded). Finally, 
the implementation of the SUSY Les Houches accord of Ref. [328] will soon be available. 

The SM and MSSM Higgs production and decay processes of relevance for the physics of 
an ILC [329] currently implemented in HERWIG are are associated with the process numbers 
IPROC = 300,400,900, 1000 and 1100 (i.e., Higgs-strahlung, vector-boson fusion, pair produc- 
tion and associated production with heavy fermions). Also the decay of SUSY particles to 
Higgs bosons is considered [323,324]. 

• ISAJET 7.73 

The event generator Isajet [326] can be used to generate SM Higgs boson production via 
the process e~^e~ ZHsm using the E+E- reaction type, and stipulating the final state to 
be 'ZO' and 'HIGGS'. The Higgs mass is input using the HMASS keyword. The Hsm decay 
modes are calculated internally at tree level. The beam polarization can be adjusted using 
the EPOL keyword, which allows input of Pl^g") and/or P£,(e"'"). The events can be generated 
with initial state photon bremsstrahlung and/or beamstrahlung using the EEBREM or EEBEAM 
keywords. Input of the beamstrahlung parameter T and bunch length is needed. Then 
an eflFective electron and positron distribution function is calculated, and the hard scattering 
events are convolved with the PDFs, to give a variable CM energy for the hard scattering. 
Fox- Wolfram QCD showers and independent hadronization are used to convert final state 
quarks and gluons into jets of hadrons. 

Isajet can also generate Higgs production in the MSSM. One may either run with weak scale 
MSSM input parameters, or with high scale parameters in supergravity models, GMSB models 
or AMSB models. The full SUSY and Higgs mass spectrum is calculated using 2-loop RGEs to 
obtain an iterative solution, and full one-loop radiative corrections are included. The various 
Higgs masses are computed using the RG improved one-loop effective potential evaluated at 
an optimized scale to account for leading two-loop terms. Agreement for the Higgs boson 
masses with FeynHiggs is in the few GeV range, typically. The production reactions include 
e+e" Zh, ZH, Ah, AH and H+R- (along with the various SM and MSSM 2^2 
processes). The sparticle mass spectrum and decay table can be output separately using the 
isasugra or isasusy subprograms. Leading QCD corrections are included for Higgs decay 
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to heavy quarks. It should be noted that the various MSSM Higgs bosons can be produced at 
large rates via sparticle cascade decays. Also, the Isajet decay table includes Higgs boson 
decays to sparticles, which can sometimes modify and even improve prospects for MSSM 
Higgs boson discovery. 

As with SM Higgs production, beam polarization and brems/beamstrahlung effects can be 
included in the event generation. While spin correlation effects are neglected, full decay 
matrix elements are used for three-body sparticle decays. Decays to left and right tau lepton 
helicity states are computed, to give on average the correct energy distributions from decays 
to tau leptons. 

• MadGraph 

MadGraph [330] is an automated tree-level matrix element generator based on helicity ampli- 
tudes as implemented in the HELAS library [331]. The associated event generator MadEvent 
[332] can then be used to generate unweighted events. For showering/hadronization, interfaces 
to PYTHIA and HERWIG are available. 

Higgs production in both the SM and MSSM [333] is handled at tree-level. In the MSSM, the 
Higgs spectrum and mixing angles are supplied from a SUSY Les Houches Accord (SLHA) 
[328] input file, which can be generated with the user's choice of spectrum calculator. Higgs 
decay widths in both the SM and MSSM are taken from HDECAY, which is integrated into the 
HELAS library. 

• PYTHIA 

PYTHIA [111,334] is a general-purpose event generator for hadronic events in e^e~, eh, and 
hh collisions (where h is any hadron or photon). For a recent brief introduction to this and 
many other codes, see [335]. The current version is always available from the PYTHIA web 
page with update notes and a number of useful examples. 

PYTHIA contains a large number of Higgs boson processes, both for production (in many cases 
with polarisation effects included, see [111], Section 8.8) and decay. The partial widths of the 
Higgs bosons are very strongly dependent on the mass, and are calculated as a function of 
the actual Higgs mass, i.e., not just at the nominal mass. Leading order expressions are used, 
with running fermion masses and, in some cases, overall if-factors, to absorb the leading 
radiative corrections. Other higher-order effects are not included. Since the Higgs is a spin-0 
particle, it decays isotropically. In decay processes such as W^W^ /Z^Z^ 4 fermions, 

angular correlations are included [336]. In decays to 2 jets, the parton shower off the final 
state is merged with the 3-jet matrix elements for hard jet emission [337], thus incorporating 
part of the NLO correction. 

The relevant physics scenarios included in PYTHIA can be categorised as follows: 

— Light Standard Model Higgs (narrow- width approximation): In e'^e~ annihilation, the 
main processes are e+e~ — > h^Z^, usually dominant close to threshold, and vector boson 
(Z^Z^ and W^W^) fusion which is important at high energies. The vector boson fusion 
process is a full 2 — > 3 body calculation. The loop-induced 77 fusion process may also be 
of interest, in particular when the effects of beamstrahlung photons and backscattered 
photons are included. For details, see [111], Section 8.5.1. 

— Heavy Standard Model Higgs: for details, see [111], Section 8.5.2. 

— SUSY/2HDM: the Higgs sector in Supcrsymmctric models is a special case of the general 
2-Higgs Doublet Model included in PYTHIA. The internal PYTHIA routines for fixing the 
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MSSM Higgs sector rely on the efFeetive potential approach of [302,303], including the 
important radiative corrections at large tan /3. Other Higgs boson spectrum calculations 
can be interfaced, either via the SUSY Les Houches Accord [328,338], or via the run- 
time interfaces to both Isasusy [339] and FeynHiggs [37,148]. Recently, an interface 
to NMSSM models via the Les Houches Accords [113,328] has also been implemented 
[340]. The additional production processes possible in a 2HDM are included, such as 
Higgs boson pair production. In production and decay, however, the effects of virtual 
superpartners are not included. For details, see [111], Sections 8.5.3-8.5.5 (2HDM) and 
8.7 (SUSY). 

^ Left-Right symmetry / Higgs Triplets: the particle content is expanded by right-handed 
and and right-handed neutrinos. The Higgs fields are in a triplet representation, 
leading to doubly-charged Higgs particles, one set for each of the two SU(2) groups. Also 
the number of neutral and singly-charged Higgs states is increased relative to the Stan- 
dard Model. PYTHIA implements the scenario of [341,342]. The main decay modes imple- 
mented are H^^ W^W^jfl^ generation indices) and ^ W^^^^ifi^. 
For details, see [111], Section 8.6.3 

— Technicolor: while not a model with a fundamental Higgs boson. Technicolor still belongs 
to the category of Higgs-related topics. Techni-rho production is possible using the 
Heavy Standard Model Higgs processes or through a vector dominance mechanism [343] ; 
see [111], Section 8.6.7. 

— Little Higgs: Standard PYTHIA does not explicitly include a Little Higgs model, but 
important features can be handled using the options for 4th generation fermions and 
possibly Z'/W' bosons. For example, single heavy top quark T production and decay 
can be simulated using the t'. See [111], Sections 8.6.1 and 8.6.2. 

• SHERPA 

SHERPA [344] is a full event generator which is capable of describing events in e+e~, pp, pp, 
and 77 collisions; for the latter the photon spectra produced through laser-backscattering are 
parametrised according to CompAZ. For the description of the perturbative aspects of such 
events, it includes a fully automated matrix element generator for the calculation of cross 
sections at the tree-level (AMEGIC++ [345]), a module for the simulation of subsequent soft 
QCD radiation through the parton shower (APACIC++ [346]) and a merging of both according 
to the prescription of [347-349]. 

All production and decay processes in the SM and the MSSM arc treated at tree-level, in- 
cluding processes involving Higgs bosons. In the case of the MSSM, the particle spectra and 
mixing angles are taken from a SUSY Les Houches Accord (SLHA) [328] input file. Additional 
interfaces to IsaSUSY, HDECAY and FeynHiggs are available. 

• WHIZARD 

WHIZARD [350] is a program system designed for the efficient calculation of multi-particle 
scattering cross sections and simulated event samples. Tree- level matrix elements are gener- 
ated automatically for arbitrary partonic processes by calling one of the external programs 
O'Mega [351], MADGRAPH or CompHEP [109]. Matrix elements obtained by alternative methods 
(e.g., including loop corrections) may be interfaced as well. The program is able to calcu- 
late numerically stable signal and background cross sections and generate unweighted event 
samples with reasonable efficiency for processes with up to six final-state particles. 



79 



Polarization is treated exactly for both the initial and final states. Final-state quark or lepton 
flavors can be summed over automatically where needed. For ILC physics, bcamstrahlung 
(CIRCE) and initial-state radiation spectra are included for electrons and photons. Exact color 
flow information is maintained, which is used when fragmenting and hadronizing events with 
the built-in PYTHIA interface. The events can be written to file in STDHEP or ASCII format. 

Currently, WHIZARD supports the Standard Model and the MSSM, optionally with anomalous 
couplings, but model extensions or completely different models can be added. Higgs pro- 
duction and decay couplings, and spectra in the case of the MSSM, are adopted from the 
corresponding matrix element generator. 
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Chapter 6 



Interplay between LHC and ILC in 
the Higgs Sector 

If a Higgs-like state is discovered at the LHC and the ILC, independent of the reahzation of the 
Higgs mechanism, important synergistic effects arise from the interplay of LHC and ILC [12]. 

6.1 Overview 

If a state resembhng a Higgs boson is detected, it is crucial to experimentally test its nature as a 
Higgs boson. To this end the couplings of the new state to as many particles as possible must be 
precisely determined, which requires observation of the candidate Higgs boson in several different 
production and decay channels. Furthermore the spin and the other CP-properties of the new 
state need to be measured, and it must be clarified whether there is more than one Higgs state. 
The LHC will be able to address some of these questions, but in order to make further progress a 
comprehensive program of precision Higgs measurements at the ILC will be necessary. While the 
ILC will provide a wealth of precise experimental information on a light Higgs boson, the LHC 
may be able to detect heavy Higgs bosons which lie outside the kinematic reach of the ILC (it is 
also possible, however, that the ILC will detect a heavy Higgs boson that is not experimentally 
observable at the LHC due to overwhelming backgrounds). Even in the case where only one scalar 
state is accessible at both colliders, important synergistic effects arise from the interplay of LHC 
and ILC. This has been demonstrated for the example of the Yukawa coupling of the Higgs boson 
to a pair of top quarks, see below. 

The LHC and the ILC can successfully work together in determining the CP properties of the 
Higgs bosons. In an extended Higgs sector with CP violation there is a non-trivial mixing between 
all neutral Higgs states. Different measurements at the LHC and the ILC (in both the electron- 
positron and the photon-photon collider options) have sensitivity to different coupling parameters. 
In the decoupling limit, the lightest Higgs boson is an almost pure CP-even state, while the heavier 
Higgs states may contain large CP-even and CP-odd components. Also in this case, high-precision 
measurements of the properties of the light Higgs boson at the ILC may reveal small deviations 
from the SM case, while the heavy Higgs bosons might only be accessible at the LHC. In many 
scenarios, for instance the MSSM, CP- violating effects are induced via loop corrections. The CP 
properties therefore depend on the particle spectrum. The interplay of precision measurements in 
the Higgs sector from the ILC and information on the SUSY spectrum from the LHC can therefore 
be important for revealing the CP structure. As an example, if CP-violating effects in the Higgs 
sector in a SUSY scenario are established at the ILC, one would expect CP-violating couplings in 
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the scalar top and bottom sector. The experimental strategy at the LHC could therefore focus on 
the CP properties of scalar tops and bottoms. 

While the most studied Higgs boson models are the SM and the MSSM, more exotic realisations 
of the Higgs sector cannot be ruled out. Thus, it is important to explore the extent to which search 
strategies need to be altered in such a case. 

A possible scenario giving rise to non-standard properties of the Higgs sector is the presence 
of large extra dimensions, motivated for instance by a "fine-tuning" and "little hierarchy" problem 
of the MSSM, see Sect. 14.151 A popular class of such models comprise those in which some or all 
of the SM particles live on 3-branes in the extra dimensions. Such models inevitably require the 
existence of a radion (the quantum degree of freedom associated with fluctuations of the distance 
between the 3-branes or the size of the extra dimension(s)). The radion has the same quantum 
numbers as the Higgs boson, and in general the two will mix. Since the radion has couplings 
that are very different from those of the SM Higgs boson, the two physical eigenstates will have 
unusual properties corresponding to a mixture of the Higgs and radion properties; the prospects 
for detecting them at the LHC and ILC must be carefully analysed. In cases of this kind the ILC 
can observe both the Higgs boson and the radion, and covers most of the parameter space where 
detection of either state at the LHC is difficult. 

A case where the LHC detects a heavy (500 GeV-1 TeV) SM-like Higgs boson rather than a 
light CP-even Higgs boson, as apparently needed to satisfy precision electroweak constraints, can 
also occur in a general two-Higgs-doublet model. The source of the extra contributions mimicking 
the effect of a light Higgs boson in the electroweak precision tests may remain obscure in this case. 
The significant improvement in the accuracy of the electroweak precision observables obtainable at 
the ILC running in the GigaZ mode and at the WW threshold will be crucial to narrow down the 
possible scenarios. 

Another challenging SUSY scenario is the MSSM with an extra singlet (NMSSM, see Sect. 14. f)) . 
especially if the extra singlet dominates the lightest Higgs. While such a state has reasonably large 
production cross sections at the LHC, it would be difficult to detect as it mainly decays hadronically. 
Such a state could be discovered at the ILC. From the measurement of its properties, the masses 
of the heavier Higgs bosons could be predicted, guiding in this way the searches at the LHC. For a 
very heavy singlet-dominated Higgs state, on the other hand, the kinematic reach of the LHC will 
be crucial in order to verify that a non-minimal Higgs sector is realised. Thus, input from both the 
LHC and the ILC will be needed in order to provide complete coverage of the NMSSM parameter 
space. 

If no light Higgs boson exists, quasi-elastic scattering processes of W and Z bosons at high 
energies provide a direct probe of the dynamics of electroweak symmetry breaking. The amplitudes 
can be measured in 6-fermion processes both at the LHC and the ILC. The two colliders are sensitive 
to different scattering channels and yield complementary information. The combination of LHC 
and ILC data will considerably increase the LHC resolving power. In the low-energy range it will 
be possible to measure anomalous triple gauge couplings down to the natural value of l/lGvr^. The 
high-energy region where resonances may appear can be accessed at the LHC only. The ILC, on 
the other hand, has an indirect sensitivity to the effects of heavy resonances even in excess of the 
direct search reach of the LHC. Detailed measurements of cross sections and angular distributions 
at the ILC will be crucial for making full use of the LHC data. 
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6.2 Higgs coupling determination 



6.2.1 Combined LHC/ILC analysis 

At the LHC, without any theory assumption, no coupling measurement of a Higgs boson can be 
performed. With "mild" theory assumptions one might hope for a ttH coupling measurement at 
the level of 12-19% [7]. At the ILC, on the other hand, the ttH coupling can be measured only at 
high energies, ~ 800 GcV. However, even in the first stage of the ILC, the -y/s = 500 GcV data 
can be combined with the LHC data on Higgs branching ratios to obtain information about the 
ttH coupling (see Ref. [352] for a first evaluation). Using information from all Higgs production 
and decay modes at the LHC [7] and from the ILC, i.e., 

• Mh, 

• crtot(e+e~ HZ), 

• o-tot(e+e- ^ HZ) x BR{H X) {X = bb, t+t', gg, WW*), and 

• (7tot(e+e- uvH) x BR(i7 bb), 

a determination of g^ifj down to the level of 11-14% can be performed without any additional 
assumptions on the underlying theory (for ruH ^ 200 GeV, and assuming SM-like BR measure- 
ments) [353]. 

In a similar way the coTipIing of the Higgs boson to photons can be constrained in a better way 
than with the LHC or the ILC results alone. In the combined analysis of LHC and ILC data [353] 
(with the ILC running at = 500 GeV and without the 77 collider option) the ^'^H coupling can 
be determined to the level of 7-8% for niH ^ 200 GeV (and assuming SM-like BR measurements). 

6.2.2 Model-Independent Analysis for the Higgs Boson Couplings 

Once the (lightest) Higgs boson is found at the LHC, its mass, production cross section, and decay 
branching ratios will be measured as precisely as possible at the LHC and the ILC in order to 
determine the Higgs boson couplings to gauge bosons and fermions. The predictions for these 
couplings differ among various models beyond the SM. By determining these couplings accurately, 
these models may be distinguished, even when no further new particle is found. 

These Higgs boson couplings may be written in a model-independent way as [10,35], 

£ = x—hbb + y(—htt + —hcc) + z—hrr + uigniwhWaW^ + ^nizhZ^Z''), 
V V V V 2 

where the four parameters x, y, z, and u represent the multiplicative factors in the Higgs-boson 
coupling constants with down-type quarks, up-type quarks, charged-leptons and massive gauge 
bosons, respectively. The SM corresponds tox = y = z = u=l. This expression is valid at 
lowest order in the SM, MSSM, NMSSM and multi Higgs doublet models without tree-level flavor 
mixing. The expected accuracy of the parameter determination has been evaluated with respect 
to the SM case. It has been shown that the u and x parameters are determined to the few-percent 
level, and y and z are constrained to less than 10%. Prom the correlation of the four parameters 
determined at the ILC, it may also be possible to distinguish various models. For example, the 
Type-I two-Higgs-doublet model, which predicts the relation x=y, and the MSSM lie in different 
regions of the x-y space. For a large tan /? value, the allowed range of the x-z space can deviate 
from the x=z line for the MSSM because of the SUSY corrections to the hbb vertex [107,354]. 
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6.3 What if Only a Light Higgs and Nothing Else is Found at the 
LHC? 



If (only) a state resembling a Higgs boson is detected, it is crucial to experimentally test its nature 
as a Higgs boson. The couplings of the supposed Higgs boson to as many particles as possible must 
measured; see also the previous subsection. Furthermore the spin and the other CP-properties of 
the new state need to be measured. The LHC will be able to address some of these questions, but 
in order to make further progress a comprehensive program of precision Higgs measurements at 
the ILC will be necessary. The significance of the precision Higgs program is particularly evident 
from the fact that many extended Higgs theories over a wide part of their parameter space have 
a lightest Higgs scalar with nearly identical properties to those of the SM Higgs boson. In this 
decoupling limit additional states of the Higgs sector are heavy and may be difficult to detect both 
at the LHC and ILC. 

6.3.1 Indirect Determination of the Heavy Higgs Boson Masses in the MSSM 

In the MSSM, information on the heavy Higgs boson masses can be extracted from the experimental 
information on branching ratios of the lightest Higgs boson [355] (see also Ref. [356] for an earlier 
study) . 

The analysis has been performed in an SPS la-based scenario [117], where Ma has been kept as 
a free parameter, but so heavy that the LHC only detects one light Higgs boson. For the parameters 
of the SPS la scenario this corresponds to the region Ma ^ 400 GeV. 

The precise measurements of Higgs branching ratios at the ILC together with accurate deter- 
minations of (parts of) the SUSY spectrum at the LHC and the ILC (see Ref. [12]) will allow in 
this case to obtain indirect information on Ma (for a discussion of indirect constraints on Ma from 
electroweak precision observables, see Ref. [53]). When investigating the sensitivity to Ma it is cru- 
cial to take into account realistic experimental errors on the other SUSY parameters that enter the 
prediction of the Higgs branching ratios. In Ref. [355] all the SUSY parameters have been varied 
according to error estimates for the measurements at LHC and ILC in this scenario. The sbottom 
masses and the gluino mass can be obtained from mass reconstructions at the LHC with ILC input; 
see Ref. [12]; the precisions have been assumed to be Arrig = ±8 GeV and Am^^ ^ « ±7.5 GeV. 
The lighter stop (which in the SPS la scenario has a mass of about 400 GeV, see Ref. [117]) will be 
accessible at the ILC, leading to an accuracy of about Am^^ = ±2 GeV. The the stop mixing angle. 
Of:, can also be determined with high accuracy [355]. For tan /3 an uncertainty of A tan (3 = 10% has 
been used (this accuracy can be expected from measurements at the ILC in the gaugino sector for 
the SPS la value of tan/3 = 10 [357]). Concerning the top quark mass an error of Amj = ±0.1 GeV 
has been assumed from the ILC, so that the parametric uncertainties on the rrih predictions be- 
come negligible, m^h is measured at the ILC with high accuracy, but a theory error from unknown 
higher-order corrections of ±0.5 GeV has been included [17]. 

The analysis is based on the comparison of the theoretical prediction [52, 358] for the ratio of 
branching ratios 

[BR{h^bb)/BR{h^ww*)]ussM .en 

[BR{h^bb)/BR{h^WW*)]gy^ ^'^ 

with its prospective experimental measurement. Even though the experimental error on the ratio 
of the two BR's is larger than that of the individual BR's, the quantity r has a stronger sensitivity 
to Ma than any single branching ratio. 

Assuming a certain precision of r, indirect bounds on Ma can be obtained. For the experimental 
accuracy of r two different values have been considered: a 4% accuracy resulting from a first phase of 
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ILC running with y/s < 500 GeV [8-10,359], and a 1.5% accuracy which can be achieved from ILC 
running at ^/s « 1 TeV [91]. Fig. l6.1l shows that a 4% accuracy on r allows to estabhsh an indirect 
upper bound on Ma for Ma values up to Ma ^ 800 GeV (corresponding to an r measurement of 
r > 1.1). With an accuracy of 1.5%, on the other hand, a precision on AMa/Ma of approximately 
20% (30%) can be achieved for Ma = 600 (800) GeV. The indirect sensitivity extends to even 
higher values of Ma- The comparison with the idealized situation where all SUSY parameters 
(except Ma) were precisely known illustrates the importance of taking into account the parametric 
errors as well as the theory errors from unknown higher-order corrections. Detailed experimental 
information on the SUSY spectrum and a precision measurement of nit are clearly indispensable 
for exploiting the experimental precision on r. 

2400 I 1 1 




Figure 6.1: The la bound on Ma, AM a, versus Ma obtained from a comparison of the precision 
measurement of r (see text) at the ILC with the MSSM prediction. The results for AMa are shown 
for a 4% accuracy of r (full line) and a 1.5% accuracy of r (dashed line). The parametric uncertain- 
ties in the prediction of r resulting from LHC /ILC measurement errors on tan /?, m^^^ ^ , m^^ , rrig , m/j , 
and nit are taken into account. Also shown is the accuracy on Ma which would be obtained if 
these uncertainties were neglected (dotted line). 



6.3.2 Distinction between the MSSM and the NMSSM 

It is possible that the MSSM and the NMSSM result in similar low-energy phenomenology. Then 
neither the LHC nor the ILC alone can identify one or the other model. Especially the Higgs 
sector does not allow the identification of the NMSSM [187] if scalar and/or pseudoscalar Higgs 
bosons with dominant singlet character escape detection due to their suppressed couplings to the 
other MSSM particles. This is illustrated in one specific example. The scenario is fixed by fi^s = 
457.5 GeV, tan /3 = 10, A = 0.5 and n = 0.2. A scan over the remaining parameters in the 
Higgs sector, the trilinear couplings Ax and A^, has been performed using NMHDECAY [181]. The 
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Figure 6.2: The masses of the two hght 
scalar Higgs bosons, mh^, m/12, and of the 
Ughtest pseudoscalar Higgs boson uia^ as 
function of the trihnear Higgs couphng in 
the NMSSM for fi^s = 457.5 GeV, tan/3 = 
10, A = 0.5 and k = 0.2, while Ax and 
Af^ have been scanned over all possible val- 
ues, hi is MSSM-like whereas /12 and oi are 
singlet-dominated Higgs particles. 
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parameter points which survive the theoretical and experimental constraints are in the region 
2740 GeV < < 5465 GeV and -553 GeV < < 0, see Fig. lO In this parameter region 
the lightest scalar Higgs boson hi is always MSSM-like whereas the second lightest scalar /12 and 
the lightest pseudoscalar oi are singlet-dominated. For —450 GeV < < — 100 GeV /12 and 
ai have very pure singlet character and masses ^ 200 GeV, which excludes their appearance in 
chargino/neutralino decays in this scenario. All other heavier Higgs bosons have masses of 0{Ax), 
hence they cannot be produced at the LHC or the ILC. In this case a combined LHC/ILC analysis 
in the neutralino sector, as carried out in Ref. [360] , is necessary. It can find the phenomenological 
differences between the two models and thus reveal also the nature of the Higgs sector. 

For A^ < — 450 GeV the smaller mass of the /i2 and a stronger mixing between the singlet 
and MSSM-like states in hi and /12 might allow a discrimination in the Higgs sector while for 

^ —100 GeV the existence of a light pseudoscalar oi may give the first hints of the NMSSM [361]. 
See Sect. 14.71 for further details. 
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Chapter 7 

Cosmological connections 



7.1 Cold dark matter and MSSM Higgs physics 

A large fraction of the total amount of matter in the universe appears to be cold dark matter (CDM). 
The density of CDM in the universe is tightly constrained by WMAP and other astrophysical and 
cosmological data [362]. Within the MSSM (or other supersymmetric extensions of the SM) it is 
well known that the lightest supersymmetric particle (LSP) is an excellent candidate for cold dark 
matter [363], with a density that falls naturally within the range [362], 

0.094 < ncBuh^ < 0.129 . (7.1) 

Assuming that the cold dark matter is composed predominantly of the lightest neutralino, a value 
of JlcDM^^ ~ 1 can easily be obtained. However, in order to exactly match eq. (|7.H) a mechanism 
in the early universe is needed that efficiently reduces the CDM density. 
One possible mechanism is the rapid Higgs pole annihilation, 

X?X?-vl^65, (7.2) 

involving the pseudo-scalar MSSM Higgs boson. This mechanism is very efficient for Ma ~ 2 m^o . 
The role of the ILC is to determine with high accuracy the masses and couplings of the involved 
particles in the process in eq. (|7.2() . The corresponding theoretical determination of JIcdm then 
has to match the experimental value of eq. (|7.H) . 

At leading order the annihilation cross section is given by 



2 Uaxx "xx 



(7.3) 



where s^j^ is the invariant mass of the two LSPs in the initial state, and VAxx coupling of the 

A to the two lightest neutralinos, given by 

VAxx = \ ^9iNii - 52iVi2) (sin/3iVi3 - cos l3Nu) , (7.4) 

with N being the unitary matrix that diagonalizes the neutralino mass matrix. 

The uncertainty of the experimental measurement of the CDM density as shown in eq. (|7.1|) 
is currently at the 10% level. How this accuracy can be matched by experimental analyses and 
theoretical evaluations will be discussed in the following subsections. 

The next generation of CDM experiments like PLANCK [364] will improve the precision of 
^CDM down to about 2%. This will require additional efforts both from the experimental as well 
as from the theoretical side. 
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7.2 Experimental issues 



The main quantities that enter the theoretical calculation of eq. (|7.2)) are the mass of the CP-odd 
Higgs boson, Ma, and its width, Ta, see eq. (|7.3() . They have to be determined with the highest 
possible accuracy. 

The ILC provides a unique tool to investigate the heavy MSSM Higgs bosons via the channels 

e+e" ^ AH^bbbb (7.5) 
e+e" ^ AH ^bb t+t' . (7.6) 

Corresponding analyses can be found in Refs. [98,365,366]. The precision that can be reached in 
Ma depends on the mass splitting between H and A. In most cases this splitting is too small to be 
disentangled at the ILC. Assuming zero splitting a mass determination better than ~ 1 GeV can 
be reached [98]. The Higgs boson width, for realistic values of Fa ^ 10 GeV could be determined 
at the 30-50% level only. 

In Ref. [366] two CMSSM points with Mh ~ Ma ~ 400 GeV have been analyzed, assuming 
1 ab~^of integrated luminosity at ^/s = 1 TeV. Here the mass average, 1/2(M^ + Mh) could be 
determined to ~ ±2.5 GeV, whereas the mass difference, 1/2(M4 — Mh) could be measured to 
about ±5 GeV. The Higgs boson width was determined better than about 25%. Assuming that 
all other relevant quantities, see Sect. 17.31 are known with negligible uncertainty, the authors used 
DarkSUSY [367] to evaluate the corresponding CDM density. The precision of I^cdm^^ was in the 
20% range, i.e., about a factor of two larger than the current experimental precision. This would 
become even worse if all other experimental uncertainties were included. 

Other experimental aspects are the neutralino properties. In order to determine Sj^-^, a mea- 
surement of the LSP mass is necessary. However, this can be achieved at the ~ 50 MeV level [8]. 
Furthermore the coupling of the lightest neutralinos to the A has to be measured. Either uaxx 
to be measured directly (which would be very difficult), or the various elements entering eq. (|7.4|) 
have to be measured. Further experimental studies in this directions are required. 

The current experimental uncertainty of ~ 10% in Ocdm^^ requires a precise determination 
of Ma and down to the level of ~ 1% for the Higgs boson mass and of about ~ 10% for the 
Higgs boson width. While the first goal seems feasible [98], the second one is very challenging. 
The situation is complicated by the strong correlation which exists between the measured widths 
of the H and A bosons and the finite mass difference between these physical states. A more 
precise determination of heavy Higgs widths would require higher integrated luminosities as well 
as improved analyses. The latter should include the study of additional signal channels, e.g., the 
Yukawa process e~^e~ — > Abb, and a more sophisticated treatment of correlation between Ta, Fh 
and Ma—Mh- It was shown [98] that by combining analyses in the HA bbbb and HA bbbT^T~ 
channels, the A boson width can be determined with a relative accuracy of about 40% for Ma ~ Mh 
= 300 GeV, assuming 500 fb^^ of integrated luminosity collected at 800 GeV centre-of-mass energy. 
Assuming a l/VA^euent dependence of the relative precision of this measurement, the accuracy can 
be improved to 20% if 2000 fb~^ of luminosity is accumulated at 800 GeV centre-of-mass energy. 
Further improvement can be achieved by inclusion of the Yukawa process e~^e~ Abb(Hbb). 
Combining the mass lineshape analysis in the e~^e~ — > Abb{Hbb) channel with the HA bbbb 
and HA — > bbbT^T~ channels would allow to improve the measurement of F^, in particular at 
high tan/3 values where the e~^e~ — > Abb process is expected to have an enhanced rate. On the 
other hand, F^ can be determined by combining mesurements of BR(A — > bb) and the rate of the 
e"'"e~ Abb process. The first measurement can be performed by exploiting the e"'"e~ HA 
process, while the second provides direct access to the partial width T^^^f^i. The total width is 
then calculated as F^ = F^_^(j5/BR(^ — > bb). Again, the precision of measurement will crucially 



88 



depend on the expected rate in the e^e~ — > Abb channel, which is governed by tan /?. Clearly, given 
the relatively poor precision on predicted by currently existing analyses, the proposed ways to 
improve the experimental accuracy on Ta should be investigated by future studies. 

In order to match the anticipated PLANCK accuracy for Ocdm^^ of about 2%, a much better 
precision in A^a and Ta will be necessary. The mass will have to be known at the ~ 0.25% level, 
whereas the width would require a ~ 3% accuracy. These goals are extremely ambitious and invite 
further investigations. 



7.3 Theoretical issues 

The theoretical calculation of ricDM^^ also requires a precise evaluation of eq. (|7.,'-{jl . The one-loop 
corrections to eq. (|7.3() are expeted at the 0(few%). Consequently, they can be omitted at the 
present level of accuracy. However, in order to match the anticipated Planck precision of about 
~ 2% a calculation of the leading one-loop corrections is necessary. Since for the 'inverse' process, 
e'^e~ XiXi! a full one-loop calculation exists [368], the tools for the required calculation appear 
to be available. 

Even if the heavy Higgs boson are not observed at the LHC and the ILC, the additional infor- 
mation from the ILC can be valuable: the calculation of r^cDM^^ can also incur large uncertainties 
due to the lack of sufficiently stringent lower bounds on the masses of A and H. We show this 
fact in Fig. 17.11 [369,370]. Here we display the variation in I^cdm^^ caused only by a variation in 
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Figure 7.1: Effect on relic density of varying the SUSY mass parameter Ma for point B' from the 
bulk region of mSUGRA (left plot) and the point LCC2 from the focus point region of mSUGRA 
(right plot). The horizontal (green-shaded) region denotes the 2a WMAP limits on the dark matter 
relic density. The red line shows the variation of the relic density as a function of Ma- The blue 
dot in each plot denotes the nominal value for Ma at point B' (left plot) and the point LCC2 (right 
plot). Taken from Refs. [369,370]. 

Ma- In the left plot, we have chosen as our nominal point one from the 'bulk' region of mSUGRA 
parameter space. The blue circle identifies the nominal values of f^cDM''-^ and Ma for the point. 
The green cross-hatched region denotes values consistent with the WMAP determination of the 
relic density. The red line illustrates how the value of Ocdm^^ changes as the parameter Ma is 
changed. Since the slope of the red line is small as it passes through the nominal value, i^cDM^^ 
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appears rather insensitive to the exact value of M4. However, it is possible that the value of Ma 
cannot be determined at a collider (if = 500 GeV for example), and only a lower bound can 
be attained. If this is the case, when calculating Ocdm^^j one must allow Ma to vary between 
its lower bound and infinity. From the plot, it is clear that this introduces large uncertainties in 
^^CDm/i^ if the lower bound is not much above 200 GeV. It is likely that this will be the best 
lower limit available from the LHC at this point [369,370]. However, a 500 GeV linear collider 
will raise this lower limit to approximately 250 GeV (and higher energies can give rise to larger 
increases). From the plot, it can be seen that the extra 50 GeV gained from the ILC will vastly 
improve the uncertainty in the calculation of ricDM^^ due to experimental uncertainty in Ma- The 
right plot shows the same story for the case of a point [371] from the 'focus point' region [372]. If 
the value of Ma can be determined, Ocdm^^ will be very insensitive to its exact numerical value. 
This is unlikely to happen since Ma ~ 3 TeV. However, VLcT>Mh'^ will be very sensitive to the lower 
experimental bound on Ma if it is around 200 GeV. 



7.4 The NMSSM case 

7.4.1 Cold dark matter and NMSSM Higgs physics 

Since the NMSSM has five neutralinos and two CP-odd Higgs bosons, there are many new ways in 
which the relic density of the x? could match the observed dark matter density; see, e.g., [201,202]. 
The latter group has made their code publicly available. In the NMSSM the uia^ and rriaj masses are 
quite unconstrained by LEP data and theoretical model structure, implying that XiXi ~^ ^1,2 X 
could be the primary annihilation mechanism for large swaths of parameter space. It is relatively 
easy to construct NMSSM models yielding the correct relic density even for a very light neutralino, 
100 MeV < m^o < 20 GeV. Even after including constraints from Upsilon decays, b — > 57, Bs — > 
ix^ ^~ and the magnetic moment of the muon, a light bino or singlino neutralino is allowed that 
can generate the appropriate relic density. 

The LSP is a mixture of the bino, neutral wino, neutral higgsinos and singlino, which is the 
superpartner of the singlet Higgs. The lightest neutralino therefore has a singlino component 
in addition to the four MSSM components; its eigenvector x? can be written in terms of gauge 
eigenstates in the form 

X? = e^&l + e^^o + ewW"^ + esB + e,5, (7.7) 

where e^j are the up-type and down-type higgsino components, e^y, es are the wino and bino 
components and is the singlet component of the lightest neutralino. Similarly, the CP-even and 
CP-odd Higgs states are mixtures of MSSM-like Higgses and singlets. For the lightest CP-even 
Higgs state we can define (see also Sect. 14.7.1(1 



, 1 



e„Re {Hi - Vu) + e^Re (i?o - v^) + 6Re {S - s)] . (7.8) 



Here, Re denotes the real component of the respective state. Lastly, the lightest CP-odd Higgs 
can be written as (a similar formula also applies for the heavier 02) 

ai = cos Oaj^ Amssm + sin 9a^ As, (7.9) 

where As is the CP-odd piece of the singlet and ^mssm is the state that would be the MSSM 
pseudoscalar Higgs if the singlet were not present. Here, 9ai is the mixing angle between these two 
states. 
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In the NMSSM context, when annihilation proceeds via one of the CP-odd Higgs bosons the 
calculation of the relic Xi density is much more flexible than in the MSSM; see Sect. 17.11 The 
general formula for the thermally averaged cross section takes the form 
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2{Xeued - KeJ 



(7.10) 



where c/ is a color factor, equal to 3 for quarks and 1 otherwise. For this result, we have assumed 
that the final state fermions are down-type. If they are instead up-type fermions, the tan^ P factor 
should be replaced by cot^ f3. For medium to large values of tan f3, neutralino dark matter can 
avoid being overproduced for any oi mass below ~ 20 - 
values of tan /3, a lower limit on can apply as well. 
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Figure 7.2: We display contours in niai - m^o parameter space for which il/i^ = 0.1. Points above 
or below each pair of curves produce more dark matter than is observed; inside each set of curves 
less dark matter is produced than is observed. These results are for a bino-like neutralino with a 
small higgsino admixture (e^ = 0.94, = 0.06). Three values of tan (3 (50, 15 and 3) have been 
used, shown as solid black, dashed red, and dot-dashed blue lines, respectively. The dotted line is 
the contour corresponding to 2m^o = niaj. For each set of lines, we have set cos^ 6a^ = 0.6. 

Fief. l7.2l shows typical NMSSM results. There, the ai composition is chosen such that cos^ 6ai = 
0.6 and the neutralino composition is specified by = 0.94 and = 0.06. These specific values 
are representative of those that can be achieved for various NMSSM parameter choices satisfying 
all constraints. For each pair of contours (solid black, dashed red, and dot-dashed blue), O/i^ = 0.1 
along the contours and the region between the lines is the space in which the neutralino's relic 
density obeys fi/i^ < 0.1. The solid black, dashed red, and dot-dashed blue lines correspond to 
tan /3=50, 15 and 3, respectively. Also shown as a dotted line is the contour corresponding to the 
resonance condition, 2m^o = nia^ . For the tan /3=50 or 15 cases, neutralino dark matter can avoid 



91 



being overproduced for any ai mass below ~ 20 — 60 GeV, as long as m^o > rub- For smaller values 
of tan (3, a lower limit on niai can apply as well. 

For neutralinos lighter than the mass of the 6-quark, annihilation is generally less efficient. This 
region is shown in detail in the right frame of Fig. 17.21 In this funnel region, annihilations to cc, 
T^T^ and ss all contribute significantly. Despite the much smaller mass of the strange quark, its 
couplings are enhanced by a factor proportional to tan/3 (as with bottom quarks) and thus can 
play an important role in this mass range. In this mass range, constraints from Upsilon and J/^l^ 
decays can be very important, often requiring fairly small values of cos6'a^. 

The above discussion focused on the case of a mainly bino LSP. If the LSP is mostly singlino, it 
is also possible to generate the observed relic abundance in the NMSSM. A number of features differ 
for the singlino-like case in contrast to a bino-like LSP, however. Most importantly, an LSP mass 
that is chosen to be precisely at the Higgs resonance, nia, — 2m^o, is not possible for this case: rua. 
is always less than 2m^o by a significant amount. Second, in models with a singlino-like LSP, the 
ai is generally also singlet-like and the product of tan^ /3 and cos^ 6a^ , to which annihilation rates 
are proportional, see Eq. 1)7. 10(1 . is typically very small. This limits the ability of a singlino-like 
LSP to generate the observed relic abundance. The result is that annihilation is too inefficient for 
an LSP that is more than 80% singlino. However, there is no problem having m^o ~ 'ma^/2 so as 
to achieve the correct relic density when the Xi is mainly bino while the ai is mainly singlet. 

Finally, for the scenarios with a low F value (see eq. (|4.12|) ). one has m^i < 2mh and the ai is 
very singlet-like, with cos^ Oa^ ^ 0.015. In this case, adequate annihilation of a very light Xi via 
XiXi ai ^ X occurs only if m^o ~ ma^/2. This requires a rather fine adjustment of the Mi 



bino soft mass relative to nia, ■ Because the oi is so light in the low fine-tuning scenarios, if m^o is 
significantly above 2mf, then consistency with relic abundance limits requires that XiXi annihilation 
proceed via one of the more conventional co-annihilation channels or via XiXi 0,2 ^ X. The 
latter case is only applicable if m^o ^ 200 GeV (with typically ^ 400 GeV). 

7.4.2 Direct dark matter detection 

In Ref. [201] it is estimated that the neutralino-proton elastic scattering cross section is on the 
order of 4 x 10~^^ cm^ (4 x 10~^ fb) for either a bino-like or a singlino-like LSP. This value may be 
of interest to direct detection searches such as CDMS, DAMA, Edelweiss, ZEPLIN and CRESST. 
To account for the DAMA data, the cross section would have to be enhanced by a local over-density 
of dark matter. 

Prospects for direct detection of dark matter in the low fine-tuning scenarios are rather con- 
strained. Since the ai is so singlet in nature, the only exchange of importance is hi exchange. In 
the low-F scenarios, hi is almost entirely Hu. In particular, the Hii component of the hi is ~ 0.1, 
and correspondingly ~ 0.99. Whether or not m^o is below nifj {m^o < rrib requiring annihilation 

via XiXi — > ai — > X), the typical composition of a bino-like Xi is such that eb > 0.8 and and Ed 
take a range of values from 0.1 up to 0.5. Keeping only the s-quark contribution and the dominant 
EdS,u piece in the external factor, we obtain 



f^elastic ~ 5 X 10-64 ( ) ( - ) ( = ) fb. (7.11 



.2 ( ^dV /lOO GeVy /tan/3^ ^ 
,025^ I rrih, ) \ l0 



Although m^o can be chosen in the range ^ 20 GeV (assuming appropriate coannihilation to 
achieve proper relic density), which is relatively optimal for experiments like ZEPLIN and CRESST, 
the predicted cross section is much below their expected sensitivity. If annihilation is via XiXi ~^ 
ai — > X, the required mass m^o ~ < mf, is not very optimal and there would be no hope of 

direct detection at such experiments. 
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7.4.3 The role of the ILC 



In the case where m^o < 20 30 GeV, and especially if m^o < mb, the best that a hadron collider 

can do will probably be to set an upper limit on m-o. Determining its composition is almost certain 

A 1 

to be very difficult. Note that the m^o — m^o mass difference should be large, implying adequate 

room for ^2 ~^ ^Xi ^ search for lepton kinematic edges and the like. (Of course, X2 ~^ ^iXi will 
also probably be an allowed channel, with associated implications for hi detection in SUSY cascade 
decays.) If the dominant annihilation channel is XiXi ai ^ X, determination of the properties 
of the ai would also be critical. A light singlet-like ai is very hard to detect. At best, it might 
be possible to bound cos^aj by experimentally establishing an upper bound on the WW aiai 
rate (proportional to cos^^aj- Thus, the ILC would be absolutely essential. Precise measurement 
of the Xi mass and composition using the standard ILC techniques should be straightforward. 
A bigger question is how best to learn about the ai at the precision level. Of course, we will 
have lots of ai's to study from hi aioi decays. The events will give precise measurements of 
glzhi'^^ihi aiai)BR(ai X)BR(ai Y), where X,Y = bb,T+T-. The problem will be 
to unfold the individual branching ratios so as to learn about the oi itself. Particularly crucial 
would be some sort of determination of cos which enters so critically into the annihilation 
rate. (One can assume that a tan /3 measurement could come from other supersymmetry particle 
production measurements and so take it as given.) There is some chance that WW — > aiai and 
Z* — ZaiQi, with rates proportional to cos^ , could be detected. The cos^ = 1 rates are very 
large, implying that observation might be possible despite the fact that the low-i*" scenarios have 
cos^ < 0.01. One could also consider whether 77 — > ai production would have an observable 
signal despite the suppression due to the singlet nature of the ai. Hopefully, enough precision 
could be achieved for the oi measurements that they could be combined with the Xi precision 
measurements so as to allow a precision calculation of the expected XiXi ^ ai ^ X annihilation 
rate. A study of the errors in the dark matter density computation using the above measurements 
as compared to the expected experimental error for the Qh^ measurement is needed. 

7.5 Electroweak Baryogenesis and quantum corrections to the 
triple Higgs boson coupling 

Baryogenesis is one of the fundamental cosmological problems. Electroweak baryogenesis [30] is 
a scenario of baryogenesis at the electroweak phase transition, so that the Higgs sector plays an 
essentially important role. In the symmetric phase (T > Tc), where Tc is the critical temperature, 
baryon number is generated by sphaleron processes. For successful electroweak baryogenesis, it is 
required that the sphaleron decouples just after the phase transition. In the SM this condition 
implies ipc/Tc ^ 1, where ipc is the vacuum expectation value at Tc. i.e., the phase transition is 
strongly first order. This is qualitatively described in the high temperature expansion, where the 
effective potential is expressed by 

Feff = DiT^ - r2)v?2 _ ^ ^ . . . _ ^7_;^2) 

One can see that (pc = 2ETc/Xt- Hence, the sphaleron decoupling condition yields E/\t ^ 0.5. 

In the SM, where E = ^i;-(2m^ + m^), this gives the upper bound for the mass of the Higgs 
boson nifi < 46 GeV. A more quantitative analysis shows that the phase transition turns into a 
smooth cross-over for m/j ^ 80 GeV. Therefore, the electroweak baryogenesis cannot be realized in 
the framework of the SM with the present experimental limit mh ^114 GeV [4]. 
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7.5.1 Strong first order phase transition by the additional bosonic loop contri- 
bution 



The condition of the sphaleron decouphng is compatible with the current experimental results when 
we consider extensions of the SM, which give larger values of the coefficient E of the cubic term in 
Eq. H7.12() than in the SM prediction. Some additional bosons can contribute to enhance E at one- 
loop level when they show the non-decoupling property: i.e., when they receive masses mainly from 
the vacuum expectation value v. In the MSSM, light stops may be a candidate of such additional 
bosons [373]. An alternative viable model may be the Two-Higgs-Doublet Model (THDM) [374]. 

The phenomenological influence of such succesful baryogenesis scenarios on the Higgs physics 
has been discussed in the THDM and MSSM [75]. It was shown that the condition on the finite 
temperature effective potential for a successful scenario exactly corresponds to the condition of a 
large non-decoupling quantum effect due to extra Higgs boson loops on the hhh coupling of the 
lightest Higgs boson h. In the THDM the leading contributions to the one-loop corrections to the 
hhh coupling is given in the SM-like limit (sin(a — /?) ~ — 1) by [74] 



A^L(THDM) 




(7.13) 



where M is the soft-breaking mass parameter of the discrete Z2 symmetry. 

At finite temperatures, in the case with <^ m\ ($ = H, A, H^), the coefficient E of the 
cubic term in Eq. 1)7. 12() is evaluated as ii^ = 12^^ i^^w + 3m^ + + m\ + 2m^±), and the 
condition of sphaleron decoupling can be satisfied for rrih ^115 GeV for large values of heavy Higgs 
boson masses. In this case, from Eq. (|7.13p the effects of the heavy Higgs boson loops on the hhh 
coupling are enhanced by m|,. These loop effects do not decouple in the large mass limit, and yield 
large deviations in the hhh coupling from the SM prediction. Such large non-decoupling effects 
can modify the hhh coupling by more than 10% [74] (maximally ~ a few times 100% within the 
requirement from perturbative unitarity [375]), so that the deviation is expected to be detected at 
ILC [77,78]. In the same scenario, the effects on the quartic Higgs boson coupling are also discussed 
in Ref. [376]. 

This connection between the condition for successful electroweak baryogenesis and the large 
non-decoupling effect on the hhh coupling is a common feature, and can be seen in some other 
models such as the MSSM with light stops. In the MSSM, the effect of additional stop loop 
contributions to the hhh coupling is expressed by using the high temperature expansion as 

AA...(MSSM) ^^ 

Am/i(SM) mfml 

where Ai?^^ is the contribution of the light stop loop to the cubic term in the finite temperature 
effective potnetial in Eq. (|7.12|) . 

Therefore, these scenarios of electroweak baryogenesis can be tested by precisely measuring the 
hhh coupling at an ILC (see Figure 173)) . 

7.5.2 Strong first order phase transition in the SM with a low cutoff scenario 

Recently the possibility of electroweak baryogenesis has been studied in the SM with a low cut-off 
scale, by adding the higher dimensional Higgs boson operators [377,378]. In contrast to previous 
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Figure 7.3: The contour of the radiative correction of the triple Higgs boson couphng constant 
overlaid with the line ^Pc/Tc = 1 in the m^-M plane for mft=100, 120, 140 and 160 GeV. Above 
the critical line, the phase transition is strong enough for a successful electroweak baryogenesis 
scenario [75]. 



approaches which rely on large cubic Higgs interactions, they allow the possibility of a negative 
quartic coupling while the stability of the potential is restored by higher dimensional operators [378]: 

= A (^c&tcD - , (7.15) 

where ^ is the SM electroweak Higgs doublet. At zero temperature, the CP-even scalar state can 
be written in terms of the zero temperature vacuum expectation value ((/?) = — 246 GeV and 
the physical Higgs boson H: ^ = ip/V2 = {H + vo)/^/2. One can minimize the zero temperature 
potential to find A = m'j^/{2vQ) and v = vq in terms of physical parameters m^^ and vq. 

The effective potential at finite temperatures is approximately obtained by adding the thermal 
mass to the zero temperature potential V{if,T) = + V^((/j, 0), where c is generated by 

the quadratic terms, T'^mf, in the high temperature expansion of the thermal potential. The 
quartic term, which is neglected here, is also induced but its contribution is a few percent in the 
physically interesting region where a strongly first-order electroweak phase transition occurs. The 
cubic interaction is also induced at finite temperatures but it has a smaller role in this scenario, thus 
neglected here. Prom the requirements of the global minimum of the vacuum at zero temprature, a 
first order phase transition, and positivity of the thermal mass term, the cutoff scale is constrained 
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as 



max 



H 



m 



< A < \/3- 



(7.16) 



H 



where rric = vq\J {Ay^ +3g'^ + g'^)/8 ~ 200 GeV, where g and g' are SU{2)l and U{1)y gauge 
couplings, and yt is the top Yukawa couphng. The left plot of Fig. 17.41 shows the contours of 
constant Vc/T^ in the A vs. rrijj plane, where is the critical temperature and Vc is the vacuum 
expectation value at Tc- The results are encouraging and motivate a full one-loop computation of 
the thermal potential embedded in a complete baryogenesis scenario. 

The value of the cutoff scale A affects the prediction of the HHH and HHHH couplings. 
Describing the zero temperature Higgs self-interations by V{H) = mj^ 
one has 



H'^/2+fiH^/3l+r]H^/Al + - 



2 S 

3—^ + 6-^. 

A2 ■ 



2 2 

11 = 3^ + 36-^ . 
vl A2 



Vq 1^ c/Q 

The right plot of Fig. 17.41 shows that the contours of /i//isM — 1 in the A vs. rrijj plane. 
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Figure 7.4: [Left] Contours of constant v^/Tc from 1 to oo. The shaded (blue) region satisfies the 
bounds of Eq. (|7.16|) . [Right] Contours of constant /i//isM — 1 in the A vs. plane. The dashed 
lines delimit the allowed region defined in Eq. (|7.16|) . For details, see [378]. 
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Chapter 8 



Impact of Machine and Detector 
Performance on Precision of 
Measurements in the Higgs Sector 



The physics processes at the ILC have been studied for about 15 years and resulted in a variety 
of studies for discovery and precision measurements. Many requirements on the detector and to a 
lesser extend on the machine development have been established. The impact of the detector (res- 
olution, material budget, and coverage) and machine aspects which directly influence the physics 
performance (beam energy and luminosity spread, background induced in the detector) are impor- 
tant aspects to be studied further. In this chapter we collect the existing results and discuss how 
the detector could be optimized in view of Higgs boson physics. Finally we present a strategy for 
future studies to facilitate the final goal of Higgs precision measurements. 

We start this section be recalling the main requirements for the ILC detector [379] . 

• Momentum resolution: ai/p '-^ 5 x 10~^ GeV~^ (a factor of 10 better than that achieved at 
LEP). Good momentum resolution is important for the reconstruction of the leptonic decays 
of Z bosons. This is particularly relevant for e~^e~ HZ, which is best suited to study the 
Higgs boson properties in detail. This works especially well by investigating the Higgs boson 
recoiling against the Z, with the clean decay Z 1^"^ l^^" ■ 

• Impact parameter resolution: Efficient b and c quark tagging is important, implying a good 
impact parameter {do) resolution: 



where p is particle momentum and 9 is the polar angle. This is a factor of three better than 
the resolution obtained at SLD. 



identify gauge bosons from their hadronic decays to jets. 

• Hermiticity: Hermiticity down to 5 mrad is needed for searches for missing energy signals 
from new physics. 




(8.1) 



• Jet energy resolution: ge/E ^ 0.3/-\/^[GcVj in order to be able to directly reconstruct and 
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• High granularity: Events at the ILC will have high local particle densities due to the boosted 
final states and the fact that many of the interesting physics processes result in final states 
with six or more jets, especially in the case of Higgs events. The ILC detector is required to 
have high granularity calorimeters to enable efficient separation of close-by showers and good 
pattern recognition in the tracking system. 

8.1 Existing Analyses 

A number of studies investigating the relation between machine and detector performance and the 
precision of measurements in the Higgs sector have been presented at this Snowmass workshop and 
documented as LC Notes or e-prints. 

8.1.1 Beam related issues 

In Ref. [56], discussing the prospects of the Higgs boson mass measurements at ILC, the impact 
of beam related systematic effects on the Higgs mass measurements have been evaluated. The 
following experimental effects have been investigated: 

• uncertainty in the beam energy measurement; 

• beam energy spread; and 

• uncertainty in the differential luminosity spectrum measurements. 

Since the analysis relies upon kinematic fits exploiting energy-momentum conservation, a shift 
in center-of-mass energy will inevitably introduce a bias in the measured Higgs mass. This is 
illustrated in Fig. 18.11 In order to keep systematic uncertainty due to a bias in the beam energy 
measurement below the statistical error on the Higgs mass, the beam energy has to be controlled 
with relative precision better than 10~^. 

To estimate the impact of the beam energy spread, Gaussian smearing of the beam energy has 
been applied at the stage of signal event generation. It has been shown that for the case of 1% 
energy spread for both electron and positron beams, the statistical precision in the Higgs boson 
mass measurement degrades from 45 to 50 MeV in the HZ bbqq channel and from 85 to 90 MeV 
in the HZ qqi~^i~ channel. For a TESLA-like machine the expected energy spread amounts to 
0.15% for the electron beam and 0.032% for the positron beam [8]. For these energy spreads no 
significant degradation of the precision in the Higgs boson mass measurement is observed. 

The energy spectra of the colliding electrons and positrons at the linear collider will be sig- 
nificantly affected by photon radiation of an electron/positron in one bunch against the coherent 
field of the opposite bunch. This effect is referred to as beamstrahlung. To have a fast simulation 
of beamstrahlung, the program CIRCE [60] has been written. It assumes that the beamstrahlung 
in the two beams is equal and uncorrelated between the beams. The spectrum is parameterized 
according to 



where x is the ratio of the colliding electron/positron energy to the initial energy of the undisrupted 
beam. The parameters depend on operational conditions of the linear collider. The normalization 
condition, J f(x)dx = 1, fixes one of these parameters, leaving only three of them independent. The 
default parameters for the ILC (taken from TESLA) operated at the centre-of-mass energy of 350 
GeV are: 





ao = 0.55, ai = 0.59, aa = 20.3, as = -0.63. 
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Figure 8.1: The spectrum of the fitted val- 
ues of the Higgs boson mass as obtained from 
200 independent signal samples in the HZ 
qqi~^i~ channel for the case when both elec- 
tron and positron beam energies are overesti- 
mated by 25 MeV (dotted histogram), when 
they are underestimated by 25 MeV (dashed 
histogram) and when no shifts are introduced 
to the beam energies (solid histogram). 



It was shown that from the analysis of the acolinearity spectrum of Bhabha scattering events, the 
parameters a, can be determined with a precision of about 1% [380]. With such a precision on 
the parameters aj, the uncertainty in the differential luminosity spectrum is expected to have a 
negligible effect on the Higgs mass measurement. 

Beamstrahlung photons can interact and among others produce hadronic final states, 77 
hadrons. At a cold machine, the average number of these interactions with the center-of-mass 
energy of colliding photons in excess of 5 GeV is expected to be about 0.25 (0.4) per bunch crossing 
in e+e~ collisions at ^/s = 500 (800) GeV. This effect gets even stronger at a 77 collider. About 
two 77 hadrons events are expected on average per bunch crossing at the ILC photon collider 
(taken from TESLA) for \/igg = 400 GeV at nominal luminosity^. These parasitic interactions 
will inevitably contaminate physics events and deteriorate the performance of physics analyses. 
The impact of the hadronic 77 background on selected Higgs analyses has been investigated in 
Reference [381]. Particles resulting from 77 — > hadrons interactions tend to be produced at low 
polar angles. Therefore, for processes like e^e~ — > HZ, where mainly central production of bosons 
is expected, the contamination from parasitic interactions can be reduced by rejecting particles 
produced at low polar angles. This allows for a partial recovery of the analysis performance. With 
special rejection criteria, allowing for suppression of 77 — > hadrons background, the precision on 
the Higgs boson mass measurement in the HZ bbqq channel at ^/s = 500 GeV with 500 fb~^ 
degrades only by a few MeV. However, for the W^PF-fusion process the effect is expected to be more 
severe. For the cold machine, the presence of the 77 — > hadrons background results in degradation 
of the signal-to-background ratio in the e~^e~ Wvv — > bbvv channel by about 15% at ^/s = 1 
TeV. 

'^This effect was taken into consideration in the dedicated analyses related to Higgs boson studies at a photon 
coUider [63, 100] 
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8.1.2 Jet identification and resolution 



Many signal signatures at ILC will involve multi-jet final states. For their analysis, accurate jet 
identification and measurement is crucial. Efficient separation of ZZ and WW decay modes of a 
heavy Higgs boson (Fig. 13. is possible with a jet energy resolution Ai^jetZ-Ejet < 0.3/-\/-E'jet- A 
recent study has investigated the impact of jet energy resolution on the measurement of the Higgs 
boson branching ratio to W bosons in the Higgs-strahlung process [382]. 
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Figure 8.2: Distribution of the reconstructed Higgs boson mass in the sample of selected Higgs- 
strahlung events with subsequent decay of the Higgs into W boson pairs and hadronic decays of the 
Z boson. Only semileptonic decays of W pairs are considered. Left plot corresponds to jet energy 
resolution of 0.3/-\/£jet, right plot corresponds to jet energy resolution of O.G/i/^^jet- Selected 
sample corresponds to an integrated luminosity of 500 fb~^ collected at \fs = 360 GeV. 

The analysis is based on selection of e~^e~ — > HZ events with subsequent semileptonic decays of 
W pairs, WW* jj^v, and hadronic decays of the Z boson, Z qq. The study is done for a light 
Higgs boson, mn = 120 GeV, assuming 500 fb^i of data collected at ^ = 360 GeV. Three scenarios 
are considered, corresponding to assumed jet energy resolution of 0.3/i/£'jct (benchmark goal for 
the ILC detector), 0.6/-^iijot (jet energy resolution for a typical LEP detector) and OA/^^E^ (some 
intermediate value between the first two scenarios). Fig. l8.2l Dresents the reconstructed Higgs boson 
mass for the first and second scenarios. A considerable degradation in the signal-to-background 
ratio from 5 to 2.5 is observed for jet energy resolution of 0.6/ ^/E^ compared to 0.3/ ^/E^. This 
results in deterioration of the precision in the H WW* branching fraction by 20% which is 
equivalent to 45% loss in the integrated luminosity. 

A similar analysis [383] has been performed which investigated the impact of jet energy res- 
olution on the Higgs boson mass measurement in the four-jet channel. It has been shown that 
a deterioration in the jet energy resolution from 30%/\/^ to 60%) /^/E results in a decrease in 
the Higgs mass accuracy from 42 MeV to 50 MeV. Like in the previous analysis, this roughly 
corresponds to an effective luminosity loss of 45%. 

In Sect. the method to determine the CP quantum number of the Higgs boson by studying 
transverse spin correlation effects in the H t^t^ decay has been described. Efficient iden- 
tification of the tau decay mode is crucial for this measurement and sets stringent requirements 
on the performance of the electromagnetic calorimeter in combination with a dedicated pattern 
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recognition algorithm [384] . Fig. 18.31 illustrates the separation of r — > vrz/ and t ^ pv modes in the 
distribution of the reconstructed tau-jet mass. The result has been obtained with dedicated recon- 
struction tools aiming to identify and separate showers resulting from the charged pion and from 
the TT^ 77 decay. The simulation is performed for a TESLA-like detector with a high granularity 
tungsten-silicon ECAL. The quality of separation between — > -k^v and — > p^v — > tt^t^^v 
modes is quantified in Tab. 18.11 These results set one of the benchmark conditions for the ECAL 
optimization studies. 



Decay Mode 


Jet mass < 0.2 GeV 


0.2 GeV < Jet mass < 2 GeV 


p^U TT^TT^U 


82% 
2% 


17% 
19% 



Table 8.1: Fraction of the r — > tt v and r — > p — > tt tt u final states with the tau-jet mass 
below 0.2 GeV and in the window 0.2 — 2 GeV. Distributions correspond to the samples in which 
the distance between Tr"*" entering calorimeter and the nearest photon is greater that 1.5 cm, one 
and a half calorimeter cell size. 



8.1.3 Identification of b- and c-quarks 

For the study of the hadronic branching fractions, an excellent h- and c-tagging capability of 
the detector is crucial. Tools have been developed which implement in detail the procedure of 
tagging b- and c-hadrons [99]. Fig. 18.41 illustrates performance of h- and c-tagging in terms of h 
(c) purity as a function of b (c) efficiency as obtained with dedicated tools for the micro-vertex 
detector configuration foreseen for a TESLA-like detector. Recently, a dedicated analysis based on 
these detailed flavor-tagging tools has been performed [70] to evaluate the ILC potential for the 
measurement of hadronic decays of a light Higgs boson with mn = 120 GeV. The study is based on 
the the selection of an inclusive sample of hadronic Higgs decays using the Higgs-strahlung process 
and exploiting all possible decay modes of the Z boson. Each selected event is assigned quantities 
that quantify the probability to contain h- or c-jets. These quantities are referred to as b- and 
c-tag variables. The branching ratios H ^ bb, H ^ cc and H ^ gg are determined from the 
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Figure 8.4: The performance of b- and c- 
tagging as envisaged for a vertex device of the 
TESLA detector. Performance is quantified in 
terms of 6/c purity as a function of b/c effi- 
ciency. 



Figure 8.5: The 2D distribution of event 6-tag 
vs. c-tag in the sample of the selected Higgs- 
strahlung events with H hadrons and Z 
qq [70]. 



fit of 2D distribution of 6-tag vs. c-tag variables (Fig. 18.5(1 with three free parameters, which are 
normalization factors quantifying the fractions oi H ^ bb, H ^ cc and H ^ gg events in the final 
selected samples. Results of this study are presented in Tab. 18.21 



Decay 


Z l+t 




Z ^ uu 


Z ^ qq 


Combined 




Relative precision (%) 


H ^bb 


3.0 




2.1 


1.5 




1.1 


H ^cc 


33.0 




20.5 


17.5 




12.1 


H^gg 


18.5 




12.3 


UA 




8.3 



Table 8.2: Relative precision on 
the measurement of the hadronic 
branching ratios of Higgs for mn 
= 120 GeV. The analysis is per- 
formed for = 350 GeV, as- 
suming an integrated luminosity of 



500 fb-^ 

Another study (reported at Snowmass) has evaluated the impact of the vertex detector parame- 
ters on the measurement of the hadronic branching fractions of the Higgs boson [385]. The analysis 
is performed for the SM Higgs boson with a mass of 127 GeV and for one representative MSSM 
parameter point giving the same mass for the light supersymmetric Higgs. The detector response is 
simulated with the SGV program [386]. Two configurations of the vertex detector are considered. 
The first configuration has five layers, located at the distance 15, 26, 37, 48, 60 mm away from the 
beam axis. In the second configuration the innermost layer is removed. For both configurations, the 
precision on T(H — > cc) and T{H — > bb) has been investigated as a function of the layer thickness 
and spatial resolution of a single layer. The analysis is based on the selection of Higgs-strahlung 
events. The prospective measurement is simulated assuming an integrated luminosity of 500 fb~^ 
collected at 350 GeV center-of-mass energy. The study revealed strong dependence of the relative 
precision of the T{H — > cc) measurement on the vertex detector characteristics as illustrated in 
Fig. 18.61 The dependence is found to be less pronounced for the H ^ bb channel. 
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Figure 8.6: The dependence of relative precision on T{H function of spatial resolution 

of single layer (left plot) and layer thickness (right plot). In the left plot layer thickness is fixed 
to 50 iim. In the right plot spatial resolution is fixed to 2/xm. Filled dots correspond to the 
vertex configuration with five layers. Open dots show results for the case when the innermost layer 
is removed. Results for the SM Higgs are indicated with circles, for MSSM Higgs with square 
markers. 

8.2 Higgs Physics and Detector Optimization 

The International Linear Collider project enters now the phase of detector design and optimiza- 
tion. At this stage, the detailed analyses studying the impact of the detector design and event 
reconstruction procedure on the precision of the measurements at ILC are of crucial importance. 
The studies discussed in previous sections have been performed using fast parametric Monte Carlo 
programs, in which the detector response and event reconstruction performance are simulated in 
an idealized manner for a specific detector configuration option. Such an approach is illustrated 
by the Simdet [387] program, which found wide acceptance within the ECFA studies. Simdet 
parameterizes the response of the ILC detector with fixed detector geometry and therefore is not 
appropriate for studies relating various detector configurations to the attainable precision of mea- 
surements in various physics channels. Simdet performs a fast Monte Carlo by smearing the 
particle momenta and energies according to subdetector resolutions obtained from studies with 
full simulation programs based on the GeantS or Geant4 packages. Parameterization of the de- 
tector response in Simdet relies on Monte Carlo studies performed with the full GeantS-based 
simulation program Brahms on single particle samples. Results of these studies are then extrap- 
olated to events with complex topology, e.g., multi-jet final states, assuming highly efficient pat- 
tern recognition in the tracking system and calorimeters in a high local particle density environ- 
ment. Clearly, the physics potential of the ILC will depend on both detector performance and 
performance of the corresponding reconstruction software, and the issue is still open whether the 
performance of the detector and event reconstruction are realistically implemented in the fast para- 
metric Monte Carlo programs. In order to answer this question and to perform detailed and reliable 
detector optimization studies, future analyses related to Higgs physics at the ILC should be done 
with a full detector simulation based on GeantS or Geant4 implementations of existing detector 
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designs and with consistent and realistic reconstruction software. 

8.2.1 Tools for Full Detector Simulation 

Currently, activities are ongoing worldwide which aim at the development of realistic simulation 
tools based on the Geant4 package as well as reconstruction software which can be used for the 
detector optimization studies. A big progress in the development of such tools has been reported 
at Snowmass. Most of these tools use the common output format realized in the LCIO data 
model [388]. This feature facilitates direct comparison between the various detector models im- 
plemented in the different simulators. The list of available full simulation programs for the linear 
collider detector is given below: 

• Brahms [389]: GeantS based package which implements the TESLA TDR geometry in a 
hard-coded manner. The program reads HEPEVT files and produces an output in the ASCII 
or LCIO format. Since the detector geometry is fixed, the Brahms package cannot be utilized 
for detector optimization studies. However, it may still be useful for revision of currently 
existing Higgs analyses with the full realistic simulation. 

• SLIC [390]: Geant4 based simulation package that uses an XML geometry input format 
called LCDD. This package includes various detector models developed within the compact 
detector concept with the silicon tracker (SiD) as well as large detector (LDC) and huge 
detector (GLD) concepts with the TPC as the main tracking device. The program reads a 
StdHep file for the event input and produces an output in the form of a LCIO file. 

• Mokka [391]: Geant4 simulator that retrieves geometry information from a MySQL database. 
It reads binary StdHep or ASCII HEPEVT files and writes a LCIO file. Although originally 
Mokka was designed for simulation of the large detector with the TPC as the main tracking 
device, the latest release of Mokka contains also models with the silicon tracker. The package 
provides an option for coherent scaling of various detector components such as the TPC, 
ECAL and HCAL, which allows to optimize the length and radius of the tracker, the depth 
and sampling frequency of calorimeters, and other geometrical parameters of the detector. 

• Jupiter [392]: Geant4 based program for simulation of the GLD detector. This program 
uses as input StdHep or HEPEVT files. The geometry of the detector is specified using 
XML interfaces. The data flow controller is based on ROOT. However, work is ongoing to 
implement output in the LCIO format to enable comparison with other detector concepts. 

8.2.2 Tools for Event Reconstruction 

The most promising strategy for event reconstruction at the ILC is based on the particle flow 
concept, implying reconstruction of the four-vectors of all particles produced in an event. The 
particle flow algorithm works best at moderate energies of individual particles, below about 100 
GeV. In this regime, the tracking system reconstructs the momentum of the charged particles with 
an accuracy better than that achievable for the energy measurements with calorimeters. This is 
illustrated in Figure IH771 Hence, in order to attain a better reconstruction of events, the charged 
particle measurement must be solely based on the tracking information. The crucial step of the 
particle flow algorithm is the correct assignment of calorimeter hits to the charged particles, and 
eflicient separation of close- by showers produced by charged and neutral particles. A clustering 
algorithm is then performed on calorimetric hits not assigned to the charged particles to reconstruct 
photons and neutral hadrons. The reconstruction of photons and neutral hadrons buried in jets is 
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Figure 8.7: The track momentum resolution 
for charged objects and the energy resolu- 
tion for photons and hadrons as a function 
of particle energy. Plotted curves correspond 
to the resolutions envisaged for the TESLA 
detector. 



nearly impossible with a typical LEP detector. This explains why LEP detectors could reconstruct 
jet energy with something like 60%/\/^ (where E represents the jet energy). The goals of the 
physics program at the linear collider requires however much better jet energy resolution in order to 
efficiently separate, e.g., ZZ and WW final states. Monte Carlo studies have shown that an ideal 
reconstruction algorithm [393], which finds each particle and measures its energy and direction 
with the detector resolution expected for single particles, could reach a jet energy resolution of 
14%/\/^- Over the years a jet energy resolution of SOYo/y/E has become accepted as a good 
compromise between the theoretically possible and practically achievable resolution. The RERECO 
package [389], developed at DESY, represents the first complete software implementation of the 
event reconstruction following the particle flow concept. A good performance of this implementation 
was demonstrated in a number of Monte Carlo studies [393,394]. The software includes pattern 
recognition in the vertex detector, TPC, and forward tracking system, followed by a track fitting 
procedure. The calorimeter clustering, track-cluster association and particle identification is done 
with the SNARK [395] package. The output is produced in the form of the reconstructed particle 
flow objects stored in the LCIO file. Despite the main advantage of being a complete and consistent 
implementation of the particle flow algorithm, the RERECO package has serious shortcomings. It 
is heavily bound to the TESLA TDR geometry implemented in Brahms and, therefore, cannot be 
used for detector optimization. Nevertheless, the revision of the Higgs analyses in the key channels 
can be done with this package for the TESLA TDR geometry. 

The MarlinReco package [396], meant as a successor of RERECO, is currently being developed 
at DESY. Presently, this package includes track reconstruction and fltting in the TPC, calorimeter 
clustering, track-cluster matching and particle identiflcation. Unlike RERECO, MarlinReco is not 
bound to any specific detector geometry and can be used for the detector optimization studies. 
Other reconstruction packages presented at the Snowmass workshop include org.lcsim [397] de- 
signed primarily for the SiD detector and the package designed for the GLD detector [398]. Both 
packages, MarlinReco and org. Icsim, use the same LCIO data model, which serves as an interface 
between the detector simulation and reconstruction software. The LCIO produced by full simula- 
tion contains a collection of hits recorded in various subdetectors. This information is used as an 
input for the reconstruction software. The reconstructed objects, namely tracks, calorimeter clus- 
ters and particle flow objects, are also stored in the LCIO flle and can be used for further analysis. 
Currently, the event reconstruction package used for the GLD detector stores information on the 
reconstructed objects in a ROOT file, but efforts are ongoing to make this package also compatible 
with the LCIO data format. 
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Existing reconstruction packages provide an option to disentangle detector effects from the 
inefficiencies of tlie reconstruction algorithm. This is done by emulating perfect pattern recognition 
in the tracking system and calorimeters and constructing true Monte Carlo tracks and calorimeter 
clusters. This approach, referred to as "perfect particle flow", is very important at the stage of 
detector optimization since it allows to trace and eliminate possible systematic effects coming solely 
from the reconstruction algorithms. 

8.2.3 Benchmarking Detector Performance with Measurements in the Higgs 
Sector 

One of the main issues that has to be addressed by future studies is the relation between detec- 
tor design and precision of measurements in the Higgs sector, going beyond the existing studies 
summarized in Sect. 18.11 We propose the following list of questions to be addressed by dedicated 
studies along with Higgs channels which can be used to benchmark the performance of the different 
detector components. 

• Vertex Detector 

How does the precision of Higgs hadronic branching ratio measurements depend on the specific 
configuration of the micro-vertex detector (number of layers, position of the innermost layer, 
etc.)? How is the analysis performance in other channels, e.g. HA bbbb, influenced by the 
micro- vertex detector configuration? The first study addressing these questions is presented 
in Ref. [385] . Efficient determination of the charge of secondary vertices is important for the 
reduction of combinatorial background in the HHZ channel with the H ^ bb decay and 
therefore crucial for the measurement of the Higgs self-coupling. How does the efficiency of 
the reconstruction of the vertex charge depend on the configuration of the vertex detector? 

• Tracking System 

What charged particle momentum resolutions can be achieved with the currently considered 
options for the main tracking device (TPC, silicon tracker)? How does the geometry of 
the tracking system (radius and length of TPC, number and position of layers in the silicon 
tracker, etc.) affect the charged particle momentum resolution? What impact will the charged 
particle momentum resolution have on the performance of analyses which utilize the inclusive 
HZ Xl^i" channel? How will the pattern recognition capabilities provided by different 
configurations of the tracking system influence the reconstruction of multi-jet final states with 
high track multiplicity, e.g. ZHH — > 6 jets or HZ WW {Z Z)qq — > 6 jets? 

• Calorimeters 

How does the precision on the Higgs photonic branching ratio measurement depend on the 
energy and angular resolution of the electromagnetic calorimeter? What impact will a specific 
option for electromagnetic and hadronic calorimeters (absorber material, active media, sam- 
pling fraction, transverse granularity) have on the efficiency of the reconstruction of multi-jet 
final states with high particle multiplicity in various Higgs channels? What jet energy resolu- 
tion can be achieved with a given option for the calorimeter system? How does the efficiency 
of r lepton identification via its hadronic decays depend on the configuration of the calorime- 
ter system. What effects will it have on the identification oi H ^ tt decays, which is crucial 
for the determination of the CP quantum numbers of Higgs bosons? 

• Muon System 
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What improvement can be achieved in the identification of high energy the muons if tracking 
and calorimeter information is complemented by the muon system? Should the muon system 
be realized in a conventional way by instrumenting the magnet return yoke or are alternative 
ways of implementing the muon detector more favorable? What active element (Resistive 
Plate Chambers, scintillating strips/tiles) is preferable for the muon system? These studies 
should be based on the analysis of final states involving high energy muons (HZ — > X^'^ iJi~ , 
e+e^ Huv ^i^vi)). 

• Overall Detector Performance and Performance of Reconstruction Software 

The overall detector and reconstruction software performance can be expressed in terms of 
the jet energy resolution. Hence, it would be desirable to investigate the dependence of jet 
energy resolution on a given detector configuration and also to study the impact of the jet 
energy resolution on the precision of measurements in the various Higgs channels involving 
multi-jet final states. 

The ultimate goal of these studies would be to establish a mapping between performance/configurations 
of different subdetectors and the precision of measurements in the Higgs sector. 

8.3 Proposed Strategy 

The detector optimization for the Higgs channels can be done in two steps. The first step is based 
on the simple MC smearing of the 4-momentum of isolated particles, jet energies and angles, and 
a parameterization of the heavy flavor tagging procedure. The resolution functions are extracted 
separately for each type of particle by running a full simulation program on single particle sam- 
ples. For charged particles it implies the study of the performance of the tracking system and the 
parameterization of the momentum resolution and reconstruction efficiency as a function of parti- 
cle momentum and polar angle 9. For neutral particles, detailed investigation of the calorimeter 
response has to be performed in order to quantify the energy, position and angular resolutions sep- 
arately for photons and neutral hadrons. The Higgs analyses relying on flavor-tagging will need as 
an input the parameterization of the vertex detector performance in terms of 6- and c- tagging effi- 
ciency versus purity. Realistic reconstruction packages can be employed to evaluate the jet energy 
resolutions for the case of perfect particle flow as well as for realistic particle flow algorithms. Once 
the main characteristics quantifying the detector response are extracted from the full simulation 
and reconstruction tools, fast Monte Carlo smearing can be applied to the generated background 
and signal samples relevant for a given Higgs channel. The program is supposed to emulate the 
measurements of the four- momenta of isolated particles and jets, heavy-flavor tagging, etc., using 
as an input the detector resolution functions extracted from the Monte Carlo studies with the full 
simulation. Results of the relevant analyses will directly reflect the potential of a specific detector 
model for measurements in the Higgs sector. 

On the other hand, with this simple parametric approach, one can study the precision of mea- 
surements as a function of the main detector performance characteristics, such as jet energy reso- 
lution, momentum resolution of the tracking system, etc., without referring to a specific detector 
model. The main advantage of the parametric approach is that it provides a fast procedure for map- 
ping the detector performance (described in terms of generic resolution functions and inefficiencies) 
into the precision of measurements in the key physics channels. 

However, this strategy has a major shortcoming. Lacking detailed simulation of the detector 
response and realistic event reconstruction, it fails to predict in an accurate way the analysis 
performance in the physics channels where details of the pattern recognition in the detector play 
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a crucial role. This situation can be illustrated by two examples. Successful determination of 
Higgs CP quantum numbers in the H — > t~^t~ channel will require efficient pattern recognition in 
the electromagnetic calorimeter in order to identify the tau decay mode and to measure precisely 
the decay products of the tau leptons. Accurate analysis in this channel is possible only with 
the full detector simulation and realistic reconstruction tools. The second example is given by 
measurements of Higgs self-couplings in the double Higgs-strahlung process, e^e~ — > HHZ. The 
large combinatorial background in this channel can be reduced by identifying the vertex charge, 
which would allow to discriminate between b and b jets. Again, reliable prediction of the vertex 
charge identification efficiency requires a detailed study of the vertex detector performance with 
the full simulation and realistic reconstruction tools. 

At this point one can conclude that it is reasonable to employ fast parametric simulation at 
the early stage of the detector optimization to reject those detector configurations that result in 
unacceptable deterioration of the precision of measurements in the Higgs sector. 

For detector configurations passing the first step of optimization with the fast parametric simu- 
lation, a more detailed study is needed that employs the full detector simulation and reconstruction. 
The large-scale analyses including background estimates represent a time-consuming procedure in 
this case. Therefore, we propose first to run full simulation programs, implementing a chosen detec- 
tor geometry, only on HEPEVT or StdHep files of signal samples. Performing event reconstruction, 
selection and analysis of only the signal samples, one can get a rough idea about the dependence 
of the analysis performance on a specific detector configuration. Variations in analyzing power of 
variables crucial for a given Higgs channel will serve as a major criterion to discriminate between 
different detector designs. An example of such an approach is illustrated in Fig. 18.81 It presents 
reconstructed Higgs boson mass spectra in the e~^e~ Hvv bbvu channel after performing 
the full simulation of the TESLA TDR detector with Mokka. The figure compares results of the 
realistic reconstruction based on the MarlinReco package and the perfect particle fiow algorithm. 
It is very important that the results of the realistic event reconstruction are cross checked with 
the simulation of perfect reconstruction. This is necessary to disentangle detector effects from the 
inefficiencies in pattern recognition algorithms and to guide the efforts of the developers of the 
reconstruction software. In this particular example, the obtained resolution on the reconstructed 
Higgs boson mass plays the role of a parameter evaluating the performance of the given detector 
model and event reconstruction package. Similar studies can be done in other channels for various 
detector options which may emerge in the future as a result of the detector R&D program for the 
ILC. 



Figure 8.8: The distribu- 
tion of the dijet mass in 
the e'^e" — > Hi/v — > bbvu 
channel. Distribution is ob- 
tained after performing full 
simulation of the TESLA 
TDR detector with Mokka. 
The result of the realistic 
event reconstruction with the 
MarlinReco package (dots) is 
compared with perfect recon- 
struction (line). 



I e^e- Hvv bbvv \ls=800Gev1 




m„(GeV) 
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